牛顿流体动量控制方程

流体动量控制方程

Equation of Motion for a Newtonian Fluid with Constant $\rho$ and $\mu$

控制方程通式：

$$\rho \dfrac {D \pmb v}{D t}=- \nabla p + \mu \nabla^2 \pmb v+\rho \pmb g$$

---

1.直角坐标系($x, y, z$)

直角坐标系Cartesian coordinates ($\textit { x,y,z }$):	NO.
$\rho \left (\dfrac{\partial v_x}{\partial t}+v_x \dfrac{\partial v_x}{\partial x}+v_y \dfrac{\partial v_x}{\partial y} + v_z \dfrac{\partial v_x}{\partial z}\right)=- \dfrac{\partial p}{\partial x} + \mu\left [\dfrac {\partial^2 v_x}{\partial x}+\dfrac {\partial^2 v_x}{\partial y}+\dfrac {\partial^2 v_x}{\partial z}\right]+\rho g_x$	1-1
$\rho \left (\dfrac{\partial v_y}{\partial t}+v_x \dfrac{\partial v_y}{\partial x}+v_y \dfrac{\partial v_y}{\partial y} + v_z \dfrac{\partial v_y}{\partial z}\right)=- \dfrac{\partial p}{\partial y} + \mu\left [\dfrac {\partial^2 v_y}{\partial x}+\dfrac {\partial^2 v_y}{\partial y}+\dfrac {\partial^2 v_y}{\partial z}\right]+\rho g_y$	1-2
$\rho \left (\dfrac{\partial v_z}{\partial t}+v_x \dfrac{\partial v_z}{\partial x}+v_y \dfrac{\partial v_z}{\partial y} + v_z \dfrac{\partial v_z}{\partial z}\right)=- \dfrac{\partial p}{\partial z} + \mu\left [\dfrac {\partial^2 v_z}{\partial x}+\dfrac {\partial^2 v_z}{\partial y}+\dfrac {\partial^2 v_z}{\partial z}\right]+\rho g_z$	1-3

---

2.圆柱坐标系（$r,\theta, z$)

圆柱坐标系Cylindrical coordinates coordinates ($\textit {r, $\theta$, z }$):	NO.
$\rho \left (\dfrac{\partial v_r}{\partial t}+v_r \dfrac{\partial v_r}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_r}{\partial \theta} + v_z \dfrac{\partial v_r}{\partial z}-\dfrac {v_\theta^2}{r}\right)=- \dfrac{\partial p}{\partial r} + \mu \left [\dfrac {\partial}{\partial r} (\dfrac{1}{r} \dfrac {\partial}{\partial r} (rv_r))+\dfrac{1}{r^2}\dfrac {\partial^2 v_r}{\partial \theta^2}+\dfrac {\partial^2 v_r}{\partial z^2}-\dfrac {2}{r^2}\dfrac{\partial v_\theta}{\partial \theta}\right]+\rho g_r$	2-1
$\rho \left (\dfrac{\partial v_\theta}{\partial t}+v_r \dfrac{\partial v_\theta}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_\theta}{\partial \theta} + v_z \dfrac{\partial v_\theta}{\partial z}+\dfrac {v_r v_\theta}{r}\right)=- \dfrac{1}{r}\dfrac{\partial p}{\partial \theta} + \mu \left [\dfrac {\partial}{\partial r} (\dfrac{1}{r} \dfrac {\partial}{\partial r} (rv_\theta))+\dfrac {1}{r^2} \dfrac {\partial^2 v_\theta}{\partial \theta^2}+\dfrac {\partial^2 v_\theta}{\partial z^2} + \dfrac {2}{r^2} \dfrac {\partial v_r}{\partial \theta}\right]+\rho g_\theta$	2-2
$\rho \left (\dfrac{\partial v_z}{\partial t}+v_r \dfrac{\partial v_z}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_z}{\partial \theta} + v_z \dfrac{\partial v_z}{\partial z}\right)=- \dfrac{\partial p}{\partial z} + \mu \left [\dfrac{1}{r}\dfrac {\partial}{\partial r} (r \dfrac {\partial v_z}{\partial r})+\dfrac{1}{r^2}\dfrac {\partial^2 v_z}{\partial \theta^2}+\dfrac {\partial^2 v_z}{\partial z^2}\right]+\rho g_z$	2-3

