Navier-stokes Equations

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.

Basic assumptions

NS equations assume that the fluid is a continuum. They are conservation equations over a control volume which we will call \Omega, which keeps the portion of the fluid which is studied. Unfortunately, \Omega 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 \Omega 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 \Omega, 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:

Equation of continuity

Conservation of mass is writen:
\frac{D\rho}{D t} + \rho\nabla\mathbf{v}= 0
with \rho the density of the fluid. If the fluid is incompressible, \rho is not a function of time. As \rho is non-nil in general, in the case of an incompressible fluid, 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
with p 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 \tau 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 \tau 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 hase 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 \mu 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}\left2\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.
  • \rho is the density,
  • u_i (i=1,2,3) the three components of velocity,
  • f_i body forces (such as gravity),
  • p the pressure,
  • \mu 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);
  • \Delta=e_{ii} is the divergence,
  • \delta_{ij} is the Kronecker delta.
If \mu 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 \mu=0 the resulting equations are known as the Euler equations; there, the emphasis is on compressible flow and shock waves).
If now in addition \rho 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 \leftv_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 \leftv_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 \leftv_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 editon 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. 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.

See also

References

  • Inge L. Rhyming Dynamique des fluides, 1991 PPUR

External links

 

<< PreviousWord BrowserNext >>
le ptomane
leonard rossiter
textile
substitution cipher
transposition cipher
geographic coordinate system
rot13
software architecture
polyurethane
westchester county, new york
nassau county, new york
jean pierre abbat
oscar niemeyer
ct
triple des
346 bc
paul ricard circuit
celestial coordinate system
eudoxus
eric
equatorial coordinate system
horizontal coordinate system
ecliptic coordinate system
eudoxus of cnidus
galactic coordinate system
phenytoin
mathematical analysis
isocyanate
analysis (philosophy)
internet backbone
ring (mathematics)
caesar cipher
367
368
369
365
364
363
362
361
360
gottlob frege
359
lucas cranach