Until now, all the investigations on fractional or generalized Navier-Stokes equations have been done under some restrictions on the different values that can take the fractional order derivative parameter . In this paper, we analyze the existence and stability of nonsingular solutions to fractional Navier-Stokes equations of type () defined below. In the case where , we show that the stability of the (quadratic) convergence, when exploiting Newton’s method, can only be ensured when the first guess is sufficiently near the solution . We provide interesting well-posedness and existence results for the fractional model in two other cases, namely, when and

1. Introduction

Navier-Stokes equations (NSEs) are considered as the foundation of fluid mechanics and, strangely enough, are rarely recorded in their entirety. They describe motions of Newtonian fluid flows irrespective of whether they are laminar or turbulent. NSEs are usually obtained according to the following considerations [1].(i)Constitutive relations governing the fluid are isotropic, which means independent of the orientation of the coordinate system axes.(ii)At rest the fluid obeys the laws of statics.(iii)The equation of the fluid is objective; that is, tensors are used. It is independent of the Galilean reference frame in which it is expressed, and independent of the observer.Although existence and characterization of smooth or nonsingular solutions of NSEs have been comprehensively investigated in numerous works [27], they remain open problems. Fractional calculus applied to differential equations has captured a huge amount of attention recently [812] and it was shown that its applications lie in modeling many problems and phenomena in applied technology, engineering, and sciences including applied mathematics, computing, physics, biology, chemistry, economic, and other domains of applications.

2. A Brief Overview on Fractional Order Derivatives

There is a huge amount of works concerning various definitions of fractional derivatives, but two of them have emerged as the most popular and most applicable on dynamical systems: the Caputo derivative and Riemann-Liouville derivative. Caputo derivative is expressed as For the fractional derivative in the sense of Riemann-Liouville, we have the following definition: In applied sciences in general and dynamical system in particular, each of these fractional derivatives presents some advantages and disadvantages [8, 13]. For example, Caputo derivatives are defined only for differentiable functions while functions that have no first order derivative might have fractional derivatives of all orders less than one in the sense of Riemann-Liouville [11]. Another example is the fact that the Riemann-Liouville derivative of a constant is not zero unlike Caputo’s derivative which gives zero when applying to a constant function. However Caputo’s derivative demands higher conditions of regularity for differentiability. Then, to compute the fractional derivative of a function in the Caputo sense, we have to calculate first its derivative. Fortunately for Riemann-Liouville derivative, the author in [14] proposed the following alternative definition:the above modified Riemann-Liouville derivative seems to have advantages of both the standard Riemann-Liouville and Caputo fractional derivatives: it is defined for arbitrary continuous (nondifferentiable) functions and the fractional derivative of a constant is equal to zero. However from its definition we do not actually give a fractional derivative of a function say but the fractional derivative of can lead to fractional derivative that is not defined at the origin for some function for which does not exist. We note that Caputo and Riemann-Liouville may have their disadvantages but still remain the best definition of the fractional derivative according to their easy applicability. We need the following definition.

Definition 1 (partial derivatives of fractional order). Assume now that is a function of variables , also of class on . We define as partial derivative of order for with respect to the function if it exists, where is the usual partial derivative of integer order .

In this paper, we consider the fractional derivative in the sense of Caputo to obtain and analyze the generalized (fractional) NSE in 3D of the formwith a real number. The system is subject to the incompressibility condition, the homogeneous no-slip boundary condition, the inflow condition and the initial condition, where , is the fluid velocity vector field with components , and at the point and time . Note that , , , , , and Here is the body force per unit mass (note that we will assume that ). and are the rigid part and the inflow part of the boundary , respectively.

3. Classical Model of Navier-Stokes Equations

3.1. Model Description and Definitions

The theory of mechanics of continuous media, also known as continuum mechanics, allows the description of the constitutive equations laws that describe the deformations of fluid medium. These laws, in combination with the general conservation principles (conservation of mass and of momentum), form the system of partial differential equations, which are equal in number to the number of unknowns of the system. Namely, for 3D motion there are four dependent variables: , and , and four independent variables: , and . This leads to the following classical Navier-Stokes equations, obtained from the generalized version (5) by putting , and of which we discuss a possible treatment using the method of finite elements technique:with the incompressibility condition, boundary condition, inflow condition, and initial condition remaining the same as (6), (7), (8), and (9), respectively. In addition we assume that does not change in time.

The finite element discretization of the Navier-Stokes problem is based on the variational formulation, and the use of Sobolev spaces is needed for the mathematical treatment of the variational formulation of the model. We use subspaces of the usual Hilbert space of square-integrable functions on , where integration is in the sense of Lebesgue:and the corresponding inner products and normsNext, for any nonnegative integer , we define the Sobolev space of square-integrable functions, all of whose derivatives of order up to are also square integrable, where denotes any and all derivatives of order . comes with the norm The following definitions can now be stated: Of particular interest is the subspace of of defined by whose elements vanish on the boundary .

For functions belonging to , the seminorm defines a norm equivalent to (16). The proof of this statement is not our aim in this dissertation, but it can be found in Braess [15]. Thus for such functions, (18) may be used instead of (16).

