BIGpedia - Navier-Stokes equations - Encyclopedia and Dictionary Online

In fluid dynamics, the Navier-Stokes equations, named after Claude-Louis Navier and George Gabriel Stokes are a set of nonlinear partial differential equations that describe the flow of fluids such as liquids and gases. For example, they model weather or the movement of air in the atmosphere, ocean currents, water flow in a pipe, as well as many other fluid flow phenomena.

It is a famous open question whether smooth initial conditions always lead to smooth solutions for all times; a $1,000,000 prize was offered in May 2000 by the Clay Mathematics Institute for the answer to this question.

Contents

1 Basic assumptions

1.1 The substantive derivative
1.2 Conservation laws

1.2.1 Equation of continuity
1.2.2 Conservation of momentum

2 The equations

2.1 General form

2.1.1 The form of the equations
2.1.2 The closure problem

2.2 Special forms

2.2.1 Newtonian fluids
2.2.2 Bingham fluids
2.2.3 Incompressible fluids

3 See also

4 References

5 External links

Basic assumptions

NS equations assume that the fluid is a continuum. They are conservation equations over a control volume which we will call $Ω$ , which keeps the portion of the fluid which is studied. Unfortunately, $Ω$ depends on time, and the conservation laws need then to be modified so as to be valid on arbitrary control volumes.

The substantive derivative

main article: substantive derivative.

Before going into the details of the Navier-Stokes equations, one must first define an operator:

$\frac{D}{Dt}(\cdot)$

This is the usual time derivative, but when following a particle of the fluid. The full form is:

$\frac{\partial(\cdot)}{\partial t} + \mathbf{v}\nabla(\cdot)$

The integral form of a conservation law L over the control volume $Ω$ is:

$\frac{d}{dt}\int_\Omega \mathbf{L}\; d\Omega = 0$

Using the particular derivative, one can swap the $\frac{d}{dt}$ and $\int$ operators:

$\int_\Omega \frac{D}{Dt}\mathbf{L}\; d\Omega = 0\quad \forall \Omega$

As this expression is valid for all $Ω$ , it simplifies to:

$\frac{D}{Dt}\mathbf{L} = 0$

Conservation laws

Main article: conservation laws.

the NS equations are derived from conservation laws using the transformation described above and need to be closed using state laws.

Over a control volume, using the transformation described above, the following quantities are deemed conserved:

Mass
Energy
Momentum
Moment of momentum aka angular momentum

Equation of continuity

Conservation of mass is written:

$\frac{D\rho}{D t} + \rho\nabla\mathbf{v}= 0$

where

ρ

is the density of the fluid.

In the case of an incompressible fluid $ρ$ is not a function of time; the equation is reduced to:

$\nabla\mathbf{v} = 0$

Conservation of momentum

Conservation of momentum is written:

$\frac{\partial\rho\mathbf{v}}{\partial t} + \nabla(\rho\mathbf{v}\otimes\mathbf{v}) = \sum\mathbf{F}$

Note that $\mathbf{v}\otimes\mathbf{v}$ is a tensor, the $\otimes$ representing the tensor product.

We can simplify it further, this becomes:

$\rho\frac{D\mathbf{v}}{D t}=\sum\mathbf{F}$

In which we recognise the usual F=ma.

The equations

General form

The form of the equations

The general form of the Navier-Stokes equations is:

$\rho\frac{D\mathbf{v}}{D t} = \nabla \mathbb{P} + \rho\mathbf{f}$

For the conservation of momentum. The tensor $\mathbb{P}$ represents the surface forces applied on a fluid particle. In general, we have the form:

$\nabla \mathbb{P} = \begin{pmatrix} \frac{\partial \sigma_{xx}}{\partial x} & \frac{\partial\tau_{yx}}{\partial y} & \frac{\partial\tau_{zx}}{\partial z} \\ \\ \frac{\partial\tau_{xy}}{\partial x} & \frac{\partial \sigma_{yy}}{\partial y} & \frac{\partial\tau_{zy}}{\partial z} \\ \\ \frac{\partial\tau_{xz}}{\partial x} & \frac{\partial\tau_{yz}}{\partial y} & \frac{\partial\sigma_{zz}}{\partial z}\end{pmatrix}$

To which we add the continuity equation:

$\frac{D\rho}{Dt} + \rho\nabla\mathbf{v} = 0$

The nature of the diagonal of $\mathbb{P}$ is known, it is the gradient of pressure, thus:

$\nabla \mathbf{\sigma} = \nabla p$

where

p is the pressure

finally, we have:

$\rho\frac{D\mathbf{v}}{D t} = \nabla p + \nabla \mathbb{T} + \rho\mathbf{f}$

where the components of $\mathbb{T}$ are the $τ$ of $\mathbb{P}$ .

The closure problem

Those equations are incomplete. To complete them, one must make hypotheses on the form of $\mathbb{P}$ . In the case of an perfect fluid $τ$ components are nil, for example.

The components of $\mathbb{P}$ are the constraints on an infinitesimal element of fluid. They represent the normal and shear constraints. $\mathbb{P}$ is therefore by definition antisymmetric.

So-called non-Newtonian fluids are simply fluids where this tensor has no special properties allowing for special solutions of the equations.

Special forms

Those are certain usual simplifications of the problem for which sometimes solutions are known.

Newtonian fluids

In Newtonian fluids the following assumption holds:

$\tau_{ij}=\mu\frac{\partial v_i}{\partial x_j}$

where

μ

is the viscosity of the fluid.

Bingham fluids

