流体力学精品课件

1、流体力学第1页，共38页，2022年，5月20日，10点54分，星期四 从微观的角度而言不论液体还是气体其分子与分子之间都是存在间隙的，例如海平面条件下，空气分子的平均自由程为 l 10-8 mm，但是这个距离与我们宏观上关心的物体（如飞行器）的任何一个尺寸 L 相比较都是微乎其微的， l / L 1 流体力学和空气动力学是从宏观上研究流体（空气）的运动规律和作用力（流体内部和流体对物体）的规律的学科，流体力学和空气动力学常用“介质”一词表示它所处理 的流体，流体包含液体和气体1.1 流体属性1.1.1 连续介质的概念第2页，共38页，2022年，5月20日，10点54分，星期四 当受到物体扰

2、动时，流体或空气所表现出的是大量分子运动体现出的宏观特性变化如压强、密度等，而不是个别分子的行为。 一旦满足连续介质假设，就可以把流体的一切物理性质如密度、压强、温度及宏观运动速度等表为空间和时间的连续可微函数，便于用数学分析工具来解决问题。 一般用努生数即分子平均自由程与物体特征尺寸之比来判断流体是否满足连续介质假设 l / L 0体积减小dv0，考虑到一定质量的流体 m=v = 常数, 其密度与体积成反比: 体积弹性模量可写为： (N/m2) 当E较大，则流体不容易被压缩，反之当E较小则流体容易被压缩。液体的E一般较大，通常可视为不可压缩流体，气体的E通常较小，且与热力过程有关，故 气体具

3、有压缩性。对具体流动问题是否应考虑空气压缩性要看流动产生的压强变化是否引起密度显著变化，一般情况下，当空气流动速度较低时，压强变化引起的密度变化很小，可不考虑空气压缩性对流动特性的影响。1.1.3 流体的压缩性与弹性、空气的流动性第7页，共38页，2022年，5月20日，10点54分，星期四 气体是流体的一种，它具有流动性。 气体受到扰动后，扰动的影响将会以波动的形式传播开去，扰动传播的速度即为声速，因此扰动的传播与气体的弹性有关。后面讲到高速流动时会证明，这里的等于声速的平方。所以气体的弹性决定于它的密度和声速。 而当飞行速度超过音速之后，扰动传播的速度仍是声速，相对于飞行速度而言，它就慢了

4、，飞机没有飞到跟前，空气微团是没有预感的，只是飞到跟前时才突然地被推开，这时流动性就很差了。 对于飞行器而言，单说空气的流动性就不够了，而必须在飞行器的飞行速度和扰动的传播速度的比值之下来谈流动性。 当飞行速度远小于音速时（低速飞行）,扰动在空气里传播速度相对于飞行速度而言是很快的，这时流动性很好。 飞行速度再大上去，到了高超音速范围，空气简直像没有流动性一样，而像固体的粒子那样向飞行器打来 。1.1.3 流体的压缩性与弹性、空气的流动性第8页，共38页，2022年，5月20日，10点54分，星期四飞行器的飞行速度 v 和扰动的传播速度 a 的比值称为马赫数 Ma：由于气体的弹性决定于声速，因

5、此马赫数的大小可以看成是气体相对压缩性的一个指标 当飞行速度远小于音速时（低速飞行），即马赫数较小时 ，可以认为此时流动的弹性影响相对较大，即压缩性影响相对较小，从而低速气体有可能被当作不可压缩流动来处理。 而当马赫数较大之后，可以认为此时流动的弹性影响相对较小，即压缩性影响相对较大，从而气体就不能被当作不可压缩流动来处理，而必须考虑流动的压缩性效应。 可以证明，近似划分气体压缩性影响的马赫数界线为 Ma 0.3 , 即当马赫数小于0.3时，气体的压缩性影响可以忽略不计，或者换言之，此时流动速度的变化不会引起气体密度的显著变化。1.1.3 流体的压缩性与弹性、空气的流动性第9页，共38页，20