We denote by the dual space consisting of bounded linear functionals on ; that is, implies that for all . A norm for is given by Since the velocity field is a vector valued function, we use the spacesFor , is equipped with the normAlternatively, for functions belonging to , we may use The inner product for functions belonging to is also given by Before stating the weak variational formulation for our model, let us discuss some preliminaries concerning the existence and uniqueness of a solution of the Navier-Stokes equations.

3.2. Existence and Uniqueness for a Solution of Navier-Stokes Equations

From the mathematical point of view, two questions concerning the Navier-Stokes equations are of main interest. Given a set of data which are sufficiently smooth,does a solution of (10), (6)–(9) exist?if a solution exists, is it unique?First, we have to clarify the notion of a solution of (10), (6)–(9). There exist several concepts of the notion of a solution of the above system, the most important of which are the classical solution and the weak solution.

Definition 2 (classical solution). A pair () is called a classical solution of the Navier-Stokes problem (10), (6)–(9) if() satisfies the Navier-Stokes problem (10), (6)–(9), and are infinitely many times differentiable with respect to space and time variables.

Then, according to John and Kaya [16], the existence of a classical solution of (10), (6)–(9) cannot yet be proven, but if a classical solution exists, it is unique.

To define a weak solution, we first need to transform (10) into a weak form by(i)multiplying (10) with a suitable vector valued function (test function),(ii)integrating over ,(iii)applying integration by parts (Green’s theorem).The last step is possible only if there are some restrictions on the domain. For the test function , one requires(i) for each time , where ,(ii) is infinitely differentiable with respect to time,(iii).This gives the weak formulation of the Navier-Stokes equations which has the following features:(i)there is no time derivative of ,(ii)there is no second order spatial derivative with respect to ,(iii)the pressure vanishes, since Green’s formula yields because on and .

Definition 3 (weak solution). A function is called weak solution of the Navier-Stokes equations if satisfies (24) for all test functions with the properties on given above, has the following regularity: where the subscript div means space of divergence-free functions; for instance, More generallyis the space of strongly measurable maps , such that and is a Banach space. Furthermore with It is obvious that all these spaces are needed for the weak formulation given in the next section.

The existence of a weak solution of (10), (6)–(9) was proved in 1934 by Leray [17]. The weak solution is unique if every other weak solution satisfies an additional regularity assumption, Serrin’s condition; see Serrin [7] or John and Kaya [16]. But it is not known in 3D if every weak solution possesses such additional condition.

According to the same article [16], the existence of a weak solution of the Navier-Stokes equations can be proven in arbitrary domains, but the uniqueness cannot yet be proven.

The answer to the question of uniqueness of the weak solution in 3D or existence of a classical solution in 3D is one of the major mathematical challenges of this century (John and Kaya [16]). There is a prize of one million US Dollars for people who can answer these questions.

3.3. The Stationary Case

We now consider our domain defined in the first chapter and assume that all the variables in the system (10), (6)–(9) are independent of time. We therefore obtain the following stationary Navier-Stokes problem:with the following boundary conditions:having the following weak formulation.

Find functions and such that discretized as follows.

Find a pair such that

For the sake of a stable approximation as , it is crucial that we relate the continuous and discrete spaces by the following hypotheses.

Hypothesis (approximation property of ). There exists an operator and an integer such that

Hypothesis (approximation property of ). There exists an operator such that

Hypothesis (uniform inf-sup condition ). For each , there exists a such thatwhere the constant is independent of , and ; is the space of linear operators from to ; and are the standard norms in and , respectively; is the standard seminorm in .

We assume that the bilinear forms and and the approximating subspaces and satisfy all the hypotheses , , required for suitable approximations. For any , the system (35) has a solution . It is well known (see Girault and Raviart [4] or Gunzburger [18]) that the solution is unique for “sufficiently small” data f or “sufficiently small” Reynolds number, Re. In fact if we set.the space of divergence-free functions, and then given , and if , then the problem (34) has a unique solution

From the previous assertion, we can state one of the fundamental properties (see Joseph [5]) of the solutions of the Navier-Stokes equations (32) and (33) which reads as follows: when the viscosity is large (or the Reynolds number is small), all solutions of the Navier-Stokes equations tend to a single basic flow. So what is the final destiny of all these uniquely determined solutions of the Navier-Stokes equations? For large values of the Reynolds number, the final set of flows which evolve from a given set of initial fields is generally “turbulent.”

As just proven, the Navier-Stokes equations have in general more than one solution, unless the data satisfies very stringent requirements. However, it can also be shown that in many practical examples these solutions are mostly isolated; that is, there exists a neighborhood in which each solution is unique. Furthermore, it can be established (Girault and Raviart [4]) that the solutions depend continuously on the Reynolds number, Re (which is inversely proportional to the kinematic viscosity). Thus as the Reynolds number varies along an interval, each solution of the Navier-Stokes equations describes an isolated branch. In particular, this means that the bifurcation phenomenon can be rare or occurs infrequently (Girault and Raviart [4]). This situation, frequently encountered in practice, is expressed mathematically by the notion of branches of nonsingular solutions. So the solutions of the problem (32) and (33) are “in general” nonsingular.

