伯努利方程及其推导过程一

伯努利方程及其推导过程一

本文主要来源于http://course.bitsde.com/web/209/ch02/16.htm，并做了部分修改和优化；

一、理想液体的运动微分方程

在微小流束上，取截面积为dA，长为ds的微元体，现研究理想液体定常流动条件下在重力场中沿流线运动时其力的平衡关系。

微元体的所受的重力为$-\rho gdAds$，压力作用在两端面上的力为$pdA-(p+\frac{\partial p}{\partial s}ds)dA$；
微元体在定常流动下的加速度为：
$$a=\frac{u}{t}=\frac{\partial u}{\partial s}\cdot \frac{ds}{dt}=u\frac{\partial u}{\partial s}$$
微元体的力平衡方程为：
$$pdA-(p+\frac{\partial p}{\partial s}ds)dA-\rho gdAds=\rho dAds\cdot u\frac{\partial u}{\partial s}$$
上式简化后可得：
$$\frac{1}{\rho }\cdot \frac{\partial p}{\partial s}+g\frac{\partial z}{\partial s}+u\frac{\partial u}{\partial s}=0$$
p、z、u都只是s的函数，进一步简化为：
$$\frac{1}{\rho }\cdot dp+gdz+udu=0$$
上式即为重力场中，理想液体沿流线作定常流动时的运动方程，即欧拉运动方程。

二、理想液体的伯努利方程

沿流线对欧拉运动方程积分得：
$$\frac{p}{\rho }+gz+\frac{u^2}{2}=C$$
上式两边同除以g得：
$$\frac{p_1}{\rho g}+z_1+\frac{u_1^2}{2g}=\frac{p_2}{\rho g}+z_2+\frac{u_2^2}{2g}$$
以上两式即为理想液体作定常流动的伯努利方程。
物理意义：
第一项为单位重量液体的压力能称为比压能（$\frac{p}{\rho g}$）；
第二项为单位重量液体的位能称为比位能（z）。
第三项为单位重量液体的动能称为比动能（ $\frac{u^2}{2g}$ ）；
由于上述三种能量都具有长度单位，故又分别称为压力水头、速度水头和位置水头。三者之间可以互相转换，但总和（H，称为总水头）为一定值。

三、实际液体流束的伯努利方程

实际液体都具有粘性，因此液体在流动时还需克服由于粘性所引起的摩擦阻力，这必然要消耗能量，设因粘性而消耗的能量为，则实际液体微小流束的伯努利方程为：
$$\frac{p_1}{\rho g}+z_1+\frac{u_1^2}{2g}=\frac{p_2}{\rho g}+z_2+\frac{u_2^2}{2g}+h_w^{\prime }$$

四、实际液体总流的伯努利方程

将微小流束扩大到总流，由于在通流截面上速度u是一个变量，若用平均流速代替，则必然引起动能偏差，故必须引入动能修正系数。于是实际液体总流的伯努利方程为：
$$\frac{p_1}{\rho g}+z_1+\frac{a_1{u}_1^2}{2g}=\frac{p_2}{\rho g}+z_2+\frac{a_2{u}_2^2}{2g}+h_w^{\prime }$$
式中：$h_w^{\prime }$—由液体粘性引起的能量损失；
$a_2$ 、$a_2$—动能修正系数，一般在紊流时取＝1，层流时取＝2。