Navier–Stokes Equation and its Fractional Approximations

We consider the Navier–Stokes equation (N-S) in dimensions two and three as limits of the fractional approximations. In 2-D the N-S problem is critical with respect to the standard $$L^2$$L2 a priori estimates and we consider its regular approximations with the fractional power operator $$(-P\Delta )^{1+\alpha }$$(-PΔ)1+α, $$\alpha >0$$α>0 small, where P is the projector on the space of divergence-free functions. In 3-D different properties of the N-S problem with respect to the standard $$L^2$$L2 a priori estimate are obtained and the 3-D regular approximating problem involves fractional power operator $$(-P\Delta )^s$$(-PΔ)s with $$s>\frac{5}{4}$$s>54. Using Dan Henry’s semigroup approach and the Giga-Miyakawa estimates we construct regular solutions to such approximations. The solutions are global in time, unique, smooth and regularized through the equation in time. Solution to 2-D and 3-D N-S equations are obtained next as a limit of such regular solutions of the approximations. Moreover, since the nonlinearity of the N-S equation is of quadratic type, the solutions corresponding to small initial data and small f are shown to be global in time and regular.


Introduction
Fifteen years ago, in the monograph [3], we were studying a direct generalization of semilinear parabolic equations, namely abstract semilinear equation with sectorial positive operator in the main part. Following the idea of Dan Henry [21], such equation will be studied using classical techniques of the theory of ordinary differential equations, modified however to cover equations with unbounded operator in a Banach space. Such approach, located inside the semigroup theory, proves its strength and utility in the study of several classical problems; some of them were reported in [3,21]. But such technique offers further possible generalizations, first to study the problems, like e.g. Korteweg-de Vries equation and its extensions [7,8,12,13], where the solutions are obtained as a limit of solutions to parabolic regularizations of such equations (the method known as vanishing viscosity technique, originated by Hopf, Oleinik, Lax in 1950th); see also [9,10]. Another possible application of Henry's technique is to study critical problems (e.g. [3,41]), not falling directly into the class of semilinear sectorial equations because the nonlinear term in it is of the same order, or value, as the main part operator. A recent paper [11] was devoted to such type problem, the quasi-geostrophic equation in R 2 (e.g. [4][5][6]26,44,45]). But certainly the most celebrated example of such critical problem is the Navier-Stokes equation in dimension two. We will obtain and study its solutions constructed as limits of solutions to sub-critical approximations (1.2) when α → 0 + where, in 2-D, the P operator is replaced with its fractional power −(−P ) 1+α , α > 0 (P is the projector on the space of divergence-free functions; see e.g. [15]). In 3-D a higher order diffusion term like −(−P ) s , s > 5 4 , is considered in the corresponding regularization of the N-S equation because the properties of the problem with respect to the natural L 2 estimate in that case are different.
The classical 3-D Navier-Stokes equation (N-S) considered here has the form: where ν > 0 is the viscosity coefficient, u = (u 1 (t, x), u 2 (t, x), u 3 (t, x)) denotes velocity, p = p(t, x) pressure, and f = ( f 1 (x), f 2 (x), f 3 (x)) external force, and is a bounded domain with C 2 boundary. It is impossible to recall even the most important results devoted to that problem, since the corresponding literature is too large; see anyway [2,[14][15][16][17]20,23,24,[29][30][31][36][37][38]40] together with the references cited there. However, we will try to recall the most important steeps in the semigroup approach to the Navier-Stokes problem, the approach used in the present paper. In 1962 and 1964, thirty years after the Leray fundamental paper [31], the semigroup treatment of the 2-D and 3-D N-S was stated by Fujita and Kato [16,25]. Existence and uniqueness of strong solutions was obtained by solving the evolution equation in X 2 space, together with basic semigroup formulation and tools. That idea was extended later by Giga and Miyakawa in a series of papers including [17,18,20] in early 1980th. Detail informations concerning analyticity of the semigroup generated by the Stokes operator with Dirichlet boundary condition on general L r type spaces, deep estimates of the nonlinearity (including (1.15), below), existence-uniqueness theory in fractional order spaces D(A γ ) were reported in that references. Here also the papers [23,43] should be mentioned, contributing to the latter results. Finally, the 1985 monograph by von Wahl, collects and extends in the direction of the global in time solvability and regularity, the existing semigroup results concerning the N-S equation. Meanwhile, in 1981 (with first version in 1974), an abstract approach to semilinear equations with sectorial operators was presented by Dan Henry in [21]. However, the N-S equation is only mentioned there. Finally, an extended discussion of the semigroup approach to the Navier-Stokes equation can be found in the review article [19].
Let describe briefly the contents of the present paper. In space dimensions 2 and 3 the N-S equation possess local in time regular solutions. Using strongly the technique of [20], we recall local in time solvability in Theorems 1.5 and 1.6. We analyze further criticality of the Navier-Stokes equation (compare [3,40,41]), in a sense that available for it L 2 ( ) a priori estimate (1.11) is not sufficient to control its nonlinearity through the viscosity term ν P u. Consequently, possible is a balance between the income from nonlinearity and the stabilizing action of the viscosity so that the local solutions can not be extended globally in time. For small initial data the decisive role is played by viscosity, while for larger initial data the nonlinear term is strong enough to unable us to obtain regularization of the solutions through the main part operator. We will try to see such effect through the estimates obtained in the paper.
In 2-D, our idea is to improving a bit the viscosity term to make the whole problem sub-critical, such that the improved viscosity together with the known L 2 ( ) a priori estimate will control the nonlinear term. The way for obtaining such effect is to replace the classical viscosity term Au = −ν P u through a bit higher fractional diffusion −A 1+α u, with small α > 0. Next, we will study the process of letting α to 0 + ; which properties/estimates of the solutions of regularized problems (with α > 0) are lost in such a limit. We also look at 3-D N-S equation in which case a stronger diffusion term, with −P operator in the power s > 5 4 , will be considered to guarantee together with the standard L 2 ( ) a priori estimate the control on the nonlinear term.
In 2-D, instead of (1.1), we consider a family of sub-critical problems, with α ∈ (0, 1 2 ]: (1. 2) The approximation proposed in 3-D is given in (2.6) with s > 5 4 . Until very recently, equations with fractional power operators were not seriously studied in the literature, therefore we hope the regularizations proposed will help to understand the difficulties faced in the original N-S equation. An analogous phenomena was studied recently [11] for the quasi-geostrophic equation in R 2 (e.g. [4][5][6]44,45]). The technique used in the present paper is similar to that of [11].
There were several tries of replacing the classical N-S equation, or the viscosity term in it, with another equation having better properties of solutions, starting with Leray α -regularization reported in paper [31], see also [14]. Modification of one factor in nonlinearity, using mollifier, was sufficient to improve properties of solutions. Another modification of the N-S equation was proposed by Lions in [33, Chapter 1, Remarque 6.11] (see also [30]), where the − operator was replaced with − + κ(− ) l , l ≥ 5 4 . Further modifications can be easily found in the literature, in particular in the closing part of [29]. In fact, in this paper we follow the idea of Lions to replace the diffusion term with a stronger one fractal diffusion term.

