Navier-Stokes Equations

ChessisGood · 2012-02-16T16:00:33-08:00

Navier-Stokes Equations: Greetings everyone, Today I would like to offer an interesting puzzle in the fields of mathematics, for discussion and attempted solving. A few days ago, I came across the website of the Clay Mathematics Institute. There, I found an interesting page listing seven so-called &

Sort:

ChessisGood

Feb 16, 2012

#61

There are two types of forces: body(mass) forces and surface forces. Body forces act on the entire control volume. The most common body force is that due to gravity. Electromagnetic phenomena may also create body forces, but this is a rather specialized situation.

Surface forces act on only surface of a control volume at a time and arise due to pressure or viscous stresses.

We find a general expression for the surface force per unit volume of a deformable body. Consider a rectangular parallelepiped with sides $dx,dy,dz$ and hence with volume $dV = dxdydz$

ChessisGood

Feb 16, 2012

#62

ChessisGood

Feb 16, 2012

#63

At the moment we assume this parallelepiped isolated from the rest of the fluid flow , and consider the forces acting on the faces of the parallelepiped.

Let the left forward top of a parallelepiped lies in a point O

To both faces of the parallelepiped perpendicular to the axis $x$ and having the area $dydz$ applied resulting stresses , equal to $\textbf{p}_{x}$ and $\textbf{p}_{x}+\frac{\partial \textbf{p}_{x}}{ \partial x}$ respectively

ChessisGood

Feb 16, 2012

#64

So we have

for $x$ - direction $\frac{ \partial \textbf{p}_{x}}{ \partial x} dxdydz$

for $y$ - direction $\frac{ \partial \textbf{p}_{y}}{ \partial y} dxdydz$

for $z$ - direction $\frac{ \partial \textbf{p}_{z}}{ \partial z} dxdydz$

$\textbf{P} = \frac{\partial \textbf{p}_{x}}{ \partial x} + \frac{\partial \textbf{p}_{y}}{ \partial y} + \frac{\partial \textbf{p}_{z}}{ \partial z}$

(6)

$\rho \frac{d \textbf{V}}{dt} = \rho \textbf{F} + \frac{\partial \textbf{p}_{x} }{ \partial x} + \frac{\partial \textbf{p}_{y} }{ \partial y} + \frac{\partial \textbf{p}_{z} }{ \partial z}$

(7)

$\left. \begin{array}{c} \rho \frac{du}{dt}= \rho F_{x} + \frac{ \partial p_{xx}}{\partial x} + \frac{ \partial p_{yx}}{\partial y} + \frac{ \partial p_{zx}}{\partial z} \\ \rho \frac{dv}{dt}= \rho F_{y} + \frac{ \partial p_{xy}}{\partial x} + \frac{ \partial p_{yy}}{\partial y} + \frac{ \partial p_{zy}}{\partial z} \\ \rho \frac{dw}{dt}= \rho F_{z} + \frac{ \partial p_{xz}}{\partial x} + \frac{ \partial p_{yz}}{\partial y} + \frac{ \partial p_{zz}}{\partial z} \\ \end{array} \right\}$

(8)

$\left. \begin{array}{c} \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) = \rho F_{x} + \frac{ \partial p_{xx}}{\partial x} + \frac{ \partial p_{yx}}{\partial y} + \frac{ \partial p_{zx}}{\partial z} \\ \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) = \rho F_{y} + \frac{ \partial p_{xx}}{\partial x} + \frac{ \partial p_{yx}}{\partial y} + \frac{ \partial p_{zx}}{\partial z} \\ \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) = \rho F_{z} + \frac{ \partial p_{xx}}{\partial x} + \frac{ \partial p_{yx}}{\partial y} + \frac{ \partial p_{zx}}{\partial z} \\ \end{array} \right\}$

(9)

ChessisGood

Feb 16, 2012

#65

The force due to the stress is the product of the stress and the area over which it acts.

$\textbf{P}_{x} = \textbf{i} \sigma_{xx} + \textbf{j} \tau_{xy} + \textbf{k} \tau_{xz}$

(10)

$\textbf{P}_{y} = \textbf{i} \tau_{yx} + \textbf{j} \sigma_{yy} + \textbf{k} \tau_{yz}$

(11)

$\textbf{P}_{z}=\textbf{i}\tau_{zx} + \textbf{j} \tau_{zy} + \textbf{k}*\sigma_{zz}$

(12)

$\Pi = \left( \begin{array}{ccc} \sigma_{xx} & \tau_{xy} & \tau_{xz} \\ \tau_{yx} & \sigma_{yy} & \tau_{yz} \\ \tau_{zx} & \tau_{zy} & \sigma_{zz} \\ \end{array} \right)$

(13)