6、22年，5月20日，10点54分，星期四 实际流体都有粘性，不过有大有小，空气和水的粘性都不算大，日常生活中人们不会理会它，但观察河流岸边的漂浮物可以看到粘性的存在。下述直匀流流过平板表面的实验突出表明了粘性的影响：yu 由于粘性影响，原来是均匀的气流流至平板后直接贴着板面的一层速度降为零，称为流体与板面间无滑移。稍外一层的气流受到层间摩擦作用速度也也下降至接近于零，但由于不紧挨板面多少有些速度，层间的互相牵扯作用一层层向外传递，离板面一定距离后，牵扯作用逐步消失，速度分布变为均匀。 1.1.4 流体的粘性第10页，共38页，2022年，5月20日，10点54分，星期四 取其中相邻的二层流体来

7、看，慢层对快层有向后的牵扯而使其有变慢的趋势，而快层对慢层有向前的牵扯使其有变快的趋势uu+du 流体相邻层间存在着抵抗层间相互错动的趋势这一特性称为流体的粘性，层间的这一抵抗力即摩擦力或剪切力，单位面积上的剪切力称为剪切应力 牛顿提出，流体内部的剪切力与流体的速度梯度成正比 dytt+dt其中比例系数是反映粘性大小的物性参数，称为动力粘性系数 1.1.4 流体的粘性第11页，共38页，2022年，5月20日，10点54分，星期四此即著名的牛顿粘性公式。其中的单位是帕：N/m2，动力粘性系数的单位是：帕秒：Ns/m2从牛顿粘性公式可以看出：1. 流体的剪应力与压强 p 无关（注意到固体摩擦力与

8、正压力有关)。2. 当 0 时， 即无论剪应力多小，只要存在剪应力，流体就会发生变形运动，因此牛顿粘性公式可看成是易流性的数学表达。 3. 当 时， 0，即只要流体静止或无变形，就不存在剪应力，换言之，流体不存在静摩擦力。 1.1.4 流体的粘性第12页，共38页，2022年，5月20日，10点54分，星期四 液体和气体产生粘性的物理原因不同，液体分子结构紧密，液体的粘性主要来自于液体分子间的内聚力，气体分子结构松散，气体粘性主要来自于气体分子的热运动，因此液体和气体的动力粘性系数随温度的变化趋势刚好相反，但粘性系数与压强基本无关 液体和气体的动力粘性系数随温度变化的关系可查阅相应表格或近似公

9、式，如气体动力粘性系数的萨特兰公式，等等。 在许多空气动力学问题里，粘性力和惯性力同时存在，在式子中和往往以（/ ）的组合形式出现，用符号表示因为量纲只包含长度和时间，为运动学量，称为运动粘性系数。液体与气体动力粘性系数随温度变化的趋势为： 液体： 温度升高，动力粘性系数变小，反之变大 气体： 温度升高，动力粘性系数变大，反之变小 空气粘性不大,初步近似可忽略其粘性作用，忽略粘性的流体称为理想流体。 1.1.4 流体的粘性第13页，共38页，2022年，5月20日，10点54分，星期四按照作用力的性质和作用方式，可分为彻体力和表面力两类彻体力：外力场作用于流体微团质量中心，大小与微团质量成正比

10、的非接触力，例如重力，惯性力和磁流体具有的电磁力等都属于彻体力，也有称为体积力或质量力的，由于彻体力按质量分布，故一般用单位质量的彻体力表示，并且往往写为分量形式：其中 是微团体积，为密度， 为作用于微团的彻体力， i 、j、 k分别是三个坐标方向的单位向量，fx 、fy 、fz 分别是三个方向的单位质量彻体力分量 。1.2 作用在流体微团上的力的分类第14页，共38页，2022年，5月20日，10点54分，星期四表面力：相邻流体或物体作用于所研究流体团块外表面，大小与流体团块表面积成正比的接触力。由于表面力按面积分布，故用单位面积上的接触力即接触应力表示，由于接触应力一般与表面法线方向并不重