Introductory Facts
Notation We are using standard notation for Sobolev spaces. Compare [39] or [3, Chapter 1] for properties of fractional order Sobolev spaces; see also [22] for Sobolev type embeddings. For r ∈ R, let r − denotes a number strictly less than r but close to it. Similarly, r + > r and r + close to r . When needed for clarity of the presentation, we mark the dependence of the solution u of (1.2) on α ∈ (0, 1 2 ], calling it u α , or u for the approximation in 3-D. Throughout the text the 'standard' or 'natural' a priori estimates mean the estimates (1.11) and (1.12).
Recall first [3,Chapter 3], that studying in a Banach space X an abstract Cauchy's problem with sectorial positive operator A and solutions varying in the phase space knowing an a priori estimate of all its potential X β solutions; u(t) Y ≤ const in another Banach space D(A) ⊂ Y , we say that the nonlinear term F is sub-critical relative to such a priori bound, if for each such X β solution u(t) an estimate is valid for all t ∈ (0, τ u 0 ), where τ u 0 is the 'life time' of that solution. In that case (e.g. [3,Chapter 3]) the X β norm of the local solution will be bounded on [0, τ u 0 ), which allows to extend such solution globally in time. Similar extendibility property holds even if θ = 1 provided that β < 1 (as we assume above), which in turn is unknown if θ = β = 1; the nonlinearity F being in the latter case called critical relative to that a priori estimate. Note that critical nonlinearities are 'of the same order' in the equation as the main part operator A (compare [40,41]). The main part operator A will not control the nonlinearity in that case, unless we find a better a priori estimate.