In Bingham fluids, we have something slightly different:

$\tau_{ij}=\tau_0 + \mu\frac{\partial v_i}{\partial x_j},\;\frac{\partial v_i}{\partial x_j}>0$

Those are fluids capable of bearing some shear before they start flowing. An example is tooth paste.

Incompressible fluids

The Navier-Stokes equations are

$\rho\frac{Du_i}{Dt}=\rho f_i-\frac{\partial p}{\partial x_i}+\frac{\partial}{\partial x_j}\left[ 2\mu\left(e_{ij}-\Delta\delta_{ij}/3\right)\right]$

for momentum conservation and

$\frac{\partial\rho u_i}{\partial x_i}=0$

for conservation of mass.

where

ρ

is the density,

u i

(

i = 1,2,3

) the three components of velocity,

f i

body forces (such as gravity),

p

the pressure,

μ

the dynamic viscosity, of the fluid at a point;

$e_{ij}=\left(\frac{\partial u_i}{\partial x_j}+\frac{\partial u_j}{\partial x_i}\right)$ ;

Δ = e i i

is the divergence,

δ i j

is the Kronecker delta.

If $μ$ is uniform over the fluid, the momentum equation above simplifies to

$\rho\frac{Du_i}{Dt}=\rho f_i-\frac{\partial p}{\partial x_i} +\mu \left( \frac{\partial^2u_i}{\partial x_i\partial x_j}+ \frac{1}{3}\frac{\partial\Delta}{\partial x_i}\right)$

(if $μ = 0$ the resulting equations are known as the Euler equations; there, the emphasis is on compressible flow and shock waves).

If now in addition $ρ$ is assumed to be constant:

$\rho \left({\partial v_x \over \partial t}+ v_x {\partial v_x \over \partial x}+ v_y {\partial v_x \over \partial y}+ v_z {\partial v_x \over \partial z}\right)= \mu \left[{\partial v_x^2 \over \partial x^2}+{\partial v_x^2 \over \partial y^2}+{\partial v_x^2 \over \partial z^2}\right]-{\partial p \over \partial x} +\rho g_x$

$\rho \left({\partial v_y \over \partial t}+ v_x {\partial v_y \over \partial x}+ v_y {\partial v_y \over \partial y}+ v_z {\partial v_y \over \partial z}\right)= \mu \left[{\partial v_y^2 \over \partial x^2}+{\partial v_y^2 \over \partial y^2}+{\partial v_y^2 \over \partial z^2}\right]-{\partial p \over \partial y} +\rho g_y$

$\rho \left({\partial v_z \over \partial t}+ v_x {\partial v_z \over \partial x}+ v_y {\partial v_z \over \partial y}+ v_z {\partial v_z \over \partial z}\right)= \mu \left[{\partial v_z^2 \over \partial x^2}+{\partial v_z^2 \over \partial y^2}+{\partial v_z^2 \over \partial z^2}\right]-{\partial p \over \partial z} +\rho g_z$

Continuity equation (assuming incompressibility):

${\partial v_x \over \partial x}+{\partial v_y \over \partial y}+{\partial v_z \over \partial z}=0$

Simplified version of the N-S equations. Adapted from Incompressible Flow, second edition by Ronald Panton

The equations are derived by considering the mass, momentum, and energy balances for an infinitesimal control volume. The Navier-Stokes equations need to be augmented by an equation of state for compressible flows. The variables to be solved for are the velocity components, the fluid density, static pressure, and temperature. The flow is assumed to be differentiable and continuous, allowing these balances to be expressed as partial differential equations. The equations can be converted to Wilkinson equations for the secondary variables vorticity and stream function. Solution depends on the fluid properties (such as viscosity, specific heats, and thermal conductivity), and on the boundary conditions of the domain of study. For a derivation of the Navier-Stokes equations, see some of the external links listed below.

Note that the Navier-Stokes equations can only describe fluid flow approximately and that, at very small scales or under extreme conditions, real fluids made out of mixtures of discrete molecules and other material, such as suspended particles and dissolved gases, will produce different results from the continuous and homogeneous fluids modelled by the Navier-Stokes equations. Depending on the Knudsen number of the problem, statistical mechanics may be a more appropriate approach. However, the Navier-Stokes equations are useful for a wide range of practical problems, providing their limitations are borne in mind.

Although the full, unsteady Navier-Stokes equations correctly describe nearly all flows of practical interest, they are too complex for practical solution in many cases and a special "reduced" form of the full equations is often used instead — these are the Reynolds-averaged Navier-Stokes (RANS) equations. The solution of the full steady Navier-Stokes equations is sufficiently accurate alone for cases where the fluid flow is laminar. For turbulent flows the Reynolds-averaged form of the equations are most commonly used. The RANS form of the equations introduce new terms that reflect the additional modelling of the small turbulent motions.

Solution of flow equations by numerical methods is called computational fluid dynamics.

References

Inge L. Rhyming Dynamique des fluides, 1991 PPUR

External links

The Clay Mathematics Institute's Navier-Stokes equation prize
- Official statement of the problem
A good derivation of the Navier-Stokes equations
Derivation of the Navier-Stokes equations
NASA page on the Navier-Stokes equations
Navier-Stokes equations (exact solutions) at EqWorld: The World of Mathematical Equations

Categories: Fluid dynamics | Partial differential equations | Equations | Unsolved problems in mathematics

The contents of this article are licensed from Wikipedia.org under the GNU Free Documentation License.
How to see transparent copy

01-04-2007 01:21:04