Navier-Stokes Equation

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.

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 modeling 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.

Basic assumptions

The equations are derived by considering the mass, momentum, and energy balances for an infinitesimal control volume.

Indeed 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}; dOmega = 0$

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

$int_Omega frac{D}{Dt}mathbf{L}; dOmega = 0quad 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{Drho}{D t} + rhonablamathbf{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:

$nablamathbf{v} = 0$

Conservation of momentum

Conservation of momentum is written:

$frac{partialrhomathbf{v}}{partial t} + nabla(rhomathbf{v}otimesmathbf{v}) = summathbf{F}$

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

We can simplify it further, this becomes:

$rhofrac{Dmathbf{v}}{D t}=summathbf{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:

$rhofrac{Dmathbf{v}}{D t} = nabla mathbb{P} + rhomathbf{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{partialtau_{yx}}{partial y} & frac{partialtau_{zx}}{partial z} frac{partialtau_{xy}}{partial x} & frac{partial sigma_{yy}}{partial y} & frac{partialtau_{zy}}{partial z} frac{partialtau_{xz}}{partial x} & frac{partialtau_{yz}}{partial y} & frac{partialsigma_{zz}}{partial z}end{pmatrix}$

To which we add the continuity equation:

$frac{Drho}{Dt} + rhonablamathbf{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:

$rhofrac{Dmathbf{v}}{D t} = nabla p + nabla mathbb{T} + rhomathbf{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 a perfect fluid $τ$ components are nil, for example. Those equations used to complete the set are equations of state.

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.

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

Main article Newtonian fluids.

In Newtonian fluids the following assumption holds:

$tau_{ij}=mufrac{partial v_i}{partial x_j}$

where

μ

is the viscosity of the fluid.

Bingham fluids

Main article Bingham plastic.

In Bingham fluids, we have something slightly different:

$tau_{ij}=tau_0 + mufrac{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.

Power-law fluid

Main article Power-law fluid.

It is an idealised fluid for which the shear stress, $τ$ , is given by

$tau = K left( frac {partial u} {partial y} right)^n$

This form is useful for approximating all sorts of general fluids.

Incompressible fluids

Main article Incompressible fluids.

The Navier-Stokes equations are

$rhofrac{Du_i}{Dt}=rho f_i-frac{partial p}{partial x_i}+frac{partial}{partial x_j}left[ 2muleft(e_{ij}-Deltadelta_{ij}/3right)right]$

for momentum conservation and

$frac{partialrho 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

$rhofrac{Du_i}{Dt}=rho f_i-frac{partial p}{partial x_i} +mu left( frac{partial^2u_i}{partial x_ipartial x_j}+ frac{1}{3}frac{partialDelta}{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 we obtain the following system:

$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^2 v_x over partial x^2}+{partial^2 v_x over partial y^2}+{partial^2 v_x 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^2 v_y over partial x^2}+{partial^2 v_y over partial y^2}+{partial^2 v_y 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^2 v_z over partial x^2}+{partial^2 v_z over partial y^2}+{partial^2 v_z 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

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.

출처 : wikipedia

기억보관소

Navier-Stokes Equation