Properties of the Stokes Operator
Familiar in the theory of the N-S equation are the following spaces: We define also the Stokes operator (N = 3) where P r denotes the projection from L r ( ) to X r given by the decomposition of L r ( ) onto the space of divergence-free vector fields and scalar-function gradient (e.g. [37] generates on X r an analytic semigroup {e −t A r } for arbitrary 1 < r < ∞. A complete description of the domains of fractional powers of the Stokes operator A r = −ν P r , D(A r ) = D(− ) ∩ X r , can be found in [20, p. 269], or in [18]. Note further, that the domains of negative powers of the operator A r are introduced through the relation (e.g. [20, p. 269]): Thus, D(A β r ), β < 0, is the completion of X r under the norm A β r · 0,r . It is also easy to see, that the resolvent of the operator A 2 fulfills an estimate: (σ (A 2 )) ≥ νλ 1 , where λ 1 is the first positive eigenvalue of − in L 2 ( ) considered with the Dirichlet boundary condition. The same estimate remains valid for the operator considered in L r ( ) with any r ∈ (1, ∞). It follows further from [20,Lem. 3.1], that the resolvent of A r is compact, also the embeddings D(A β r ) ⊂ D(A α r ) are compact when 0 < α < β ([21, Th. 1.4.8]). In fact the operator A = A 2 (we skip the subscript further for simplicity) is self-adjoint in the Hilbert space L 2 ( ); see e.g. [17,20]. For such type operators the powers of the order (1 + α) have similar properties; in particular they are also sectorial operators.
Consequently, the operators A r , 1 < r < ∞, are sectorial positive. Fractional powers of the order 1 + α for such operators are introduced through the Balakrishnan formula ( [27,34]): We recall next the, specific for the N-S equation, standard a priori estimate. Calculating on regular solutions (which extends to weak solutions), it is obtained multiplying (1.1) in [L 2 ( )] 3 through u, to get: since the nonlinear component vanish in that calculation due to condition div u = 0: The term ∇ p · udx is transformed as follows: Consequently, for divergence-free f , an L 2 ( ) estimate is obtained: (1.10) thanks to the Poincaré inequality. Integrating the above, we obtain: where c P denotes the constant in the Poincaré inequality. Having already the last estimate one can return to (1.10), to see that for arbitrary T > 0. These are the strongest natural a priori estimates known for solutions of the N-S equation.

Local in Time Solvability of the 3-D and 2-D N-S Problems
We will rewrite N-S equation in a form of an abstract parabolic equation with sectorial positive operator and solve it using Dan Henry's approach ( [3,21]). I want to make a comment here, that using such approach we have fairly large choice of phase spaces (= space where solution varies). In fact, calling here paper [35], as far as we consider the corresponding to (1.3) linear Cauchy problem with sectorial positive operator: we can 'set it' at any level of the fractional power scale X β = D((−A) β ), β ∈ R, corresponding to −A (see [1,Section V.2] for an extension of that idea). When we move to the semilinear problem (1.3) with nonlinearity F subordinated to A, it is an art to choose the proper level at that scale to be the phase space for semilinear problem. For that, we need to consider a priori estimates available for the specific equation, usually of physical origin, e.g. following from energy decay or conservation of mass valid in the process described through the equation. The full semilinear problem will be next written abstractly as The standard way to set the problem (1.1) in the above setting, in L 2 ( ) (see e.g. [17,20,37]), is to apply to the equation the projector P = P 2 : The pressure term disappears then from the equation. The realization A of the Stokes operator acts from D(A) → H . We also introduce the energy space V = {u ∈ [H 1 0 ( )] 3 ; divu = 0}, and the simplified notation for the nonlinearity: F(u) = −P(u · ∇)u.
We will rewrite the classical N-S equation in an equivalent form, using the property of the divergence-free functions. We have: We need to emphasize here, that all the technical tools needed for local in time solvability of the N-S equation (in standard spaces), like sectoriality of the Stokes operator, description of its domain and needed estimates, were reported in [17,18,20]. We recall next an estimate, important for further calculations, borrowed from [20, Lemma 2.1]. A similar observation was given also in [21, p. 18] in dimension one, see also the review of earlier results in [20, p. 271].
Consequently, the following estimate holds: (1.14) Observation 1. 4 We have also the following estimates used further in the text. From [20, Lemma 2.1] we get, for all N ∈ N: Now, for any δ ∈ (0, 1 2 ), using the theory of interpolation; In a similar way, starting from the estimates: for any δ ∈ (0, 1 2 ), we get: for the space dimension N = 2, and in the base space X − 1 8 for the space dimension N = 3. The corresponding phase spaces will be; X 1 2 [20,Proposition 1.4]). Note that, in both cases, the phase spaces are contained in the space [L ∞ ( )] N . Note also, there is another possible choice of the phase spaces (e.g. [20]), if we decide to work in the spaces [L r ( )] N , N = 2, 3, with r > N .
We will formulate now the corresponding local existence results for N = 2, 3. Case N = 3. The main tool is the estimate taken from [20, Lemma 2.2] (with δ = 1 8 , θ = ρ = 3 4 ): Since the form above is bi-linear, we have also the following consequences of the last estimate: Consequently, the nonlinear term F(u) + P f = −P(u · ∇)u + P f acts from According to abstract theory in [3,21] and thanks to the properties of the Stokes operator as reported in [17,20], this suffices to obtain a local in time solution of the 3-D equation (1.1) (also of 3-D equation (2.6)), more precisely we have: Here τ > 0 is the 'life time' of that local in time solution. Moreover, the Cauchy formula is satisfied: where e −At denotes the linear semigroup corresponding to the operator A. We need also to mention that the considered here mild solutions have additional regularity properties, as described in particular in [3, p. 218 ). This property is used in the calculations below.
Case N = 2. We will use a version of the estimate in [20, Lemma 2.2] (with δ = 1 4 , θ = ρ = 1 2 ): Since the form above is bi-linear, we have also the following consequences of the last estimate: Consequently, the nonlinear term F(u) + P f = −P(u · ∇)u + P f acts from ). According to [3,21], this suffices to obtain a local in time solution of the 2-D equation (1.1) (also of 2-D equation (1.2)), more precisely: Here τ > 0 is the 'life time' of that local in time solution. Moreover, the Cauchy formula is satisfied: where e −At denotes the linear semigroup corresponding to the operator A. Remark 1.7 Another choices of the pairs base/phase spaces are possible [20]. In particular, setting for N = 2:

Properties of 2-D and 3-D N-S Equation with Respect to the Known a Priori Bound
We will consider now criticality of the N-S equation (1.1) in the sense stated e.g. in [3,40,41], that means 'the nonlinear term is of the same order as the main part operator with respect to the given a priori estimate'. To exhibit this concept better we will assume in the definition below that (1.3) is stated in a Banach 'base' space denoted Z , A is a sectorial positive operator in Z and F is a map Lipschitz continuous on bounded sets from Z β into Z where β ∈ [0, 1). In particular we will consider in this setting Z β solutions as in [3, p. 55].
Definition 1. 8 The problem (1.3) will be called sub-critical with respect to the given a priori estimate in Y , if for all its possible Z β solutions an estimate holds The problem (1.3) will be called critical with respect to the given a priori estimate in Y if, instead of (1.25), for all its possible Z β solutions an estimate holds for all t ∈ (0, τ u 0 ), the life time of the considered Z β solution.
Using the abstract setting as in Definition 1.8 we now indicate that the 3-D N-S can manifest a different behavior; namely, In case of the Navier-Stokes equation (1.1) the mentioned above a priori estimate will be the Y = [L 2 ( )] N estimate (1.11).
Following first (1.19), thus choosing X − 1 8 and X 3 4 + to play a role of Z and Z β respectively, we will see that (1.27) holds for (1.1) with respect to the standard L 2 estimate. Indeed, the estimate (1.18) written for the local solution u = u(t) obtained in Theorem 1.5 extends for N = 3, with the use of the Nirenberg-Gagliardo type estimates, to: which gives the exponent above equal 9 7 .
On the other hand, using the estimate of Corollary 1.3, we will find for the 3-D N-S equation that 29) in which case = 3 2 . Finally, for the further use, we will check how large the exponent s > 1 should be, for the 'strengthen diffusion' of the form A s replacing the usual operator A in (1.1) to make the L 2 ( ) estimate critical. Using again Corollary 1.3 we get that (1.30) We find such critical value of s from the condition; 3 2(2s−1) = 1. Consequently, s = Case N = 2. Using Corollary 1.3 it is easy to see criticality of the 2-D Navier-Stokes equation. More precisely, to verify that its nonlinearity is critical with respect to the standard L 2 ( ) a priori estimates, as a map from X 31) and no better estimate (with exponent smaller than 1) seems possible.

