View source for Symbol of a differential operator

The intent of this article is to highlight the important points of the '''derivation of the [[Navier–Stokes equations]]''' as well as the application and formulation for different families of [[fluid]]s.

==Basic assumptions==

The Navier–Stokes equations are based on the assumption that the fluid, at the scale of interest, is a continuum, in other words is not made up of discrete particles but rather a continuous substance. Another necessary assumption is that all the fields of interest like [[pressure]], [[velocity]], [[density]], [[temperature]] and so on are [[differentiable]], [[Weak derivative|weakly]] at least.

The equations are derived from the basic principles of [[conservation of mass]], [[Conservation of momentum|momentum]], and [[Conservation of energy|energy]]. For that matter, sometimes it is necessary to consider a finite arbitrary volume, called a [[control volume]], over which these principles can be applied. This finite volume is denoted by <math>\Omega</math> and its bounding surface <math>\partial \Omega</math>. The control volume can remain fixed in space or can move with the fluid.

==The material derivative==
{{main|material derivative}}

Changes in properties of a moving fluid can be measured in two different ways. One can measure a given property by either carrying out the measurement on a fixed point in space as particles of the fluid pass by, or by following a parcel of fluid along its [[Streamlines, streaklines and pathlines|streamline]]. The derivative of a field with respect to a fixed position in space is called the ''Eulerian'' derivative while the derivative following a moving parcel is called the ''advective'' or ''material'' ("Lagrangian" <ref>Munson, Bruce R. Fundamentals of Fluid Mechanics. 7th. Jefferson City: John Wiley and Sons, Inc., 2013. Print.</ref>) derivative.

The material derivative is defined as the operator:

:<math>\frac{D}{Dt} \ \stackrel{\mathrm{def}}{=}\  \frac{\partial}{\partial t} + \mathbf{v}\cdot\nabla </math>

where <math>\mathbf{v}</math> is the velocity of the fluid. The first term on the right-hand side of the equation is the ordinary Eulerian derivative (i.e. the derivative on a fixed reference frame, representing changes at a point with respect to time) whereas the second term represents changes of a quantity with respect to position (see [[advection]]). This "special" derivative is in fact the ordinary derivative of a function of many variables along a path following the fluid motion; it may be derived through application of the [[chain rule]] in which all independent variables are checked for change along the path (i.e. the [[total derivative]]).