---

3.球坐标系（$r, \theta, \phi$)

球坐标系Spherical coordinates（$\textit {r, $\theta$, $\phi$ }$):	NO.
$\rho \left (\dfrac{\partial v_r}{\partial t}+v_r \dfrac{\partial v_r}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_r}{\partial \theta} + \dfrac{v_\phi}{r sin\theta} \dfrac{\partial v_r}{\partial \phi}-\dfrac {v_\theta^2+v_\phi^2}{r}\right)=- \dfrac{\partial p}{\partial r} + \mu \left [\dfrac{1}{r^2}\dfrac {\partial^2}{\partial r^2} (r^2 v_r)+\dfrac{1}{r sin\theta}\dfrac {\partial}{\partial \theta}(sin\theta \dfrac {\partial v_r}{\partial \theta})+\dfrac{1}{r^2 sin^2\theta}\dfrac {\partial^2 v_r}{\partial \phi^2}\right]+\rho g_r$	3-1
$\rho \left (\dfrac{\partial v_\theta}{\partial t}+v_r \dfrac{\partial v_\theta}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_\theta}{\partial \theta} + \dfrac{v_\phi}{r sin\theta} \dfrac{\partial v_\theta}{\partial \phi} + \dfrac {v_r v_\theta - v_\phi^2 cot \theta}{r}\right)=- \dfrac {1}{r}\dfrac{\partial p}{\partial \theta} + \mu \left [\dfrac{1}{r^2}\dfrac {\partial}{\partial r} (r^2 \dfrac {\partial v_\theta}{\partial r})+\dfrac{1}{r^2}\dfrac {\partial}{\partial \theta} (\dfrac {1}{sin \theta} \dfrac {\partial}{\partial \theta}(v_\theta sin \theta))+ \dfrac{1}{r^2 sin^2\theta}\dfrac {\partial^2 v_\theta}{\partial \phi^2}+\dfrac {2}{r^2} \dfrac{\partial v_r}{\partial \theta} - \dfrac {2cot\theta}{r^2 sin \theta} \dfrac{\partial v_\phi}{\partial \phi}\right]+\rho g_\theta$	3-2
$\rho \left (\dfrac{\partial v_\phi}{\partial t}+v_r \dfrac{\partial v_\phi}{\partial r}+\dfrac {v_\theta}{r} \dfrac{\partial v_\phi}{\partial \theta} + \dfrac{v_\phi}{r sin\theta} \dfrac{\partial v_\phi}{\partial \phi} + \dfrac {v_\phi v_r + v_\theta v_\phi cot \theta}{r}\right)=- \dfrac {1}{r sin\theta}\dfrac{\partial p}{\partial \phi} + \mu \left [\dfrac{1}{r^2}\dfrac {\partial}{\partial r} (r^2 \dfrac {\partial v_\phi}{\partial r})+\dfrac{1}{r^2}\dfrac {\partial}{\partial \theta} (\dfrac {1}{sin \theta} \dfrac {\partial}{\partial \theta}(v_\phi sin \theta))+ \dfrac{1}{r^2 sin^2\theta}\dfrac {\partial^2 v_\phi}{\partial \phi^2} + \dfrac {2}{r^2 sin \theta} \dfrac{\partial v_r}{\partial \phi} + \dfrac {2cot\theta}{r^2 sin \theta} \dfrac{\partial v_\theta}{\partial \phi}\right]+\rho g_\phi$	3-3

---

参考文献

1. R. Byron Bird, Warren E. stewart, Edwin N. Lightfoot.* Transport phenomena:Revised second edition* John Wiely &Sons, Inc.