32)
where the Nirenberg-Gagliardo type estimates were used: (1.33) No better estimates seems possible.
Note further, that in 2-D the nonlinearity is also critical with respect to the L 2 ( ) a priori estimate as a map between X 1 = D(A) ⊂ [H 2 ( )] 2 and X . Indeed, the following estimate holds, (1.34) and no smaller exponent on the H 2 ( ) norm seems possible.

Global in Time Solutions in 3-D: Small Data
As well known (e.g. [41], [ To get estimate of the solution u in [H 1 ( )] 3 multiply (1.2) by Au to obtain: with the standard use of Young's inequality, and the embedding 3 , we arrive at the differential inequality (e.g. [42]): Analyzing its right hand side, real function 3 , we see that g(0) = C ν P f 2 [L 2 ( )] 3 > 0, g (0) < 0, and g has a minimum for the argument z min = 3c ν 14C ν 3 4 , with g(z min ) < 0 when the 'free term' C ν P f 2 [L 2 ( )] 3 is small. More precisely, to keep the value of y(t) bounded for all positive times, we need to assume the smallness hypothesis: Let the data: u(0) 2 [H 1 ( )] 3 and P f 2 [L 2 ( )] 3 be so small, that: Note that ν

Regularization of the 3-D N-S Equation
It was shown above that in the 3-D case the N-S equation (1.1) may behave with respect to the standard estimate according to (1.27). Global in time extendibility of the local solutions constructed in Theorem 1.5 is in general unknown (unless for small data). The viscosity term in the equation (1.1) together with the standard a priori estimate are not strong enough to control the nonlinearity. In the present subsection we discuss an approximation/regularization (with stronger diffusion term) of the 3-D N-S equation, for which the L 2 ( ) estimates (same as for the original equation) are sufficient to make such problems sub-critical. The idea of such regularization was proposed first by Lions in [33, Chapter 1, Remarque 6.11].
We consider approximation of the original 3-D Navier-Stokes equation having global in time, unique and regular solutions. Consider namely the approximation/regularization of (1.1) having the form: It is clear how to determine the proper, sufficiently large, value of the exponent s > 1 to guarantee, together with the L 2 ( ) estimates valid also for solutions of (2.6), that the nonlinear term is subordinated to the main part operator −(A + A s ). We need to compare the bound obtained for the nonlinearity with the income from the improved viscosity term. Estimate of the nonlinearity, obtained from (1.16), reads Let s > 5 4 . With the use of the Nirenberg-Gagliardo type estimate we obtain the subordination condition where θ = 2 2s− 1 2 < 1. Consequently, the standard L 2 ( ) a priori estimate is sufficient to assure the global in time extendibility of such local solutions u s when s > 5 4 . We will study such approximation next. 5 4 proposed for regularization of (1.1) in [33, Chapter 1, Remarque 6.11] is the same (also, as in (1.30)).

3-D N-S Equation
Let us study the introduced above regularization (2.6) of the original 3-D N-S equation with exponent s > 5 4 . We formulate first the corresponding existence result. Here τ > 0 is the 'life time' of that local in time solution. Also, the corresponding Cauchy formula is satisfied. Furthermore, together with the standard [L 2 ( )] 3 a priori estimate valid for solutions of (2.6) uniformly in > 0, the obtained above local solution will be extended globally in time in the class (2.9).
The proof, similar as for Theorem 1.5, is omitted.
Passing to the limit → 0 + in (2.10) We will describe next shortly the process of passing to the limit, as → 0 + , in such approximations. The idea is similar as in [7,8,33], using the parabolic regularization technique. We write the corresponding approximating equation: with s > 5 4 (but close to 5 4 ), and denote its solution by u . Applying to it the operator A −s+ 1 2 we get Using to nonlinear term the estimate of [20, Lemma 2.2] with δ = s − 1 2 , = s − 1, we get where 1 z = 1 2 + 2 3 , consequently z = 6 3+4(s−1) < 3 2 (but close). Due to the standard L 2 estimate, all the right hand side components in (2.11) belong to L 2 (0, T ; [L 2 ( )] 3 ) (in particular − A 1 2 u; note that 7 2 − 2s < 1 in (2.12)). Consequently, Thus the Lions-Aubin compactness lemma (e.g. [33,Chapter I,5]) will be used to conclude that any sequence {u n } has a subsequence convergent almost everywhere in (0, T ) × , as n → 0 + . Also, any sequence {u n } has a subsequence convergent in 3 ) of u , over a sequence n → 0 + . Remember that from the standard a priori estimate written for the problem (2.6) it follows that the D(A with const independent on > 0. We recall a weak formulation of the approximating equations (2.10): (2.14) where v ∈ D(A 1 2 ) arbitrary. Note that, due to (2.13), the second right hand side component vanishes when → 0 + . For the nonlinear component we have (2.15) due to the estimate Finally, thanks to [37, Lemma 1.1] and the convergence u n → U in L 2 (0, T ; (2.17) in the sense of 'scalar distributions' (e.g. [33,37]); the derivative d dt in the sense of distributions.
Consequently, we were able to pass to the limit in the weakly formulated regularization (2.10) of the 3-D N-S equation. We obtain where v ∈ D(A 1 2 ) was arbitrary, and the convergence of the first right hand side component is in the weak sense. We obtain a weak solution to the 3-D N-S equation in the spirit of the Leray original definition (e.g. [40, p. 139]). Such solution is global in time, as was first observed in [24].
Note that the weak solution obtained above is global in time, while eventually not unique, since it depends on the chosen subsequence. For regular initial data it must coincide (for small times) with the unique local strong solution described in Theorem 1.5, since the last exists on certain time interval t ∈ [0, τ ) and fulfills (2.18) for T < τ.

2-D Critical N-S Equation as a Limit of Sub-Critical Approximations
We will describe now the convergence of the solutions of the 2-D fractal approximations (1.2) to the solution of the limiting 2-D N-S equation. Precisely as in Theorem 1.6, the solutions u α of (1.2) (the superscript is added for clarity) will be constructed in the class: Further, the L 2 ( ) a priori estimates are satisfied for u α uniformly in α ∈ (0, 1 2 ]. Consequently, we claim: Theorem 3.1 The local in time solution u α (t) of the approximating problem (1.2) constructed as in Theorem 1.6 will be extended globally in time in the above class. Moreover, the standard L 2 ( ) a priori estimate (1.11) is valid for it uniformly in α ∈ (0, 1 2 ]. All the nice properties of the abstract semilinear sectorial equation, like elegant theory of existence and uniqueness, regularization of the solution for positive times, are valid for solutions u α (t) with any α > 0. Typical regularization goes from the phase space X β , β < 1, to the space X 1 = D(A); see e.g. [3,21]. However, for problems (1.2), estimates of norms better that L 2 ( ) will depend on α > 0, possibly blowing up when α → 0 + , and therefore can not be extended in general to a limit solution u.

2-D (1.1) as a Limit of (1.2) When α → 0 +
A precise description of letting α → 0 + in the equation (1.2) is given next. In this section we consider the solutions u α of (1. Solutions u α fulfill also, uniformly in α ∈ (0, 1 2 ], estimate (1.11) in [L 2 ( )] 2 . More precisely, for such solutions of (1.2) we have an estimate: This is the main information allowing us to let α → 0 + in the equation (1.2). Passing in (1.2) to the limit α → 0 + . We look at (1.2) as an equation in [L 2 ( )] 2 , and 'multiply' by the test function We will discuss now the convergence of the terms in (3.3) one by one. Note that when α → 0 + then, by Lemma 4.1, , we obtain: where u is the weak limit of u α in [L 2 ( )] 2 as α → 0 + (over a sequence {α n } convergent to 0 + ; various sequences may lead to various weak limits). We return to (3.3) to see that letting α → 0 + over a sequence {α n }, where u denotes weak limit in [L 2 ( )] 2 of such sequence, we have since the right hand side is convergent. Consequently, the left hand side has a limit as Note that, for the 'test functions' A −1 φ varying in a separable Banach space, passing countable many times to a subsequence, we can chose a common subsequence proper for all test functions in dense subset of the space. Consequently, the equation below will be fulfilled in the whole space. We obtain which is a weak form of the limiting equation.

Separation of terms
The two terms of [A −α u α t + P(u α · ∇)u α ] will be separated when letting α → 0 + . More precisely we have:

Remark 3.2 Since the approximating solutions u α satisfy (in particular)
the time derivative d dt and the convergence are understood in D (0, T ; [L 2 ( )] 2 ) (space of the 'scalar distributions' [33]). Consequently, Convergence of the nonlinear term F(u α ) = P(u α ·∇)u α will be discussed next. As seen from (1.12), the approximating solutions u α are bounded in L 2 (0, T ; [H 1 0 ( )] 2 ) uniformly in α > 0. Estimate (1.15): together with the consequence of Eq. (1.2) show that u α t are bounded in L 2 (0, T ; D(A −1 )) uniformly in α ∈ (0, 1 2 ]. By Lions-Aubin compactness lemma [33,Theorem 5.1], the family {u α } α∈(0, 1 2 ] as bounded in the space . In particular any sequence {u α n }, α n → 0 + , has a subsequence convergent almost everywhere in (0, T ) × . Moreover, the above compactness allows to pass to a limit (using estimates of Observation 1.4) in the nonlinear term in (3.3). Indeed, if u α n → u, α n → 0 + in the above sense then, using a consequence of [20,Lemma 2.2]; for small > 0 The construction presented above allows us to formulate the following theorem: (3.13) Due to denseness of the set D(A 3 4 ) in X 2 , the right hand side of (3.13) defines a unique element in X 2 . The left hand side ω φ is defined in (3.6) and discussed in Remark 3.2. 2 )) is unique. Indeed, if u 1 , u 2 are two such (weak) solutions, applying projector P, taking the difference of the equations and multiplying the result by w = u 1 − u 2 , we get: (3.14) The last term vanishes for divergence-free functions. The earlier term, using an equivalent form of the nonlinearity in (1.13) and the Nirenberg-Gagliardo inequality, is estimated as follows: Inserting the last estimate into (3.14), using Cauchy's inequality, we obtain Since , due to the Gronwall lemma. (3.16) Further, due to (1.16), Nirenberg-Gagliardo and Young inequalities, For small > 0, α > 0, we thus have 1+5α < 4. It follows from the natural a priori estimates and Gronwall inequality that the family {u α } α∈(0, 1 2 ] is bounded in the norm Consequently, it is precompact in with arbitrary p > 1. The convergence u α → u as α → 0 + is thus verified in the above space.

Remark 3.6
To understand better the method used in the study of global extendibility of solutions we will discuss shortly an example 1 of the Burgers type system in 3-D, obtained by neglecting the viscosity in homogeneous 3-D N-S equation: It is easy to see that each component of the sufficiently regular (that is, varying in [L ∞ ( )] 3 ) solution of (3.20) fulfills Maximum Principle: We are thus given a natural a priori estimate in Y = [L ∞ ( )] 3 for such system. The nonlinear term is as in the N-S equation, and we have the estimates: which shows the nonlinearity is sub-critical in that case. Consequently, the local solution will be extended globally in time.
The above example shows that our approach is sensitive (through subordination condition (1.8)) not only on the form of nonlinearity, but also another specific properties of the problem. It also indicates the role of the pressure in the classical N-S equation, the term which seems responsible for the delicate properties of solutions of that equation.

Remark 3.7
It is evident from the considerations above that the phenomenon of loosing regularity by local smooth solutions of 3-D N-S is possible only if they enter the supercritical range of the nonlinearity (compare (1.28)). It would be therefore interesting to consider local solutions corresponding to initial data u 0 fulfilling where 0 < ε < 2 7 , as eventual candidates for such phenomenon.

Conclusion
The PROBLEM connected with the 3-D Navier-Stokes equation is not solved in the paper, while we are explaining its reason. The viscosity in the classical N-S equation together with the standard a priori estimate could not be used in general to prevent better norms of its solutions from blowing-up in a finite time, especially in 3-D case. In order, for small initial data we enjoy the standard property of problems with quadratic nonlinearity, that near zero function the estimates of nonlinearity obtained from the equation are proportional to the square of the norm, therefore there is a ball centered at zero the solutions originated in it will never leave that ball. For such solutions, the nonlinear term in N-S is subordinated to the main part operator. Of course, if a better estimate was known, for example the estimate in L N ( ) (L N + ( )), it would make the N-D Navier-Stokes equation critical (sub-critical) due to the following calculation Consequently, the left hand side tends to zero as 0 < β → 0 + .
Proof Our task is, for fixed φ ∈ X and β near 0 + , to estimate the expression: In the estimates we are using the following properties; taken from [34, p. 62] equality valid for η ∈ (0, 1) the simple formula: Returning to the proof, we split the integral in (4.2) into (0, L) and (L , ∞) and estimate the first part, where L > 0 will be chosen later. Note that letting β → 0 + the result of the estimate above is bounded by | sin(πβ) π |2L(1 + M) φ X → 0, for any fixed L > 0.
Next using (4.3), the integral over (L , ∞) is, for φ ∈ X , estimated as follows: where due to (4.4) we see that > 0 arbitrary fixed. Consequently we obtain: for sufficiently large value of L ≥ 1, as specified in (4.7). Note that letting β → 0 + in the resulting estimate we have: sin(πβ) πβ L −β + sin(πβ) π L −β−1 β + 1 M φ X ≤ + sin(πβ)L −1 M φ X , (4.9) for chosen large value of L. For such L we get a final estimate of the integral in (4.2) having the form: where > 0 was arbitrary. The right hand side of (4.10) will be made small when we let β near 0 + , noting was an arbitrary positive number.

Properties of the Fractional Powers Operators
Recall first the Balakrishnan definition of fractional power of non-negative operator (e.g. [27, p. 299], [28]). Let A be a closed linear densely defined operator in a Banach space X , such that its resolvent set contains (−∞, 0) and the resolvent satisfies: Then, for η ∈ (0, 1), Note that there is another definition, through singular integrals, of the fractional powers of the (− ) −α in L p (R N ) frequently used in the literature. See [34,Chapter 2.2] for the proof of equivalence of the two definitions for 1 < p < N 2 α ; see also [11, section 4.3]. See further [32] for the commutator estimates.
Moment inequality We recall here, the moment inequality estimate valid for fractional powers of non-negative operators. They are suitable to compare various fractional powers. Recalling [46, p. 98], we have the following version of the moment inequality with precise constant; for 0 ≤ α < β < γ ≤ 1, where φ ∈ D(A γ ).
We are using the following bound, valid for positive operators, in the calculations below: (4.14) We thus have: (4.15) The second integral over (L , ∞) is estimated next: (4.16) Minimizing with respect to L > 0 we get L min = M A 2 φ X (M+1) Aφ X , which leads to the final estimate: as α → 0 + .