In this article, a well known model is to be investigated that represents the flow of an incompressible fluid in both bounded and unbounded domains of . This model is commonly called the Navier-Stokes equation following the French engineer Navier who was the first to propose this model. This model was investigated later by Poisson and de Saint Venant. However, Stokes was the one who justified the model based on the principles of continuum mechanics. By advent of 1930, the interest in this model increased significantly and outstanding results were obtained by Leray, Hopf, Ladyzhenskaya and Finn.
This equation describes the flow of what is so called the Newtonian fluid. These are the fluids that exhibit shearing stress due to the presence of frictional forces. Frictional forces within fluids are consequences of its viscosity. Also, the gradient of the velocity represents a measure for the relative motion of the fluid’s particles. Moreover, deformation of fluids is commonly associated with internal and external body forces; the internal force is what we refer to as the pressure of the fluid. The derivation of the Navier-Stokes equation is a natural application of Newton’s second law of motion, the balance of momentum and the mass conservation, which eventually leads to the definition of the Cauchy stress tensor. In Newtonian fluids, this stress tensor is a function in the pressure, the viscosity and the gradient of the velocity. For a convenient physical background about the basics of continuum mechanics and how we derive the Navier-Stokes equation we propose   . Also a very interesting work from a physical point of view can be found in   . In particular, the work of Kambe in  was the source of inspiration for the ideas in this article.
This model poses a serious challenge when it comes to proving the existence and the smoothness of its solution. This problem was perfectly addressed by Ladyzhenskaya in two dimensional spaces among many other issues in higher dimensional spaces  . However, a decisive answer in the three-dimensional space or higher remains unavailable. It is almost impossible to enlist all the results obtained for this equation. Therefore, we suggest for the interested reader to review the monographs    and the references within for much more details.
Recently, the interest in this equation is not fading at all. There are persistent efforts to clarify the properties of the solution, especially its smoothness. Among many respectful results, we mention the outstanding analysis by Tao in  , the work of Constantin in    . A very interesting result for partial regularity of suitable weak solutions was obtained by Caffarelli in  .
In this article the idea is simple. A class of possible solutions is proposed and then it is proved that it indeed represents the unique classical solution of the problem. Most of the results are obtained by considering standard theories of partial differential equations. Some of the results in the monograph  are also used repeatedly. In the next section a statement of the problem is introduced along with some definitions, notations and employed functional spaces. After- wards, the proofs of the main results are established.
2. Statement of the Problem
The spatial domain is which is either a bounded region in or the whole of and this point shall be specified explicitly. For the sake of conciseness, the notation is used to denote . Clearly, such notation should not be taken to imply a moving boundary. The main model equation is in the form
where the last equation in the above model is what many authors commonly refer to as the incompressibility condition or the solenoidal condition. The first term in the first equation is the acceleration of the fluid’s flow in time, the second is the convective term that represents the acceleration of the flow in space, the third represents the diffusion scaled by the kinematic viscosity constant , the fourth is the pressure, and the last one represents the total of the external body forces. The initial profile is denoted by and the boundary datum is denoted by . The solution is the vector field representing the velocity of the flow in each direction, and its rotation is the vorticity. Note that in by compatibility.
The well known Lebesgue spaces will be used repeatedly to represent the functions with bounded mean of order q. The Sobolev space is used to represent functions with bounded derivatives such that for a vector field one has for every and . This motivates the usage of the space , which is a well known space of functions in the theory of incompressible fluids as a representative for divergence free (solenoidal) bounded vector fields such that .
By laws of classical mechanics, the energy generated by a moving object is proportional to the square of its velocity. Hence, the energy generated by the flow is defined as follows
Recall that the above integral represents the norm of in the Lebesgue space . The smoothness of is such that
The smoothness of the boundary datum is such that
Finally, the forcing term is smooth in space and time such that
Note that the intersections in the above conditions are not really required in the case of bounded domain since boundedness of the domain and continuity of the functions are enough to imply boundedness in the sense of the mean. However, these requirements are of significant importance in the case of unbounded domain as will be shown later.
The target is to define a class of possible solutions to Model Equation (1) from which , and can be concluded. Once is obtained, then can easily be recovered from the main model. The validity of this solution as a meaningful physical solution will be investigated when inserted in the main model. A meaningful solution is a unique and smooth solution that vanishes as , and in the case of unbounded domain it vanishes as as well.
Remark. The curl operator or the rotation of a vector field has a physical meaning only in three dimensional space. However, it will be used in for the sake of generality. Most of the results depend on the curl operator in the sense of a differential operator without direct exposure to its definition. Note that the main interest is to find the velocity, which means that any use of the rotation, in spite of its significance in this work, is nothing more than a transient step. It can always be assumed that the space is three dimensional when needed, and a generalization becomes possible by reverting back to the results of the velocity. In particular, some of the vorticity ideas introduced in  are adopted in this study.
3. Main Results
In this section a class of possible solutions is proposed and the insertion of these solutions in the main model is investigated to check where they will lead to. This is initiated by the statement of the following claim.
Claim 1. The unique solution of Model Equation (1) is in the form where is a scalar field and is a vector field such that, at least, and .
An important question in the theory of Navier-Stokes equation is the ability to verify the compatibility condition on the boundary with minimum restrictions on the flux passing through the boundary especially if is divided into several parts (  , p. 4-8). This condition is a natural consequence of the incompressibility of the flow. Hence, it takes the form
where is the outward unit vector normal to . This motivates the introduction of the following lemma.
Lemma 1 (Tangential flow). Let be an arbitrary domain, be any scalar field such that and let be any vector field independent of . The Compatibility Condition (6) is satisfied for every divergence free vector field in the form . In particular, on every part of , and its rotation are tangents to such that .
Proof. Given that is divergence free such that leads to
Identity (7) can be used to deduce that . Upon dot product by one obtains
and when integrated over for any provides
which is a true identity for every arbitrary , and and for every . On the other hand, given that , one can use the identity
where , which upon dot product by provides
where we used also that . Integrate over for any to get
which is true for arbitrary , and . Now, since one also has
Identity (10) implies that is orthogonal to the space spanned by and . Combine Identities (9), (11) and (12), and exclude the cases and by the arbitrariness of the choice to deduce that it is necessary that ; that is to say is tangential to every part of . It can also be deduced that either on every part of (this actually means on ), or on every part of . Both cases imply that on every part of . Hence, both and are tangential to the boundary.
Remark. As pointed out, the question of verifying the Compatibility Condition (6) on the boundary of is an important open question in the mathematical theory of Navier-Stokes equation. Some results were obtained to justify the validity of such compatibility under certain restrictions on the flux of the flow in terms of the viscosity of the fluid, for details on this issue refer to  . In the present case, the compatibility is naturally achieved given, of course, that the solution of Model Equation (1) is indeed in the form proposed in Claim 1. This shall be verified by the statements of the following theorems.
Theorem 2 (Bounded Domain). Let be a bounded domain with sufficiently smooth boundaries and let . Suppose , and satisfy Conditions (3), (4) and (5) respectively. If is in the form proposed in Claim 1, then Model Equation (1) has a classical solution with bounded energy such that , and . In particular, the exact solution is given by solving the following system
where can be defined uniquely in terms of the values of and on the boundary. Moreover, if then .
Proof. The proof is quite simple and it depends mostly on classical results and the standard theory of linear parabolic and elliptic equations of second order. If , then by virtue of Identity (7) one has . Hence, the main equation takes the form
Apply the divergence operator to get
where the incompressibility condition is used. By the standard theory of elliptic equations, if satisfies Condition (5), then . However, if , then is actually , for details on such equation see (  , pp. 326-343). This concludes one part of System (13), however, a further discussion on a unique definition of will be introduced at the end of this proof.
Now, revert to Equation (14) and apply the curl operator to get
which is the first equation in System (13). Finally, given the incompressibility of and applying a simple vector identity lead to the third equation in System (13) that is
The fundamental solution to Equation (16) with initial profile and force is given formally by
As explained in (  , Chapter 7), because of the smoothing property of the Gaussian kernel, it is enough to have a contentious data under the integral sign to guarantee that , which is what has been already assumed. Given the assumptions on the growth of in time and the form of Formula (18), the solution as . Moreover, the continuity of the data, the clear decay in time and the assumed boundedness of imply that in , which implies actually that for every . Uniqueness of as per Expression (18) is not clear unless we insert the boundary datum explicitly in Expression (18). This can be done by introducing the auxiliary variable for which one obtains a homogeneous heat equation. Anyhow, the presence of any form of boundary conditions guarantees the uniqueness of Moreover, since is actually a derivative of then it suffices to show that is unique, which is our main concern, to conclude the uniqueness of
Now, go back to Equation (17). By virtue of the results obtained above for and the standard theory of elliptic equations one directly concludes that , for more details see (  , pp. 326-343). Such regularity, the boundedness of and the global decay of in time imply that is bounded in for every , which in turn implies the boundedness of the energy of the flow as defined by Expression (2). Since , then (  , Proposition 2.16]) can be used to conclude that a formal solution for Equation (17) takes the form
where is the area of the unit sphere in . The solution can be enforced to take the values on the boundary in a standard manner by introducing Dirichlet Green’s function. We refrain from discussing these details being highly dependent on the choice of the domain. The uniqueness of follows by the presence of the boundary condition , for details see   .
Finally, go back to Equation (15) to solve for . In this case one only needs to calculate from the main model by knowing the values of and on the boundary. This provides a form of boundary conditions for which consequently guarantees its uniqueness up to a constant, for details see  . The infinite differentiability of also follows from the main model and the fact that provided, of course, that as well. If , then is certainly . This completes the proof.
Remark. As explained in the proof of Theorem 2, the uniqueness of follows by the ability to define uniquely. The order of solving the equations in System (13) is not really important since none of the quantities induces the other; they act simultaneously. Another way of solvability can be introduced by which one can obtain the same results. This can be a topic for a future study.
The problem of proving the existence of regular and smooth enough solutions for the Navier-Stokes equation in bounded domain was exhaustively investigated as pointed out in the introduction. The real problem was to prove the boundedness of the solution in an unbounded domain, clearly because of the unboundedness of the domain itself. This fact manifests the need to show that the solution’s support is bounded in , or equivalently to show that the solution decays rapidly as .
The solution obtained in Theorem 2 represents a perfect candidate as a solution for Model Equation (1) in unbounded domains also except for one issue. One needs to prove the boundedness of in for every so that the boundedness of the energy can be claimed, and also to ensure that the solution does indeed vanish as . For to be bounded in , it should attain a rate of decay, at least, for any . One can argue that some of the results in the literature require a rate of decay higher than that for the surface integrals to vanish; these restrictions can be dropped because these integrals already vanish by virtue of Lemma 1, (see  , Lemma 1.5]). However, this does not mean that rapid rates of decay are not achievable, they are achievable as demonstrated next.
In order to derive such an estimate one goes back to Formula (19) that represents the fundamental solution for . If then it is expected that the outcome of this integral will provide nothing less than a linear rate of growth for , which is a bad answer to the problem in hand. Therefore, Formula (18) shall be used to help us estimate some rates of decay for and consequently for so that boundedness in can be proved for every
Theorem 3 (The Domain RN). Suppose all the conditions of Theorem 2 are satisfied for . Then there exists a classical solution for Model Equation (1) represented by the System (13) and defined as
where , and is the surface area of a unit
sphere. The pressure can be defined from Model Equation (1) up to a constant where can be specified uniquely in terms of the values of and on the boundary. Moreover, the energy of the flow is bounded for every where grows at most as .
Proof. The proof of smoothness and uniqueness is identical to the one introduced in the proof of Theorem 2 and it follows by the standard theory of linear second order elliptic and parabolic equations. The focus here will be on proving the boundedness in for both and , which necessarily entails an estimation of appropriate decay rates as pointed out in the preceding discussion. Since , then . By assumption, which implies that and since for meaningful physical interpretation, then at least . Hence, as per the standard theory of linear second order parabolic equations, for details refer to (  , Theorem 6, p. 386). Now, it is needed to show that is bounded in for every . There are two ways to show this; in one of them a rough estimate will be provided for the minimum rate of decay of given the assumptions on the data.
The first direction depends on the results in (  , Theorems 3.4 and 3.6). In these theorems a regularization technique by mollifiers along with energy estimates were used to prove global in time existence. In particular, (  , Theorem 3.4) states that if and , then there exists a unique continuous solution locally in time such that this is true up to , which coincides with the assumptions on . The local existence of a unique continuous solution was extended to a global in time existence in (  , Theorem 3.6) given that
and such that . Since and are bounded in for every , then it follows that there exists a ball such that the supports of and are entirely inside where . Hence, the integrals in Formula (20) can be restricted to the ball , apply Hölder’s inequality, maximize the exponential term which is bounded for every and for every , and with some estimation procedures it becomes easy to find an estimate of the form in the whole of and for every . Since is a decreasing function in time then by assuming a simple scale of time one obtains
It follows that all the conditions of (  , Theorem 3.4, Theorem 3.6) are satisfied such that exists globally in time and such that for every , which implies the boundedness of the energy for every as well.
The second way is trying to get an estimate for in terms of to confirm the boundedness in . Consider the following argument: fix in Expression (20), calculate , which is in for every because . Take the absolute value of both sides, perform some manipulation to the integrands, use and maximize the time integral (the integrand is a decreasing function in time) so that one finally gets the term with the highest power for as follows
where is the collection of every possible appearance of any power of after the integration. Since is in as pointed out, then to obtain a decreasing integrand when calculating the -norm of this derivative it is necessary that is decreasing in such that it is at least
for some . But is in as well, and the integrands of the above estimate are the same except for the terms with positive powers of . That is to say that . Incorporate this estimate in Expression (21) of the solution and one can readily see that as desired, which in turn confirms the boundedness of the energy for every and implies the decay of as .
However, since , then we actually have for every . This means that we can differentiate Formula (20) four times and repeat the same argument as above to conclude that we actually have
which in turn leads to and this
completes the proof.
Remark. Better estimates for the decay of can be obtained by repeating the same procedure described above for higher order derivatives of which, of course, entails higher assumptions on the data so that we can claim the boundedness of the considered derivative apriori. On the other hand, the assumptions on the smoothness of the forcing term can be relaxed. It is enough to assume that so that one obtains , which is sufficient to conclude that . The boundedness of in is enough to deduce a sufficient decay estimate for as shown above. The introduced assumptions were chosen for consistency with the standard theory and the Embedding Theorem, and in the same time to illustrate the estimation procedure. For more information review (  , pp. 382-386).
We managed to prove that our solution is indeed a classical solution of Model Equation (1). Here come some important questions, what if there exists another solution in a more general form? Moreover, does the choice of the domain or the choice of the boundary data play a role in the uniqueness of the solution? The answer to these questions is addressed by the statement of the next theorem.
Theorem 4 (Uniqueness). Let be an arbitrary domain. Suppose , and are satisfying Conditions (3), (4) and (5) respectively. Then Claim 1 is true and the unique classical solution of Model Equation (1) is in the form . This solution is defined as per Theorems 2 and 3.
Proof. The proof here depends on our results in Theorems 2 and 3, and also on (  , Theorems 3.4 and 3.6). By virtue of the assumptions on the data and our results in Theorems 2 and 3, it is clear that our solution satisfies all the conditions in  for every .
Assume that there exists a more general solution than the proposed one and denote it by . Such solution should definitely inherit the smoothness proved in  as well. That is to say that and such that is the triad solution of Model Equation (1) with boundary datum and initial profile . Let and let . Hence, has zero boundary and initial data and it obeys the equation
Dot product the above equation by and integrate by parts over any arbitrary domain to get
where we used the Divergence theorem in the right hand side, the facts that and that on the boundary, and the results of Lemma 1. Now, recalling that and using the vector identity
one can dot product Equation (22) by and integrate as above to get
where Lemma 1 is used again. Identities (23) and (25) imply one of three possibilities. Either which is excluded for being trivial, or almost everywhere. The third possibility is being orthogonal to the space spanned by and , which by Equation (22) implies that is orthogonal to and .
Let us start with almost everywhere. Multiply this equation by , integrate by parts over , employ the Divergence theorem and recall that on to get
This readily implies that almost everywhere. But by the results obtained for and , one concludes that at least which implies that is identically zero. Hence, the solution is the unique solution for Model Equation (1) in this case.
Now, if and are orthogonal to , then by Identity (24) we deduce that is orthogonal to the space spanned by and everywhere, which is equivalent to the nature of the solution . This means that and are parallel to each other on the boundary. Since both solutions coincide at the boundary and both are extended continuously to the interior of the domain, then it is not hard to conclude that everywhere, which is the aim of this proof. Assume not. Let so that the difference equation for takes the form
where . Use the vector identity
and given the incompressibility of one also has
Use these two identities, bearing in mind the incompressibility, to find that
Now recall that we are discussing the possibility in which is orthogonal to and such that . Following the same steps as above one can write
Now, dot product Equation (26) by , use Identity (27), integrate over and use the Divergence theorem to reach
Also dot product Equation (26) by , use Identity (28), integrate over and use the Divergence theorem to get
Now, reverse the Divergence theorem in the surface integrals that include values of and and reverse all the steps made to conclude and to obtain identities in the form
Now, there are three other possibilities. The trivial solution; that is and it is excluded. Another one is and this one is equivalent to because the uniqueness of implies the uniqueness of and vice versa. It remains that is orthogonal to the space spanned by and . But we already have orthogonal to and , which implies that and its rotation (curl) are parallel to each other. By definition, the curl operator is the unique vector field for which for every vector field (  , p. 32). If , then let to get which in turn implies that then as well. Hence, it is necessary that which implies that . This completes the proof.
Remark. In the proof of uniqueness one can argue that the statement of the proof was given in the sense of classical solutions, and such that there still exists another solution in a weaker form. While this may sound true, but in fact it is not. The basic idea of the proof is based on using the coincidence on the boundary and then moving back to the interior. Assuming the existence of a weaker solution does not change the fact that it is going to coincide with the proposed one on the boundary. However, further investigation on this specific point will be introduced in a future study.
4. Conclusions and Suggestions
A class of possible solutions for the incompressible Navier-Stokes equation in bounded and unbounded regions of is proposed. It was demonstrated that for such class of vector fields, the flux of the energy of the flow is orthogonal to the space spanned by and its associated rotation . It was also proved that and are tangential to the boundary such that .
An investigation of the validity of this class of vector fields as a candidate for a solution to the incompressible Navier-Stokes equation was carried out in both bounded and unbounded domains. Given plausible assumptions on the data and a forcing term in the case of unbounded domains, it turned out that this class of solutions represents perfectly a classical solution of the problem. Verification was established for the infinite differentiability, the uniqueness and the boundedness of the energy in appropriate spaces in light of well known and standard theories. An appropriate estimate was also given for the minimum rate of decay of the solution as . Moreover, global existence in time and the corresponding rate of decay were quite obvious in the deduced formulas.
Finally, it was proved that this class of solutions represents actually the unique classical solution of the incompressible Navier-Stokes equation. In light of the non-convective nature of the proposed solution and the uniqueness argument under arbitrary settings, the incompressible Navier-Stokes equation can safely be reduced to a linear equation. This point is quite interesting and it motivates further investigations on a possible relation between incompressibility and convection in fluid mechanics.
Sincere thanks to the members of JAMP for their professional performance, and special thanks to managing editor Hellen Xu for a rare attitude of high quality.
 Ladyzhenskaya, O. (1969) The Mathematical Theory of Viscous Incompressible Flow, Translated from the Russian by Richard A. Silverman, Mathematics and its Applications. 2nd Edition, Gordon and Breach.
 Constantin, P. and Fefferman, C. (1993) Direction of Vorticity and the Problem of Global Regularity for the Na-vier-Stokes Equations. Indiana University Mathematics Journal, 42, 775-789.
 Caffarelli, L., Kohn, R. and Nirenberg, L. (1982) Partial Regularity of Suitable Weak Solutions of the Navier-Stokes Equations. Communication on Pure & Applied Mathematics, 35, 771-831.