11、合，故又可以将接触应力分解为法向应力和切向应力ATPn法向应力即静压强，切向应力即摩擦应力或剪切应力：上述画出的表面力对整个流体而言是内力，对所画出的流体团块来说则是外力。1.2 作用在流体微团上的力的分类第15页，共38页，2022年，5月20日，10点54分，星期四 流体内任取一个剖面一般有法向应力和切向应力，但切向应力完全是由粘性产生的，而流体的粘性力只有在流动时才存在，静止流体是不能承受切向应力的 流体中的法向应力称为压强 p（注），其指向沿着表面的内法线方向，压强的量纲是力/长度2，单位为（N/m2）或 （帕：Pa） 在理想（无粘）流体中，不论流体静止还是运动，尽管一般压强是位置的函

12、数 p=p(x,y,z), 但在同一点处压强不因受压面方位不同而变化，这个结果称为理想流体内压强是各向同性的。 注：关于有粘的运动流体，严格说来压强指的是三个互相垂直方向的法向力的平均值，加负号 。1.3 理想流体内一点的压强及其各向同性第16页，共38页，2022年，5月20日，10点54分，星期四 如讨论P点处压强，在周围取如图微元4面体ABCO,作用在各表面的压强如图所示，理想流体无剪切应力，由于dx、dy、dz的取法任意，故面ABC的法线方向n方向也是任意的yxzdxdydzpzpxpypnnABCoP分别沿 x、y、z三个方向建立力的平衡关系：x方向合外力质量加速度（x方向）方程左端

13、等于：方程右端等于：三阶小量0，由此可得： 因为图中的n方向为任取，故各向同性得证同理可得：即：1.3 理想流体内一点的压强及其各向同性第17页，共38页，2022年，5月20日，10点54分，星期四下面我们来研究压强在静止流体中的分布规律。在平衡流体（静止或相对静止）中取定一笛卡儿坐标系oxyz，坐标轴方位任意。在流体内取定一点P（x ,y ,z）,然后以该点为中心点沿坐标轴三个方向取三个长度dx,dy,dz,划出一微元六面体作为分析对象:假设：六面体体积：d=dxdydz中心点坐标： x ,y ,z中心点压强：p = p（x,y ,z）中心点密度： =（x,y,z）中心点处沿三个方向的单位

14、质量彻体力 fx, fy, fz微元六面体的表面力可以用中心点处压强的一阶泰勒展开表示,如图为x方向表面力，其他方向同理可得。由于流体静止故无剪应力。xyzPdxdydz1.4 流体静平衡微分方程第18页，共38页，2022年，5月20日，10点54分，星期四x方向的表面力为：x方向的彻体力为：流体静止，则x方向的合外力为零：两边同除以 d=dxdydz 并令 d 趋于零，可得x方向平衡方程：同理可得y,z方向的平衡方程：流体平衡微分方程1.4 流体静平衡微分方程第19页，共38页，2022年，5月20日，10点54分，星期四这三个式子表明当流体平衡时，若压强在某个方向有梯度的话，必然是由于彻

15、体力在该方向有分量的缘故将上三个式子分别乘以dx，dy，dz，然后相加起来，得到：此式左端是个全微分：如果右端的三个彻体力分量 fx ，fy ，fz ，符合下列关系的话：则第一式右端括号也是某个函数的全微分。记该函数为（x,y,z），称为彻体力的势函数，或称彻体力为有势力。重力、惯性力和电磁力都为有势力。1.4 流体静平衡微分方程第20页，共38页，2022年，5月20日，10点54分，星期四当彻体力为有势力时，则平衡微分方程可写为：设彻体力与势函数的关系为：等压面的概念：流场中压强相等的空间点组成的几何曲面或平面p=c等压面在等压面上满足：上式积分后为一几何曲面或平面，该曲面上满足 dp=0