$\left. \begin{array}{c} \rho \frac{du}{dt} = \rho F_{x} - \frac{\partial p}{\partial x} + \left( \frac{\partial \sigma_{x} '}{\partial x} + \frac{\partial \tau_{xy}}{\partial y} + \frac{\partial \tau_{xz}}{\partial z} \right) \\ \rho \frac{dv}{dt} = \rho F_{y} - \frac{\partial p}{\partial y} + \left( \frac{\partial \tau_{xy} }{\partial x} + \frac{\partial \sigma_{y} '}{\partial y} + \frac{\partial \tau_{yz}}{\partial z} \right) \\ \rho \frac{dw}{dt} = \rho F_{z} - \frac{\partial p}{\partial z} + \left( \frac{\partial \tau_{xz} '}{\partial x} + \frac{\partial \tau_{yz}}{\partial y} + \frac{\partial \sigma_{z}'}{\partial z} \right) \end{array} \right\}$

(14)

ChessisGood

Feb 16, 2012

#66

deformation and rotation

ChessisGood

Feb 16, 2012

#67

$\left. \begin{array}{c} u = u_{0} + \left( \frac{\partial u}{\partial x} \right)_{0} \left( x - x_{0} \right) + \left( \frac{\partial u}{\partial y} \right)_{0} \left( y - y_{0} \right) + \left( \frac{\partial u}{\partial z} \right)_{0} \left( z - z_{0} \right) \\ v = v_{0} + \left( \frac{\partial v}{\partial x} \right)_{0} \left( x - x_{0} \right) + \left( \frac{\partial v}{\partial y} \right)_{0} \left( y - y_{0} \right) + \left( \frac{\partial v}{\partial z} \right)_{0} \left( z - z_{0} \right) \\ w = w_{0} + \left( \frac{\partial w}{\partial x} \right)_{0} \left( x - x_{0} \right) + \left( \frac{\partial w}{\partial y} \right)_{0} \left( y - y_{0} \right) + \left( \frac{\partial w}{\partial z} \right)_{0} \left( z - z_{0} \right) \\ \end{array} \right\}$

(15)

$\left. \begin{array}{c} u= u_{0} + \omega_{y} \left( z - z_{0} \right) - \omega_{z} \left( y - y_{0} \right) \\ v= v_{0} + \omega_{z} \left( x - x_{0} \right) - \omega_{x} \left( z - z_{0} \right) \\ w= w_{0} + \omega_{x} \left( y - y_{0} \right) - \omega_{y} \left( x - x_{0} \right) \\ \end{array} \right\}$

(16)

$\left. \begin{array}{c} \omega_{x} = \frac{1}{2} \left( \frac{\partial w}{\partial y} - \frac{\partial v}{\partial z} \right) \\ \omega_{y} = \frac{1}{2} \left( \frac{\partial u}{\partial z} - \frac{\partial w}{\partial x} \right) \\ \omega_{z} = \frac{1}{2} \left( \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \right) \\ \end{array} \right\}$

(18)

$\left. \begin{array}{c} u_{solid} = u_{0} + \frac{1}{2} \left( \frac{\partial u}{\partial z} - \frac{\partial w}{\partial x} \right)_{0} \left( z - z_{0} \right) - \frac{1}{2} \left( \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \right)_{0} \left( y - y_{0} \right) \\ v_{solid} = v_{0} + \frac{1}{2} \left( \frac{\partial v}{\partial x} - \frac{\partial u}{\partial y} \right)_{0} \left( x - x_{0} \right) - \frac{1}{2} \left( \frac{\partial w}{\partial y} - \frac{\partial v}{\partial z} \right)_{0} \left( z - z_{0} \right) \\ w_{solid} = w_{0} + \frac{1}{2} \left( \frac{\partial w}{\partial y} - \frac{\partial v}{\partial z} \right)_{0} \left( y - y_{0} \right) - \frac{1}{2} \left( \frac{\partial u}{\partial z} - \frac{\partial w}{\partial x} \right)_{0} \left( x - x_{0} \right) \\ \end{array} \right\}$

(19)

$\left. \begin{array}{c} u = u_{solid} + u_{def} \\ v = v_{solid} + v_{def} \\ w = w_{solid} + w_{def} \\ \end{array} \right\}$

(20)