For example, the measurement of changes in wind velocity in the [[Earth's atmosphere|atmosphere]] can be obtained with the help of an [[anemometer]] in a weather station or by observing the movement of a weather balloon. The anemometer in the first case is measuring the velocity of all the moving particles passing through a fixed point in space, whereas in the second case the instrument is measuring changes in velocity as it moves with the fluid.

==Conservation laws==
The Navier–Stokes equation is a special case of the (general) [[continuity equation]]. It, and associated equations such as mass continuity, may be derived from [[conservation laws|conservation principles]] of:

* [[Mass]]
* [[Momentum]]
* [[Energy]].

This is done via the [[Reynolds transport theorem]], an integral solution relation stating that the sum of the changes of some [[intensive property]] (call it <math>L</math>) defined over a control volume <math>\Omega</math> must be equal to what is lost (or gained) through the boundaries of the volume plus what is created/consumed by sources and sinks inside the control volume. This is expressed by the following integral equation:

:<math>\frac{d}{dt}\int_{\Omega} L \ dV = -\int_{\partial\Omega} L\mathbf{v\cdot n} \ dA - \int_{\Omega} Q \ dV</math>

where '''v''' is the velocity of the fluid and <math>Q</math> represents the sources and sinks in the fluid, taking the sinks as positive. Recall that <math>\Omega</math> represents the control volume and <math>\partial \Omega</math> its bounding surface.

The [[divergence theorem]] may be applied to the [[surface integral]], changing it into a [[volume integral]]:

:<math>\frac{d}{d t} \int_{\Omega}  L \ dV = -\int_{\Omega} \nabla \cdot ( L\mathbf{v}) \ dV - \int_{\Omega} Q \ dV.</math>

Applying [[Leibniz integral rule|Leibniz's rule]] to the integral on the left and then combining all of the integrals:

:<math>\int_{\Omega} \frac{\partial L}{\partial t} \ dV = - \int_{\Omega}\nabla \cdot (L\mathbf{v}) \ dV - \int_{\Omega} Q \ dV
\qquad \Rightarrow \qquad
\int_{\Omega} \left( \frac{\partial L}{\partial t} + \nabla \cdot (L\mathbf{v}) + Q\ \right) dV = 0.</math>

The integral must be zero for '''any''' control volume; this can only be true if the integrand itself is zero, so that:

:<math>\frac{\partial L}{\partial t} + \nabla \cdot (L\mathbf{v}) + Q = 0.</math>

From this valuable relation (a very generic [[continuity equation]]), three important concepts may be concisely written: conservation of mass, conservation of momentum, and conservation of energy. Validity is retained if <math>L</math> is a vector, in which case the vector-vector product in the second term will be a [[Dyadic product|dyad]].

===Conservation of momentum===
The most elemental form of the Navier–Stokes equations is obtained when the conservation relation is applied to momentum. Writing momentum as <math>\rho \mathbf{v}</math> gives:

:<math>\frac{\partial}{\partial t}(\rho \mathbf{v}) + \nabla \cdot (\rho \mathbf{v} \mathbf{v}) + \mathbf{Q} = 0</math>

where <math>\mathbf{v} \mathbf{v}</math> is a [[Dyadic product|dyad]], a special case of [[tensor product]], which results in a second rank tensor; the [[divergence]] of a second rank tensor is again a vector (a first rank tensor).<ref name="TA">
{{cite book
  | last = Lebedev
  | first = Leonid P.
  | authorlink =
  | coauthors =
  | title = Tensor Analysis
  | publisher = World Scientific
  | year = 2003
  | location = 
  | pages = 
  | url =
  | doi =
  | id =  
  | isbn = 981-238-360-3}}
</ref> Noting that a [[body force]] (notated <math>\mathbf{b}</math>) is a source or sink of [[momentum]] (per volume) and expanding the derivatives completely:

:<math>\mathbf{v} \frac{\partial \rho}{\partial t} + \rho \frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \mathbf{v} \cdot \nabla \rho + 
\rho \mathbf{v} \cdot \nabla \mathbf{v} + \rho \mathbf{v} \nabla \cdot \mathbf{v} = \mathbf{b}</math>

Note that the [[gradient]] of a vector is a special case of the [[covariant derivative]], the operation results in second rank tensors;<ref name="TA"/> except in Cartesian coordinates, it's important to understand that this isn't simply an element by element gradient. Rearranging and recognizing that <math>\mathbf{v} \cdot \nabla \rho + \rho \nabla \cdot \mathbf{v} = \nabla \cdot (\rho \mathbf{v})</math>:

:<math>\mathbf{v} \left(\frac{\partial \rho}{\partial t} + \mathbf{v} \cdot \nabla \rho + \rho \nabla \cdot \mathbf{v}\right) + \rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = \mathbf{b}</math>

:<math>\mathbf{v} \left(\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v})\right) + \rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = \mathbf{b}</math>

The leftmost expression enclosed in parentheses is, by mass continuity (shown in a moment), equal to zero. Noting that what remains on the right side of the equation is the [[convective derivative]]:

:<math>\rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = \mathbf{b}
\qquad \Rightarrow \qquad
\rho\frac{D \mathbf{v}}{D t} = \mathbf{b}</math>

This appears to simply be an expression of [[Newton's second law]] ('''F''' = m'''a''') in terms of [[body force]]s instead of point forces. Each term in any case of the Navier–Stokes equations is a body force. A shorter though less rigorous way to arrive at this result would be the application of the [[chain rule]] to acceleration:

:<math>\rho \frac{d}{d t}(\mathbf{v}(x, y, z, t)) = \mathbf{b}

\qquad \Rightarrow \qquad

\rho \left(
\frac{\partial \mathbf{v}}{\partial t} + 
\frac{\partial \mathbf{v}}{\partial x}\frac{d x}{d t} + 
\frac{\partial \mathbf{v}}{\partial y}\frac{d y}{d t} + 
\frac{\partial \mathbf{v}}{\partial z}\frac{d z}{d t} 
\right) = \mathbf{b} \qquad \Rightarrow</math>

:<math>\rho \left(
\frac{\partial \mathbf{v}}{\partial t} + 
u \frac{\partial \mathbf{v}}{\partial x} + 
v \frac{\partial \mathbf{v}}{\partial y} + 
w \frac{\partial \mathbf{v}}{\partial z} 
\right) = \mathbf{b}

\qquad \Rightarrow \qquad

\rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = \mathbf{b}</math>

where <math>\mathbf{v} = (u, v, w)</math>. The reason why this is "less rigorous" is that we haven't shown that picking <math>\mathbf{v} = \left(\frac{d x}{d t}, \frac{d y}{d t}, \frac{d z}{d t}\right)</math> is correct; however it does make sense since with that choice of path the derivative is "following" a fluid "particle", and in order for [[Newton's second law]] to work, forces must be summed following a particle. For this reason the [[convective derivative]] is also known as the particle derivative.

===Conservation of mass===

Mass may be considered also. Taking <math>Q = 0</math> (no sources or sinks of mass) and putting in density:

:<math>\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0</math>

where <math>\rho</math> is the mass density (mass per unit volume), and <math>\mathbf{v}</math> is the velocity of the fluid. This equation is called the '''mass continuity equation''', or simply "the" continuity equation. This equation generally accompanies the Navier–Stokes equation.

In the case of an [[incompressible fluid]], <math>\rho</math> is a constant and the equation reduces to:

:<math>\nabla\cdot\mathbf{v} = 0</math>

which is in fact a statement of the conservation of volume.

==General form of the equations of motion==

The generic body force <math>\mathbf{b}</math> seen previously is made specific first by breaking it up into two new terms, one to describe forces resulting from stresses and one for "other" forces such as gravity. By examining the forces acting on a small cube in a fluid, it may be shown that

:<math>\rho\frac{D\mathbf{v}}{D t} = \nabla \cdot \boldsymbol{\sigma} + \mathbf{f}</math>

where <math>\boldsymbol{\sigma}</math> is the [[Cauchy stress tensor]], and <math>\mathbf{f}</math> accounts for other body forces present. This equation is called the [[Cauchy momentum equation]] and describes the non-relativistic momentum conservation of ''any'' continuum that conserves mass. <math>\boldsymbol{\sigma}</math> is a rank two symmetric tensor given by its covariant components:

:<math>\sigma_{ij} = \begin{pmatrix}
\sigma_{xx} &  \tau_{xy} & \tau_{xz} \\
\tau_{yx} &  \sigma_{yy} & \tau_{yz} \\
\tau_{zx} &  \tau_{zy} & \sigma_{zz}
\end{pmatrix}</math>

where the <math>\sigma</math> are [[normal stress]]es and <math>\tau</math> [[shear stress]]es. This tensor is split up into two terms:

:<math>\sigma_{ij} = \begin{pmatrix}
\sigma_{xx} &  \tau_{xy} & \tau_{xz} \\
\tau_{yx} &  \sigma_{yy} & \tau_{yz} \\
\tau_{zx} &  \tau_{zy} & \sigma_{zz}
\end{pmatrix}
=
-\begin{pmatrix}
\pi &0&0\\
0&\pi &0\\
0&0&\pi 
\end{pmatrix}
+ 
\begin{pmatrix}
\sigma_{xx}+\pi &  \tau_{xy} & \tau_{xz} \\
\tau_{yx} & \sigma_{yy}+\pi  & \tau_{yz} \\
\tau_{zx} &  \tau_{zy} & \sigma_{zz}+\pi 
\end{pmatrix}
= -\pi I + \mathbb{T}
</math>

where <math>I</math> is the 3 x 3 [[identity matrix]] and <math>\mathbb{T}</math> is the [[deviatoric stress tensor]]. Note that the mechanical [[pressure]] ''π'' is equal to minus the mean normal stress:{{sfn| Batchelor|2000|pp=141}}

:<math>\pi = -\frac{1}{3} \left( \sigma_{xx} + \sigma_{yy} + \sigma_{zz} \right).</math>

The motivation for doing this is that pressure is typically a variable of interest, and also this simplifies application to specific fluid families later on since the rightmost tensor <math>\mathbb{T}</math> in the equation above must be zero for a fluid at rest. Note that <math>\mathbb{T}</math> is [[traceless]]. The Navier–Stokes equation may now be written in the most general form:

:<math>\rho\frac{D\mathbf{v}}{D t} = -\nabla \pi + \nabla \cdot\mathbb{T} + \mathbf{f}</math>

This equation is still incomplete. For completion, one must make hypotheses on the forms of <math>\mathbb{T}</math> and <math>\pi</math>, that is, one needs a constitutive law for the stress tensor which can be obtained for specific fluid families and on the pressure; additionally, if the flow is assumed compressible an equation of state will be required, which will likely further require a conservation of energy formulation.

==Application to different fluids==

The general form of the equations of motion is not "ready for use", the stress tensor is still unknown so that more information is needed; this information is normally some knowledge of the viscous behavior of the fluid. For different types of fluid flow this results in specific forms of the Navier–Stokes equations.

===Newtonian fluid===
{{main|Newtonian fluid}}

====Compressible Newtonian fluid====
The formulation for Newtonian fluids stems from an observation made by [[Isaac Newton|Newton]] that, for most fluids,

:<math>\tau \propto \frac{\partial u}{\partial y}</math>

In order to apply this to the Navier–Stokes equations, three assumptions were made by Stokes:

:* The stress tensor is a linear function of the strain rates.
:* The fluid is isotropic.
:* For a fluid at rest, <math>\nabla \cdot \mathbb{T}</math> must be zero (so that [[hydrostatic pressure]] results).

Applying these assumptions will lead to:

:<math>\mathbb{T}_{ij} = \mu\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}+\delta_{ij}\frac{\lambda}{\mu}\frac{\partial u_k}{\partial x_k}\right)</math>

That is, the deviatoric of the deformation rate tensor is identified to the deviatoric of the stress tensor, up to a factor ''μ''.

<math>\delta_{ij}</math> is the [[Kronecker delta]]. ''μ'' and ''λ'' are proportionality constants associated with the assumption that stress depends on strain linearly; ''μ'' is called the first coefficient of [[viscosity]] (usually just called "viscosity") and ''λ'' is the second coefficient of viscosity (related to [[bulk viscosity]]). The value of ''λ'', which produces a viscous effect associated with volume change, is very difficult to determine, not even its sign is known with absolute certainty. Even in compressible flows, the term involving ''λ'' is often negligible; however it can occasionally be important even in nearly incompressible flows and is a matter of controversy. When taken nonzero, the most common approximation is ''λ''&nbsp;≈&nbsp;-&nbsp;⅔&nbsp;''μ''.{{sfn |Batchelor|2000|pp=144}}

A straightforward substitution of <math>\mathbb{T}_{ij}</math> into the momentum conservation equation will yield the '''Navier–Stokes equations for a compressible Newtonian fluid''':

:<math> \rho \left(\frac{\partial u}{\partial t} + u \frac{\partial u}{\partial x} + v \frac{\partial u}{\partial y} + w \frac{\partial u}{\partial z}\right) =  -\frac{\partial p}{\partial x} + 
\frac{\partial}{\partial x}\left(2 \mu \frac{\partial u}{\partial x} - \frac{2\mu}{3} \nabla \cdot \mathbf{v}\right) + 
\frac{\partial}{\partial y}\left(\mu\left(\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}\right)\right) + 
\frac{\partial}{\partial z}\left(\mu\left(\frac{\partial u}{\partial z} + \frac{\partial w}{\partial x}\right)\right) + 
\rho g_x</math>

:<math> \rho \left(\frac{\partial v}{\partial t} + u \frac{\partial v}{\partial x} + v \frac{\partial v}{\partial y}+ w \frac{\partial v}{\partial z}\right) = -\frac{\partial p}{\partial y} + 
\frac{\partial}{\partial x}\left(\mu\left(\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y}\right)\right) + 
\frac{\partial}{\partial y}\left(2 \mu \frac{\partial v}{\partial y} - \frac{2\mu}{3} \nabla \cdot \mathbf{v}\right) + 
\frac{\partial}{\partial z}\left(\mu\left(\frac{\partial v}{\partial z} + \frac{\partial w}{\partial y}\right)\right) + 
\rho g_y</math>

:<math> \rho \left(\frac{\partial w}{\partial t} + u \frac{\partial w}{\partial x} + v \frac{\partial w}{\partial y}+ w \frac{\partial w}{\partial z}\right) = -\frac{\partial p}{\partial z} + 
\frac{\partial}{\partial x}\left(\mu\left(\frac{\partial w}{\partial x} + \frac{\partial u}{\partial z}\right)\right) + 
\frac{\partial}{\partial y}\left(\mu\left(\frac{\partial w}{\partial y} + \frac{\partial v}{\partial z}\right)\right) + 
\frac{\partial}{\partial z}\left(2 \mu \frac{\partial w}{\partial z} - \frac{2\mu}{3} \nabla \cdot \mathbf{v}\right) + 
\rho g_z</math>

or, more compactly in vector form,

:<math>\rho \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = 
-\nabla p + \nabla \cdot \left(\mu (\nabla \mathbf{v} + (\nabla \mathbf{v})^T)\right) + \nabla \left( - \frac{2\mu}{3}\nabla \cdot \mathbf{v}\right)
+ \rho \mathbf{g}</math>

where the [[matrix transpose|transpose]] has been used. Gravity has been accounted for as "the" body force, i.e. <math>\mathbf{f} = \rho \mathbf{g}</math>. The associated mass continuity equation is:

:<math>\frac{\partial \rho}{\partial t} + \nabla \cdot (\rho \mathbf{v}) = 0</math>

In addition to this equation, an [[equation of state]] and an equation for the conservation of energy is needed. The equation of state to use depends on context (often the [[ideal gas law]]), the conservation of energy will read:

:<math>\rho \frac{D h}{D t} = \frac{D p}{D t} + \nabla \cdot (k \nabla T) + \Phi</math>

Here, <math>h</math> is the [[enthalpy]], <math>T</math> is the [[temperature]], and <math>\Phi</math> is a function representing the dissipation of energy due to viscous effects:

:<math>\Phi = \mu \left(2\left(\frac{\partial u}{\partial x}\right)^2 + 2\left(\frac{\partial v}{\partial y}\right)^2 + 2\left(\frac{\partial w}{\partial z}\right)^2 + \left(\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y}\right)^2 + \left(\frac{\partial w}{\partial y} + \frac{\partial v}{\partial z}\right)^2 + \left(\frac{\partial u}{\partial z} + \frac{\partial w}{\partial x}\right)^2\right) + \lambda (\nabla \cdot \mathbf{v})^2</math>

With a good equation of state and good functions for the dependence of parameters (such as viscosity) on the variables, this system of equations seems to properly model the dynamics of all known gases and most liquids.

====Incompressible Newtonian fluid====
For the special (but very common) case of incompressible flow, the momentum equations simplify significantly. Taking into account the following assumptions:
*Viscosity <math>\mu</math> will now be a constant
*The second viscosity effect <math>\lambda = 0</math>
*The simplified mass continuity equation <math>\nabla \cdot \mathbf{v} = 0</math>
then looking at the viscous terms of the <math>x</math> momentum equation for example we have:
:<math>\begin{align}
 &\frac{\partial}{\partial x}\left(2 \mu \frac{\partial u}{\partial x} + \lambda \nabla \cdot \mathbf{v}\right) + 
 \frac{\partial}{\partial y}\left(\mu\left(\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}\right)\right) + 
 \frac{\partial}{\partial z}\left(\mu\left(\frac{\partial u}{\partial z} + \frac{\partial w}{\partial x}\right)\right) \\  \\
 & = 
 2 \mu \frac{\partial^2 u}{\partial x^2} + 
 \mu \frac{\partial^2 u}{\partial y^2} + \mu \frac{\partial^2 v}{\partial y \, \partial x} + 
 \mu \frac{\partial^2 u}{\partial z^2} + \mu \frac{\partial^2 w}{\partial z \, \partial x} \\  \\
 & = 
 \mu \frac{\partial^2 u}{\partial x^2} + 
 \mu \frac{\partial^2 u}{\partial y^2} + 
 \mu \frac{\partial^2 u}{\partial z^2} + 
 \mu \frac{\partial^2 u}{\partial x^2} + \mu \frac{\partial^2 v}{\partial y \, \partial x} + \mu \frac{\partial^2 w}{\partial z \,  \partial x} \\  \\
 & = \mu \nabla^2 u + \mu \frac{\partial}{\partial x} \cancelto{0}{\left(\frac{\partial u}{\partial x} + \frac{\partial v}{\partial y} + \frac{\partial w}{\partial z}\right)} = \mu \nabla^2 u
\end{align}</math>
Similarly for the <math>y</math> and <math>z</math> momentum directions we have <math>\mu \nabla^2 v</math> and <math>\mu \nabla^2 w</math>.

===Non-Newtonian fluids===
{{main|Non-Newtonian fluid}}
A non-Newtonian fluid is a [[fluid]] whose flow properties differ in any way from those of [[Newtonian fluids]]. Most commonly the [[viscosity]] of non-Newtonian fluids is not independent of [[shear rate]] or shear rate history. However, there are some non-Newtonian fluids with shear-independent viscosity, that nonetheless exhibit normal stress-differences or other non-Newtonian behaviour. Many [[salt]] solutions and molten [[polymers]] are non-Newtonian fluids, as are many commonly found substances such as [[ketchup]], [[custard]], [[toothpaste]], starch suspensions, [[paint]], [[blood]], and [[shampoo]]. In a Newtonian fluid, the relation between the [[shear stress]] and the [[shear rate]] is linear, passing through the origin, the constant of proportionality being the coefficient of viscosity. In a non-Newtonian fluid, the relation between the shear stress and the shear rate is different, and can even be time-dependent. The study of the non-Newtonian fluids is usually called [[rheology]]. A few examples are given here.

===Bingham fluid===
{{main|Bingham plastic}}

In Bingham fluids, the situation is slightly different:

:<math>
\frac{\partial u}{\partial y} = \left\{
\begin{matrix} 
0 &, \quad \tau < \tau_0 \\ 
(\tau - \tau_0)/ {\mu} &, \quad \tau \ge \tau_0 
\end{matrix}\right.</math>

These are fluids capable of bearing some shear before they start flowing. Some common examples are [[toothpaste]] and [[clay]].

===Power-law fluid===
{{main|Power-law fluid}}

A power law fluid is an idealised [[fluid]] for which the [[shear stress]], <math>\tau</math>, is given by

:<math>\tau = K \left(\frac{\partial u}{\partial y}\right)^n </math>

This form is useful for approximating all sorts of general fluids, including shear thinning (such as latex paint) and shear thickening (such as corn starch water mixture).

==Stream function formulation==

In the analysis of a flow, it is often desirable to reduce the number of equations or the number of variables being dealt with, or both. The incompressible Navier-Stokes equation with mass continuity (four equations in four unknowns) can, in fact, be reduced to a single equation with a single dependent variable in 2D, or one vector equation in 3D. This is enabled by two [[vector calculus identities#Combinations of multiple operators|vector calculus identities]]:

:<math>\nabla \times (\nabla \phi) = 0</math>

:<math>\nabla \cdot (\nabla \times \mathbf{A}) = 0</math>

for any differentiable scalar <math>\phi</math> and vector <math>\mathbf{A}</math>. The first identity implies that any term in the Navier-Stokes equation that may be represented as the gradient of a scalar will disappear when the [[Curl (mathematics)|curl]] of the equation is taken. Commonly, pressure and gravity are what eliminate, resulting in (this is true in 2D as well as 3D):

:<math>\nabla \times \left(\frac{\partial \mathbf{v}}{\partial t} + \mathbf{v} \cdot \nabla \mathbf{v}\right) = \nu \nabla \times (\nabla^2 \mathbf{v})</math>

where it's assumed that all body forces are describable as gradients (true for gravity), and density has been divided so that viscosity becomes [[kinematic viscosity]].

The second vector calculus identity above states that the divergence of the curl of a vector field is zero. Since the (incompressible) mass continuity equation specifies the divergence of velocity being zero, we can replace the velocity with the curl of some vector <math>\vec \psi</math> so that mass continuity is always satisfied:

:<math>\nabla \cdot \mathbf v = 0 \quad \Rightarrow \quad \nabla \cdot (\nabla \times \vec \psi) = 0 \quad \Rightarrow \quad 0 = 0</math>

So, as long as velocity is represented through <math>\mathbf v = \nabla \times \vec \psi</math>, mass continuity is unconditionally satisfied. With this new dependent vector variable, the Navier-Stokes equation (with curl taken as above) becomes a single fourth order vector equation, no longer containing the unknown pressure variable and no longer dependent on a separate mass continuity equation:

:<math>\nabla \times \left(\frac{\partial}{\partial t}(\nabla \times \vec \psi) + (\nabla \times \vec \psi) \cdot \nabla (\nabla \times \vec \psi)\right) = \nu \nabla \times (\nabla^2 (\nabla \times \vec \psi))</math>

Apart from containing fourth order derivatives, this equation is fairly complicated, and is thus uncommon. Note that if the cross differentiation is left out, the result is a third order vector equation containing an unknown vector field (the gradient of pressure) that may be determined from the same boundary conditions that one would apply to the fourth order equation above.

===2D flow in orthogonal coordinates===

The true utility of this formulation is seen when the flow is two dimensional in nature and the equation is written in a general [[orthogonal coordinates|orthogonal coordinate system]], in other words a system where the basis vectors are orthogonal. Note that this by no means limits application to [[Cartesian coordinates]], in fact most of the common coordinates systems are orthogonal, including familiar ones like [[Cylindrical coordinates|cylindrical]] and obscure ones like [[Toroidal coordinates|toroidal]].

The 3D velocity is expressed as (note that the discussion has been coordinate free up till now):

:<math>\mathbf v = v_1 \mathbf e_1 + v_2 \mathbf e_2 + v_3 \mathbf e_3</math>

where <math>\mathbf e_i</math> are basis vectors, not necessarily constant and not necessarily normalized, and <math>v_i</math> are velocity components; let also the coordinates of space be <math>(x_1, x_2, x_3)</math>.

Now suppose that the flow is 2D. This doesn't mean the flow is in a plane, rather it means that the component of velocity in one direction is zero and the remaining components are independent of the same direction. In that case (take component 3 to be zero):

:<math>\mathbf v = v_1 \mathbf e_1 + v_2 \mathbf e_2</math>

:<math>\frac{\partial v_1}{\partial x_3} = \frac{\partial v_2}{\partial x_3} = 0</math>

The vector function <math>\vec \psi</math> is still defined via:

:<math>\mathbf v = \nabla \times \vec \psi</math>

but this must simplify in some way also since the flow is assumed 2D. If orthogonal coordinates are assumed, the [[Curl (mathematics)|curl]] takes on a fairly simple form, and the equation above expanded becomes:

:<math>v_1 \mathbf e_1 + v_2 \mathbf e_2 = 
\frac{\mathbf{e}_{1}}{h_{2} h_{3}} 
\left[
\frac{\partial}{\partial x_{2}} \left( h_{3} \psi_{3} \right) - 
\frac{\partial}{\partial x_{3}} \left( h_{2} \psi_{2} \right)
\right] + 
\frac{\mathbf{e}_{2}}{h_{3} h_{1}} 
\left[
\frac{\partial}{\partial x_{3}} \left( h_{1} \psi_{1} \right) - 
\frac{\partial}{\partial x_{1}} \left( h_{3} \psi_{3} \right)
\right] + 
\frac{\mathbf{e}_{3}}{h_{1} h_{2}} 
\left[
\frac{\partial}{\partial x_{1}} \left( h_{2} \psi_{2} \right) - 
\frac{\partial}{\partial x_{2}} \left( h_{1} \psi_{1} \right)
\right]
</math>

Examining this equation shows that we can set <math>\psi_1 = \psi_2 = 0</math> and retain equality with no loss of generality, so that:

:<math>v_1 \mathbf e_1 + v_2 \mathbf e_2 = 
\frac{\mathbf{e}_{1}}{h_{2} h_{3}} \frac{\partial}{\partial x_{2}} \left( h_{3} \psi_{3} \right)
- \frac{\mathbf{e}_{2}}{h_{3} h_{1}} \frac{\partial}{\partial x_{1}} \left( h_{3} \psi_{3} \right)
</math>

the significance here is that only one component of <math>\vec \psi</math> remains, so that 2D flow becomes a problem with only one dependent variable. The cross differentiated Navier–Stokes equation becomes two 0 = 0 equations and one meaningful equation.

The remaining component <math>\psi_3 = \psi</math> is called the '''[[stream function]]'''. The equation for <math>\psi</math> can simplify since a variety of quantities will now equal zero, for example:

:<math>
\nabla \cdot \vec \psi = \frac{1}{h_{1} h_{2} h_{3}} \frac{\partial}{\partial x_3} \left(\psi h_1 h_2\right) = 0
</math>

if the [[orthogonal coordinates|scale factors]] <math>h_1</math> and <math>h_2</math> also are independent of <math>x_3</math>. Also, from the definition of the [[vector Laplacian]]

:<math>
\nabla \times (\nabla \times \vec \psi) = \nabla(\nabla \cdot \vec \psi) - \nabla^2 \vec \psi = -\nabla^2 \vec \psi
</math>

Manipulating the cross differentiated Navier–Stokes equation using the above two equations and a variety of identities<ref>{{cite web
| url = http://mathworld.wolfram.com/VectorDerivative.html
| title = Vector Derivative
| author = [[Eric W. Weisstein]]
| publisher = [[MathWorld]]
| accessdate = 2008-06-07
}}</ref> will eventually yield the 1D scalar equation for the stream function:

:<math>
\frac{\partial}{\partial t}(\nabla^2 \psi)
 + (\nabla \times \vec \psi) \cdot \nabla(\nabla^2 \psi) = \nu \nabla^4 \psi</math>

where <math>\nabla^4</math> is the [[biharmonic operator]]. This is very useful because it is a single self-contained scalar equation that describes both momentum and mass conservation in 2D. The only other equations that this [[partial differential equation]] needs are initial and boundary conditions.

:{| class="toccolours collapsible collapsed" width="60%" style="text-align:left"
!Derivation of the scalar stream function equation
|-
|
Distributing the [[curl (mathematics)|curl]]:

:<math>
\frac{\partial}{\partial t}(\nabla \times(\nabla \times \vec \psi))
 + \nabla \times\left((\nabla \times \vec \psi) \cdot \nabla (\nabla \times \vec \psi)\right)
 = \nu \nabla \times (\nabla^2 (\nabla \times \vec \psi))</math>

Replacing curl of the curl with the Laplacian and expanding convection and viscosity:

:<math>
-\frac{\partial}{\partial t}(\nabla^2 \vec \psi)
 + \nabla \times\left(
\nabla \left( \frac{(\nabla \times \vec \psi) \cdot (\nabla \times \vec \psi)}{2} \right) + \left(\nabla \times (\nabla \times \vec \psi) \right) \times (\nabla \times \vec \psi)
\right)
 = \nu (\nabla^2(\nabla(\nabla \cdot \vec \psi)) - \nabla^4 \vec \psi)</math>

Above, the curl of a gradient is zero, and the divergence of <math>\vec \psi</math> is zero. Negating:

:<math>
\frac{\partial}{\partial t}(\nabla^2 \vec \psi)
 + \nabla \times \left(\nabla^2 \vec \psi \times (\nabla \times \vec \psi)\right)
 = \nu \nabla^4 \vec \psi</math>

Expanding the curl of the cross product into four terms:

:<math>
\frac{\partial}{\partial t}(\nabla^2 \vec \psi) + 
(\nabla \times \vec \psi) \cdot \nabla(\nabla^2 \vec \psi)
- (\nabla^2 \vec \psi) \cdot \nabla (\nabla \times \vec \psi)
+ (\nabla^2 \vec \psi)(\nabla \cdot (\nabla \times \vec \psi))
- (\nabla \times \vec \psi)(\nabla \cdot (\nabla^2 \vec \psi))
 = \nu \nabla^4 \vec \psi</math>

Only one of four terms of the expanded curl is nonzero. The second is zero because it is the dot product of orthogonal vectors, the third is zero because it contains the divergence of velocity, and the fourth is zero because the divergence of a vector with only component three is zero (since it's assumed that nothing (except maybe <math>h_3</math>) depends on component three).

:<math>
\frac{\partial}{\partial t}(\nabla^2 \vec \psi) + (\nabla \times \vec \psi) \cdot \nabla(\nabla^2 \vec \psi)
 = \nu \nabla^4 \vec \psi</math>

This vector equation is one meaningful scalar equation and two 0 = 0 equations.
|}

The assumptions for the stream function equation are listed below:

* The flow is incompressible and Newtonian.
* Coordinates are [[orthogonal coordinates|orthogonal]].
* Flow is 2D: <math>v_3 = \frac{\partial v_1}{\partial x_3} = \frac{\partial v_2}{\partial x_3} = 0</math>
* The first two scale factors of the coordinate system are independent of the last coordinate: <math>\frac{\partial h_1}{\partial x_3} = \frac{\partial h_2}{\partial x_3} = 0</math>, otherwise extra terms appear.

The [[stream function]] has some useful properties:

* Since <math>-\nabla^2 \vec \psi = \nabla \times (\nabla \times \vec \psi) = \nabla \times \mathbf v</math>, the [[vorticity]] of the flow is just the negative of the Laplacian of the stream function.
* The [[level curves]] of the stream function are [[Streamlines, streaklines, and pathlines|streamlines]].

==The stress tensor==
The derivation of the Navier-Stokes equation involves the consideration of forces acting on fluid elements, so that a quantity called the [[stress (physics)|stress tensor]] appears naturally in the [[Cauchy momentum equation]]. Since the divergence of this tensor is taken, it is customary to write out the equation fully simplified, so that the original appearance of the stress tensor is lost.

However, the stress tensor still has some important uses, especially in formulating boundary conditions at [[Capillary surface|fluid interfaces]]. Recalling that <math>\sigma = -\pi I + \mathbb{T}</math>, for a Newtonian fluid the stress tensor is:

:<math>
\sigma_{ij} = -p\delta_{ij}+ \mu\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right) + \delta_{ij} \lambda \nabla \cdot \mathbf{v}.
</math>

If the fluid is assumed to be incompressible, the tensor simplifies significantly:

:<math>
\begin{align}

\sigma &= 
-\begin{pmatrix}
p&0&0\\
0&p&0\\
0&0&p
\end{pmatrix} +

\mu \begin{pmatrix}
2 \displaystyle{\frac{\partial u}{\partial x}} & \displaystyle{\frac{\partial u}{\partial y} + \frac{\partial v}{\partial x}} &\displaystyle{ \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x}} \\
\displaystyle{\frac{\partial v}{\partial x} + \frac{\partial u}{\partial y}} & 2 \displaystyle{\frac{\partial v}{\partial y}} & \displaystyle{\frac{\partial v}{\partial z} + \frac{\partial w}{\partial y}} \\
\displaystyle{\frac{\partial w}{\partial x} + \frac{\partial u}{\partial z}} & \displaystyle{\frac{\partial w}{\partial y} + \frac{\partial v}{\partial z}} & 2\displaystyle{\frac{\partial w}{\partial z}}
\end{pmatrix} \\

&= -p I + \mu (\nabla \mathbf{v} + (\nabla \mathbf{v})^T) = -p I + 2 \mu e\\

\end{align}
</math>

<math>e</math> is the [[strain rate]] tensor, by definition:

:<math>e_{ij} = \frac{1}{2}\left(\frac{\partial u_i}{\partial x_j} + \frac{\partial u_j}{\partial x_i}\right).</math>

==Notes==
{{Reflist}}

==References==
*{{Cite book
 | first = G.K.
 | last = Batchelor
 | authorlink = George Batchelor
 | title = An Introduction to Fluid Dynamics
 | year = 2000
 | location = New York
 | publisher = Cambridge University Press
 | isbn = 978-0-521-66396-0
 | postscript = <!-- Bot inserted parameter. Either remove it; or change its value to "." for the cite to end in a ".", as necessary. -->
}}
*{{Cite book
  | last = White
  | first = Frank M.
  | title = Viscous Fluid Flow
  | edition = 3rd
  | publisher = McGraw Hill
  | year = 2006
  | location = New York, NY
  | pages =
  | isbn = 0-07-240231-8
  | postscript = <!-- Bot inserted parameter. Either remove it; or change its value to "." for the cite to end in a ".", as necessary. -->
}}
*[http://web.mit.edu/1.63/www/Lec-notes/Surfacetension/ Surface Tension Module], by John W. M. Bush, at [[MIT OCW]].
*[http://www.hrpub.org/download/20131107/UJAM1-12600416.pdf/ Navier –Stokes First Exact Transformation]http://www.hrpub.org/download/20131107/UJAM1-12600416.pdf

{{DEFAULTSORT:Derivation Of The Navier-Stokes Equations}}
[[Category:Articles with inconsistent citation formats]]
[[Category:Concepts in physics]]
[[Category:Equations of fluid dynamics]]
[[Category:Aerodynamics]]
[[Category:Partial differential equations]]