But when the Reynolds number Re is large (small viscosity), compared to the other parameters of the fluid, there arises a boundary layer in the neighborhood of where the viscosity predominates while it is negligible in the interior of . At the same time, the transition to turbulence occurs. Thus the solutions of the Navier-Stokes equations are seriously discontinuous at certain values of the Reynolds number, Re.

In the coming sections, we show that, under hypotheses , , , the problem (35) possesses a branch of nonsingular solutions that, as , converges to a given branch of nonsingular solutions of (34). This convergence depends on the Reynolds number Re. We also use a finite element method, the Newton method, to show that the stability of such a convergence is guaranteed for small variations of Re.

3.4. Nonsingular Solutions and Approximation

Let and be two Banach spaces and a compact interval of the real line . We are given a -mapping () and we want to solve the equation that is, we want to find pairs which are solutions of (42).

Let be a branch of solutions of (42). This means that Moreover, we suppose that these solutions are nonsingular in the sense thatAs an immediate consequence of (44), it follows from the implicit function theorem (see [4]) that is a -function from into .

Let us show that our problem for the Navier-Stokes equations (34) fits into the above framework. We first setand we introduce the intermediate space Next we define a linear operator as follows: given , we denote by the solution of the Dirichlet problem for the Stokes equations:Finally, with the data , we associate a -mapping from into , defined for a , by and we set It is clear that and (48) becomes or simply Now we may state the lemma.

Lemma 4. The pair is a solution of problem (32) and (33) if and only if , with , is a solution of (42), where the spaces and are defined by (45) and the compound mapping is defined by (49) and (50).

Proof. If () is a solution of the problem (32) and (33) then From both (47) and , applied to and , we can write or Then, (51) yields and (49) yields which means From Lemma 4, it is clear that if () is a solution of the problem (32) and (33), then , where , is a nonsingular solution of (42).

We may now state the following proposition.

Proposition 5. Assume that hypotheses , , and hold. Let be a branch of nonsingular solutions of the Navier-Stokes problem (34). Then there exists a neighborhood of the origin in and for sufficiently small a unique branch of nonsingular solutions of problem (35) such thatMoreover, we have the convergence propertyIn addition, if the mapping is continuous from into for some integer with , we have for all Our goal is not to prove this proposition, but to use the results. A similar proof of this proposition can be found in Girault and Raviart [4].

It is also possible to derive an -estimate for the velocity. But the following regularity must be satisfied: where is a closed subspace of .

Proposition 6. We retain the hypotheses of Proposition 5 and we assume that (61) holds. If the mapping is continuous from into for some integer , then we have the following -estimate for all : These propositions show that the convergence of finite element approximations is guaranteed, and the stability constants depend on .

3.5. Stability of Newton’s Method

The method discussed here is one of the finite element methods that is intended to solve the Navier-stokes equations. We saw in Lemma 4 that it suffices to investigate equations of type (42): where is defined by (49) and is a -mapping () defined on (with a compact interval of the real line and . For simplicity, let us set and fix and assume that is a nonsingular solution of (63). Then where is a Banach space.

Since is an isolated solution of (63) and since is at least differentiable, an efficient way to approximate is by the Newton method. The Newton algorithm reads as follows.

Starting from an initial guess , construct the sequence in by or equivalently As is a linear operator, each step of Newton’s method requires the solution of a different problem relative to . If this is too costly, the simplest alternative is to replace (66) by or equivalently The drawback of Newton’s method is that the stability of its (quadratic) convergence (see [4]) can only be ensured when the first guess is sufficiently near the solution . If this solution is part of a branch of nonsingular solutions, and if we know the solution at a neighboring point, say for an adequate increment , then we can derive from this value the first guess to start Newton’s algorithm.

Since is a branch of nonsingular solution of (63), then is a -mapping (), so is the mapping and we can differentiate both sides of (63):that is, we find a first order differential equation of the form where The simplest way to solve (71) is to use the one-step, explicit Euler’s method; this brings us to choose In other words is defined by Let us estimate the difference . From (71), we infer thatSubtracting (73) from this equality yields Hence Thus is and if is small enough, the solutions and stay close, characterizing the stability, and defined by (73) is an adequate starting value for Newton’s algorithm.

4. Existence Results: Discussion and Concluding Remark

We have analyzed existence of nonsingular solution to NSEs of type (5), with , and used Newton’s method to study its stability. It happens that the stability of its (quadratic) convergence, when exploiting Newton’s method, can only be ensured when the first guess is sufficiently near the solution Now what happens for the system (5) with a general Does a nonsingular solution exist? If yes, in which spaces? We know [19] that there is global existence for the fractional NSEs (5), with only and the system is well-posed in some supercritical Besov spaces as well as in the largest critical spaces It has also been shown [20] that any classical solution to NSEs (5) with is always global in time. Thus, it is clear that investigations on the system of NSEs of type (5) are far from being complete and the work developed here will contribute to future possible investigation in order to study, for a general parameter , the existence and well-posedness of the fractional Navier-Stokes equations in some Banach spaces.

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.