$\left. \begin{array}{c} u_{def} = \left( \frac{ \partial u}{ \partial x} \right)_{0} \left( x - x_{0} \right) + \frac{1}{2} \left( \frac{ \partial v}{ \partial x} + \frac{ \partial u}{ \partial y} \right)_{0} \left( y - y_{0} \right) + \frac{1}{2} \left( \frac{ \partial u}{ \partial z} + \frac{ \partial w}{ \partial x} \right)_{0} \left( z - z_{0} \right)\\ v_{def} = \frac{1}{2} \left( \frac{ \partial v}{ \partial x} + \frac{ \partial u}{ \partial y} \right)_{0} \left( x - x_{0} \right) + \left( \frac{ \partial v}{ \partial y} \right)_{0} \left( y - y_{0} \right) + \frac{1}{2} \left( \frac{ \partial u}{ \partial z} + \frac{ \partial w}{ \partial x} \right)_{0} \left( z - z_{0} \right)\\ w_{def} = \frac{1}{2} \left( \frac{ \partial u}{ \partial z} + \frac{ \partial w}{ \partial x} \right)_{0} \left( x - x_{0} \right) + \frac{1}{2} \left( \frac{ \partial w}{ \partial y} + \frac{ \partial v}{ \partial z} \right)_{0} \left( y - y_{0} \right) + \left( \frac{ \partial w}{ \partial z} \right)_{0} \left( z - z_{0} \right)\\ \end{array} \right\}$

(21)

$\dot{\epsilon}_{ij} \equiv \left( \begin{array}{ccc} \dot{\epsilon}_{x} & \dot{\epsilon}_{xy} & \dot{\epsilon}_{xz} \\ \dot{\epsilon}_{yx} & \dot{\epsilon}_{y} & \dot{\epsilon}_{yz} \\ \dot{\epsilon}_{zx} & \dot{\epsilon}_{zy} & \dot{\epsilon}_{z} \\ \end{array} \right) \equiv$

(11)

$\equiv \left( \begin{array}{ccc} \frac{\partial u}{ \partial x} & \frac{1}{2} \left( \frac{\partial v}{ \partial x} + \frac{\partial u}{ \partial y} \right) & \frac{1}{2} \left( \frac{\partial w}{ \partial x} + \frac{\partial u}{ \partial z} \right) \\ \frac{1}{2} \left( \frac{\partial u}{ \partial y} + \frac{\partial v}{ \partial x} \right) & \frac{\partial u}{ \partial x} & \frac{1}{2} \left( \frac{\partial v}{ \partial x} + \frac{\partial u}{ \partial y} \right) \\ \frac{1}{2} \left( \frac{\partial u}{ \partial y} + \frac{\partial v}{ \partial x} \right) & \frac{1}{2} \left( \frac{\partial u}{ \partial y} + \frac{\partial v}{ \partial x} \right) & \frac{\partial u}{ \partial x} \end{array} \right)$

(22)

ChessisGood

Feb 16, 2012

#68

ChessisGood

Feb 16, 2012

#69

Newtonian Fluids

Newton came up with the idea of requiring the stress $\tau$ to be linearly proportional to the time rate at which at which strain occurs. Specifically he studied the following problem. There are two flat plates separated by a distance $h$ . The top plate is moved at a velocity $V$ , while the bottom plate is held fixed.

ChessisGood

Feb 16, 2012

#70

ChessisGood

Feb 16, 2012

#71

Newton postulated (since then experimentally verified) that the shear force or shear stress needed to deform the fluid was linearly proportional to the velocity gradient:

$\tau \propto \frac{V}{h}$

(2)

The proportionality factor turned out to be a constant at moderate temperatures, and was called the coefficient of viscosity, $\mu$ . Furthermore, for this particular case, the velocity profile is linear, giving $v / h= \partial u / \partial y$ .

Therefore, Newton postulated:

$\tau = \mu \frac{\partial u}{\partial y}$

(2)

Fluids that have a linear relationship between stress and strain rate are called Newtonian fluids. This is a property of the fluid, not the flow. Water and air are examples of Newtonian fluids, while blood is a non-Newtonian fluid.

ChessisGood

Feb 16, 2012

#72

Stokes Hypothesis

Stokes extended Newton's idea from simple 1-D flows (where only one component of velocity is present) to multidimensional flows. He developed the following relations, collectively known as Stokes relations

$\sigma_{x} = 2 \mu \frac{ \partial u}{ \partial x} + \lambda \left( \frac{ \partial u}{ \partial x} + \frac{ \partial v}{ \partial y} + \frac{ \partial w}{ \partial z} \right)$

(12)

$\sigma_{y} = 2 \mu \frac{ \partial v}{ \partial y} + \lambda \left( \frac{ \partial u}{ \partial x} + \frac{ \partial v}{ \partial y} + \frac{ \partial w}{ \partial z} \right)$

(12)

$\sigma_{z} = 2 \mu \frac{ \partial w}{{\partial}z} + \lambda \left( \frac{ \partial u}{ \partial x} + \frac{ \partial v}{ \partial y} + \frac{ \partial w}{ \partial z} \right)$

(12)

$\tau_{xy} = \tau_{yx} = \mu \left( \frac{\partial u}{\partial y} + \frac{\partial v}{\partial x} \right)$

(12)

