Lieb's square ice constant: Difference between revisions
en>David Eppstein remove unexplained and stale expert tag |
en>David Eppstein rewrite for clarity |
||
Line 1: | Line 1: | ||
The '''Cauchy momentum equation''' is a vector [[partial differential equation]] put forth by [[Cauchy]] that describes the non-relativistic momentum transport in any [[Continuum mechanics|continuum]]:<ref>{{cite book | |||
| last = Acheson | |||
| first = D. J. | |||
| title = Elementary Fluid Dynamics | |||
| publisher = [[Oxford University Press]] | |||
| year = 1990 | |||
|page = 205 | |||
| isbn = 0-19-859679-0}}</ref> | |||
:<math>\rho \frac{D \mathbf{v}}{D t} = \nabla \cdot \boldsymbol{\sigma} + \mathbf{f}</math> | |||
or, with the [[convective derivative|material derivative]] expanded out, | |||
:<math>\rho \left[\frac{\partial \mathbf{v}}{\partial t} + (\mathbf{v} \cdot \nabla) \mathbf{v}\right] = \nabla \cdot \boldsymbol{\sigma} + \mathbf{f}</math> | |||
where <math>\rho</math> is the [[density]] of the continuum, <math>\boldsymbol{\sigma}</math> is the [[Cauchy stress tensor|stress tensor]], and <math>\mathbf{f}</math> contains all of the [[body force]]s per unit volume (often simply [[density]] times [[gravity]]). <math>\mathbf{v}</math> is the velocity vector field, which depends on time and space. | |||
The stress tensor is sometimes split into pressure and the deviatoric stress tensor: | |||
:<math>\boldsymbol{\sigma} = -p\mathbb{I} + \mathbb{T}</math> | |||
where <math>\scriptstyle \mathbb{I}</math> is the <math>\scriptstyle 3 \times 3</math> [[identity matrix]] and <math>\scriptstyle \mathbb{T}</math> the deviatoric stress tensor. The divergence of the stress tensor can be written as | |||
:<math>\nabla \cdot \boldsymbol{\sigma} = -\nabla p + \nabla \cdot\mathbb{T}.</math> | |||
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]] (<math>i^{th}</math> component) to a [[control volume]] in the continuum being modeled gives: | |||
:<math>m a_i = F_i\,</math> | |||
:<math>\rho \int_{\Omega} \frac{d u_i}{d t} \, dV = \int_{\Omega} \nabla_j\sigma_i^j \, dV + \int_{\Omega} f_i \, dV</math> | |||
:<math> \int_{\Omega} (\rho \frac{d u_i}{d t} - \nabla_j\sigma_i^j - f_i )\, dV = 0</math> | |||
:<math> \rho \dot{u_i} - \nabla_j\sigma_i^j - f_i = 0 </math> | |||
where <math>\Omega</math> 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 [[tensor derivative|derivative]] of the stress tensor is one of the forces that constitutes <math>F_i</math>. | |||
===Cartesian coordinates=== | |||
:<math>\begin{align} | |||
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{align} | |||
</math> | |||
===Cylindrical coordinates=== | |||
:<math> | |||
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_{z r}}}{\partial z} - \frac {\tau_{\phi \phi}}{r} + \rho g_r</math> | |||
:<math> | |||
\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_{z r}}}{\partial z} + \rho g_{\phi}</math> | |||
:<math> | |||
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_{z z}}}{\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.</math> | |||
By expressing the shear stress in terms of [[viscosity]] and fluid [[Shear velocity|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 == | |||
*[[Navier–Stokes equations]] | |||
*[[Computational fluid dynamics]] | |||
==References== | |||
{{Reflist}} | |||
[[Category:Concepts in physics]] | |||
[[Category:Continuum mechanics]] | |||
[[Category:Partial differential equations]] |
Revision as of 01:58, 15 September 2012
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]
or, with the material derivative expanded out,
where is the density of the continuum, is the stress tensor, and contains all of the body forces per unit volume (often simply density times gravity). is the velocity vector field, which depends on time and space.
The stress tensor is sometimes split into pressure and the deviatoric stress tensor:
where is the identity matrix and the deviatoric stress tensor. The divergence of the stress tensor can be written as
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 ( component) to a control volume in the continuum being modeled gives:
where 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 .
Cartesian coordinates
Cylindrical coordinates
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
43 year old Petroleum Engineer Harry from Deep River, usually spends time with hobbies and interests like renting movies, property developers in singapore new condominium and vehicle racing. Constantly enjoys going to destinations like Camino Real de Tierra Adentro.
- ↑ 20 year-old Real Estate Agent Rusty from Saint-Paul, has hobbies and interests which includes monopoly, property developers in singapore and poker. Will soon undertake a contiki trip that may include going to the Lower Valley of the Omo.
My blog: http://www.primaboinca.com/view_profile.php?userid=5889534