16、，上方程称为等压面方程如果我们知道某一点的压强值 pa 和彻体力势函数 a 的值,则任何其它点的压强和势函数之间的关系便可表为：或：1.4 流体静平衡微分方程第21页，共38页，2022年，5月20日，10点54分，星期四等压面方程还可写为：其中： 为彻体力向量为等压面上的向径等压面上式表明：等压面处处与彻体力相正交例如：1.在重力场下静止液体等压面必然为水平面gaa3. 在水平向右加速容器中的液体除了受到重力之外，还受到向左的惯性力，二者合成的彻体力向左下方，因此等压面是向右倾斜的平面2. 在加速上升电梯中的液体除了受到重力之外，还受到向下的惯性力，二者合成的彻体力均为向下，因此等压面也是水

17、平面1.4 流体静平衡微分方程第22页，共38页，2022年，5月20日，10点54分，星期四1.5 重力场静止液体中的压强分布规律设封闭容器自由面处压强为p0，如图建立坐标系，考虑距水平轴高度为 y 处的某单位质量流体，其彻体力可表示为：积分得（注意g）：此式称为平衡基本方程其中g为重力加速度，代入平衡微分方程得：或：上式表明，在平衡流体中p/与y之和为常数。显然，静止流体中等压面为水平面 ycp0。xygy第23页，共38页，2022年，5月20日，10点54分，星期四对于不同高度上的1、2两点，平衡基本方程可以写为:平衡流体中不同高度处，压力水头与高度水头可以互相转换，但总水头保持不变。

18、 的几何意义为:y -代表所研究流体质点在坐标系中所处高度，称为高度水头 p/-代表所研究流体质点在真空管中上升高度，称为压力水头H-由于方程量纲为高度，该积分常数代表上述二高度之和称为总水头，如图所示：y2。11yxp0。yH真空1.5 重力场中静止液体中的压强分布规律第24页，共38页，2022年，5月20日，10点54分，星期四 的物理意义为:y -代表单位重量流体的重力势能简称势能 p/-代表单位重量流体的压力势能简称压力能H -代表平衡流体中单位重量流体的总能量平衡基本方程 表明:平衡流体中势能与压力能可以互相转换，但总能量保持不变假设自由液面距水平轴距离为H，则自由面与 y 处流体

19、满足：。xygp0yHh其中 h = H-y 是所论液体距自由面的深度1.5 重力场中静止液体中的压强分布规律第25页，共38页，2022年，5月20日，10点54分，星期四式 表明：平衡流体中距自由面深 h 处的压强来自于两部分的贡献： 一是上方单位面积上的液重h，因此压强随距自由面的淹没深度而线性增加 二是自由面上的压强贡献 P0，而该贡献处处相同与深度无关当自由面为大气压 pa 时，距自由面深h处的压强可表为：压强的计量：以真空为压强参考值计量的压强称为绝对压强，如上式中的 p以大气压pa为参考压强，高出大气压部分的压强称为相对压强 pb= p-pa以大气压pa为参考压强，不足大气压部分

20、的压强称为真空度 pv= pa-p对于同一个压强值 p，其相对压强 pb 与其真空度 pv 之间的关系为pb= -pv 1.5 重力场中静止液体中的压强分布规律第26页，共38页，2022年，5月20日，10点54分，星期四湿式大气压力计 例：湿式大气压力表的工作原理 有一种大气压力表是用汞柱的高度来表达大气压的数值的。一根上端封闭的长玻璃管和一个盛汞的底盒，玻管竖立。玻管中有汞与底盒中的汞连通。玻管中汞柱的 上端是真空的 。参看右图，把坐标平面xy放在管中汞柱的上表面，该处的按式 ，玻管下面与盒中汞面等高的A处(距上表面的深度为h)的压强 pA 是 而 pA 和大气压 pa 是相等的，即这样

