Cauchy momentum equation

From Wikipedia, the free encyclopedia

The Cauchy momentum equation is a vector partial differential equation put forth by Cauchy that describes the non-relativistic momentum transport in any continuum:[1]

\rho {\frac {D{\mathbf {v}}}{Dt}}=\nabla \cdot {\boldsymbol {\sigma }}+{\mathbf {f}}

or, with the material derivative expanded out,

\rho \left[{\frac {\partial {\mathbf {v}}}{\partial t}}+({\mathbf {v}}\cdot \nabla ){\mathbf {v}}\right]=\nabla \cdot {\boldsymbol {\sigma }}+{\mathbf {f}}

where \rho is the density of the continuum, {\boldsymbol {\sigma }} is the stress tensor, and {\mathbf {f}} contains all of the body forces per unit volume (often simply density times gravity). {\mathbf {v}} is the velocity vector field, which depends on time and space.

The stress tensor is sometimes split into pressure and the deviatoric stress tensor:

{\boldsymbol {\sigma }}=-p{\mathbb {I}}+{\mathbb {T}}

where \scriptstyle {\mathbb {I}} is the \scriptstyle 3\times 3 identity matrix and \scriptstyle {\mathbb {T}} the deviatoric stress tensor. The divergence of the stress tensor can be written as

\nabla \cdot {\boldsymbol {\sigma }}=-\nabla p+\nabla \cdot {\mathbb {T}}.

All non-relativistic momentum conservation equations, such as the Navier–Stokes equation, can be derived by beginning with the Cauchy momentum equation and specifying the stress tensor through a constitutive relation.

Derivation

Applying Newton's second law (i^{{th}} component) to a control volume in the continuum being modeled gives:

ma_{i}=F_{i}\,
\rho \int _{{\Omega }}{\frac {du_{i}}{dt}}\,dV=\int _{{\Omega }}\nabla _{j}\sigma _{i}^{j}\,dV+\int _{{\Omega }}f_{i}\,dV
\int _{{\Omega }}(\rho {\frac {du_{i}}{dt}}-\nabla _{j}\sigma _{i}^{j}-f_{i})\,dV=0
\rho {\dot {u_{i}}}-\nabla _{j}\sigma _{i}^{j}-f_{i}=0

where \Omega represents the control volume. Since this equation must hold for any control volume, it must be true that the integrand is zero, from this the Cauchy momentum equation follows. The main challenge in deriving this equation is establishing that the derivative of the stress tensor is one of the forces that constitutes F_{i}.

Cartesian coordinates

{\begin{aligned}x:\;\;\rho \left({\frac {\partial u_{x}}{\partial t}}+u_{x}{\frac {\partial u_{x}}{\partial x}}+u_{y}{\frac {\partial u_{x}}{\partial y}}+u_{z}{\frac {\partial u_{x}}{\partial z}}\right)&=-{\frac {\partial P}{\partial x}}+{\frac {\partial \tau _{{xx}}}{\partial x}}+{\frac {\partial \tau _{{xy}}}{\partial y}}+{\frac {\partial \tau _{{xz}}}{\partial z}}+\rho g_{x}\\y:\;\;\rho \left({\frac {\partial u_{y}}{\partial t}}+u_{x}{\frac {\partial u_{y}}{\partial x}}+u_{y}{\frac {\partial u_{y}}{\partial y}}+u_{z}{\frac {\partial u_{y}}{\partial z}}\right)&=-{\frac {\partial P}{\partial y}}+{\frac {\partial \tau _{{yx}}}{\partial x}}+{\frac {\partial \tau _{{yy}}}{\partial y}}+{\frac {\partial \tau _{{yz}}}{\partial z}}+\rho g_{y}\\z:\;\;\rho \left({\frac {\partial u_{z}}{\partial t}}+u_{x}{\frac {\partial u_{z}}{\partial x}}+u_{y}{\frac {\partial u_{z}}{\partial y}}+u_{z}{\frac {\partial u_{z}}{\partial z}}\right)&=-{\frac {\partial P}{\partial z}}+{\frac {\partial \tau _{{zx}}}{\partial x}}+{\frac {\partial \tau _{{zy}}}{\partial y}}+{\frac {\partial \tau _{{zz}}}{\partial z}}+\rho g_{z}.\end{aligned}}

Cylindrical coordinates

r:\;\;\rho \left({\frac {\partial u_{r}}{\partial t}}+u_{r}{\frac {\partial u_{r}}{\partial r}}+{\frac {u_{{\phi }}}{r}}{\frac {\partial u_{r}}{\partial \phi }}+u_{z}{\frac {\partial u_{r}}{\partial z}}-{\frac {u_{{\phi }}^{2}}{r}}\right)=-{\frac {\partial P}{\partial r}}+{\frac {1}{r}}{\frac {\partial {(r{\tau _{{rr}})}}}{\partial r}}+{\frac {1}{r}}{\frac {\partial {\tau _{{\phi r}}}}{\partial \phi }}+{\frac {\partial {\tau _{{zr}}}}{\partial z}}-{\frac {\tau _{{\phi \phi }}}{r}}+\rho g_{r}
\phi :\;\;\rho \left({\frac {\partial u_{{\phi }}}{\partial t}}+u_{r}{\frac {\partial u_{{\phi }}}{\partial r}}+{\frac {u_{{\phi }}}{r}}{\frac {\partial u_{{\phi }}}{\partial \phi }}+u_{z}{\frac {\partial u_{{\phi }}}{\partial z}}+{\frac {u_{r}u_{{\phi }}}{r}}\right)=-{\frac {1}{r}}{\frac {\partial P}{\partial \phi }}+{\frac {1}{r}}{\frac {\partial {\tau _{{\phi \phi }}}}{\partial \phi }}+{\frac {1}{r^{2}}}{\frac {\partial {(r^{2}{\tau _{{r\phi }})}}}{\partial r}}+{\frac {\partial {\tau _{{zr}}}}{\partial z}}+\rho g_{{\phi }}
z:\;\;\rho \left({\frac {\partial u_{z}}{\partial t}}+u_{r}{\frac {\partial u_{z}}{\partial r}}+{\frac {u_{{\phi }}}{r}}{\frac {\partial u_{z}}{\partial \phi }}+u_{z}{\frac {\partial u_{z}}{\partial z}}\right)=-{\frac {\partial P}{\partial z}}+{\frac {\partial {\tau _{{zz}}}}{\partial z}}+{\frac {1}{r}}{\frac {\partial {\tau _{{\phi z}}}}{\partial \phi }}+{\frac {1}{r}}{\frac {\partial {(r{\tau _{{rz}})}}}{\partial r}}+\rho g_{z}.

By expressing the shear stress in terms of viscosity and fluid velocity, and assuming constant density and viscosity, the Cauchy momentum equation will lead to the Navier–Stokes equations. By assuming inviscid flow, the Navier–Stokes equations can further simpify to the Euler equations.


See also

References

  1. ↑ Acheson, D. J. (1990). Elementary Fluid Dynamics. Oxford University Press. p. 205. ISBN 0-19-859679-0. 
This article is issued from Wikipedia. The text is available under the Creative Commons Attribution/Share Alike; additional terms may apply for the media files.