$\tau_{xz} = \tau_{zx} = \mu \left( \frac{\partial u}{\partial z} + \frac{\partial w}{\partial x} \right)$

(12)

$\tau_{zy} = \tau_{yz} = \mu \left( \frac{\partial w}{\partial y} + \frac{\partial v}{\partial z} \right)$

(12)

The quantity $\mu$ is called molecular viscosity, and is a function of temperature.

The coefficient $\lambda$ was chosen by Stokes so that the sum of the normal stresses $\sigma_{x}$ , $\sigma_{y}$ and $\sigma_{z}$ are zero.

Then

$\lambda = - \frac{2}{3}\mu$

(12)

ChessisGood

Feb 16, 2012

#73

substitution

$\left. \begin{array}{c} \rho \frac{du}{dt} = \rho F_{x}- \frac{\partial p}{ \partial x} + 2 \frac{\partial}{ \partial x} \left( \mu \frac{ \partial u }{ \partial x } \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] - \frac{2}{3} \frac{\partial}{\partial x}\left( \mu div \textbf{V}\right)\\ \rho \frac{dv}{dt} = \rho F_{y} - \frac{\partial p}{ \partial y} + \frac{\partial }{ \partial x} \left[ \mu \left( \frac{\partial u}{ \partial y} + \frac{\partial v}{ \partial x} \right)\right] + 2 \frac{\partial}{ \partial y} \left( \mu \frac{ \partial v }{ \partial y } \right) + \frac{\partial }{ \partial z} \left[ \mu \left( \frac{\partial v}{ \partial z} + \frac{\partial w}{ \partial y} \right)\right] - \frac{2}{3} \frac{\partial }{ \partial y} \left( \mu div \textbf{V} \right) \\ \rho \frac{dw}{dt} = \rho F_{z} - \frac{\partial p}{ \partial z}+ \frac{\partial }{ \partial x} \left[ \mu \left( \frac{\partial u}{ \partial z} + \frac{\partial w}{ \partial x} \right)\right] + \frac{\partial }{ \partial y} \left[ \mu \left( \frac{\partial v}{ \partial z} + \frac{\partial w}{ \partial y} \right)\right] + 2 \frac{\partial}{ \partial z} \left( \mu \frac{ \partial w }{ \partial z } \right) - \frac{2}{3} \frac{\partial }{ \partial z} ( \mu div \textbf{V} ) \\ \end{array} \right\}$

ChessisGood

Feb 16, 2012

#74

Existence and uniqueness

The existence and uniqueness of classical solutions of the 3-D Navier-Stokes equations is still an open mathematical problem and is one of the Clay Institute's Millenium Problems. In 2-D, existence and uniqueness of regular solutions for all time have been shown by Jean Leray in 1933. He also gave the theory for the existence of weak solutions in the 3-D case while uniqueness is still an open question.

However, recently, Prof. Penny Smith submitted a paper, Immortal Smooth Solution of the Three Space Dimensional Navier-Stokes System, which may provide a proof of the existence and uniqueness.(It has a serious flaw, so the author withdrew the paper)

ChessisGood

Feb 16, 2012

#75

History

Claude Louis Marie Henri Navier’s name is associated with the famous Navier-Stokes equations that govern motion of a viscous fluid. He derived the Navier-Stokes equations in a paper in 1822. His derivation was however based on a molecular theory of attraction and repulsion between neighbouring molecules. Euler had already derived the equations for an ideal fluid in 1755 which did not include the effects of viscosity. Navier did not recognize the physical significance of viscosity and attributed the viscosity coefficient to be a function of molecular spacing.

The equations of motion were rederived by Cauchy in 1828 and by Poisson in 1829. In 1843 Barre de Saint-Venant published a derivation of the equations that applied to both laminar and turbulent flows. However the other person whose name is attached with Navier is the Irish mathematician-physicist George Gabriel Stokes. In 1845 he published a derivation of the equations in a manner that is currently understood.

ChessisGood

Feb 16, 2012

#76

References

C. L. M. H. Navier (1822), "Memoire sur les lois du mouvement des fluides", Mem. Acad. Sci. Inst. France, 6, 389-440.

ChessisGood

Feb 16, 2012

#77

External links

Huskie99

Feb 16, 2012

#78

The answer is 6 - come on, give us a tough one!!

sapientdust

Feb 16, 2012

#79

Bravo! You have truly mastered all the subtleties of the art of cut-and-paste.

nameno1had

Feb 16, 2012

#80

There are a lot of variables to calculate in fluid dynamics. The compounded nature of some of these variables as they effect one another in various environments would make it difficult to formulate certain scenarios. It isn't as simple as setting up dominos and predicting where they will fall. Does anyone have a super computer they aren't using and I'll get right on these problems... I should be done in no time.