21、，要计算大气压的值的话，只要把气压表上读下来的汞柱高度米乘以汞的重度就是了，大气压的读数往往只说汞柱高就行了，一个标准气压是760毫米汞柱。 1.5 重力场中静止液体中的压强分布规律第27页，共38页，2022年，5月20日，10点54分，星期四1.6 液体的相对平衡问题在以匀加速运动或匀角速度转动的相对平衡流体中，如果将坐标系固连在以匀加速运动或匀角速度转动的容器上，对液体引入惯性力(达朗伯原理)，则同样可以利用平衡微分方程求解问题。如图圆筒作匀角速转动 ， 求其中液体的等压面形状和压强分布规律。 yxr2r2y2x将坐标系固连于转筒，并建如图坐标系。考虑距底壁为 z ,半径为 r 处单位质

22、量流体，会受到一个向下的彻体力大小为 g ,此外还受到一个向外的惯性力大小为2r。在直角坐标系中，三个方向的彻体力可表为 yzg第28页，共38页，2022年，5月20日，10点54分，星期四求等压面：由等压面方程可得：积分得：即：为旋转抛物面族特别地，设自由面最低点距坐标原点高 H 时（如图），可定出自由面对应的常数：r = 0 时，c = z = H故自由面方程为其中 称为超高，即液面高出抛物线顶点的部分。yzgH1.6 液体的相对平衡问题第29页，共38页，2022年，5月20日，10点54分，星期四求压强分布：由平衡微分方程方程可得：积分得：由自由面条件，可定出积分常数：x = y =

23、 0 ， z = H 时, p = pa，定得积分常数 c = pa+g H， 带入上述积分结果，得如果令方括号等于H，则上式可以写为：其中 H 即为从自由面向下的淹没深度，等于超高加上距顶点的深度。yzgHH上述压强分布表明，在旋转平衡液体中，压强随深度线性增加，随半径呈平方增加1.6 液体的相对平衡问题第30页，共38页，2022年，5月20日，10点54分，星期四即在上图中，A点处压强大于自由面顶点处压强，而B点处压强又大于A点处压强， C点处压强又大于B点处压强。此外压强分布还与旋转角速度的平方 2 成正比，这是因为离心力是以旋转角速度的平方 2 成正比的，如旋转角速度很大，这个彻体力

24、可以很大 ，从而一定半径处的压强会很大。由于随半径不同各处的惯性离心力不同，因此合成的惯性力方向随半径而变化，在外侧惯性力较大故合彻体力方向趋于水平 ，在圆心附近惯性力较小故合彻体力方向趋于垂直，这是旋转平衡液体的等压面成为抛物面形状的原因。旋转液体的特点在在工程中也有很重要的应用，例如旋转铸造或离心铸造等，对于铸造薄壁容器、列车车轮等有重要意义右图为旋转液体压强分布演示：1.6 液体的相对平衡问题第31页，共38页，2022年，5月20日，10点54分，星期四1.7 标准大气 气象条件逐日都有些变化，更不用说不同的季节了,并且不同地区气象也不相同。无论做飞行器设计，还是做实验研究，都要用到大

25、气的条件，为了便于比较，工程上需要规定一个标准大气。这个标准是按中纬地区的平均气象条件定出来的。这样做计算时，都依此标准进行计算；做实验时，也都换算成标准条件下的数据。 标准大气规定在海平面上，大气温度为 15 或 T0 = 288.15K ，压强 p0 = 760 毫米汞柱 = 101325牛/米2，密度0 = 1.225千克/米3 从基准面到11 km的高空称为对流层，在对流层内大气密度和温度随高度有明显变化，温度随高度增加而下降，高度每增加1km，温度下降 6.5 K，即： 从 11 km 到 21km 的高空大气温度基本不变，称为同温层，在同温层内温度保持为 216.5 K。普通飞机主要在对流层和平流层里活动。 第32页，共38页，2022年，5月20日，10点54分，星期四 我们可以用静平衡微分方程把压强随高度而下降的变化规律推导出来。出发点只能是平衡微分方程，因为大气的密度是变量 。y（km）T（k）020406080100120160200240280320360400 高度大于 21km 以上时大气温度随高度的变化参见下图，大气温度随高度变化的原因复杂，主要因素有：地表吸收太阳热量、臭氧吸热与电离放热、空气或宇宙尘埃受短波辐射升温等。1.7 标准大气第33页，