不可压缩粘性流体的运动微分方程.pptx

第八章第八章 粘性流体绕过物体的流动粘性流体绕过物体的流动不可压缩粘性流体的运动微分方程及其解析解边界层理论圆柱体绕流物体的阻力第1页/共89页第2页/共89页第一节第一节 不可压缩粘性流体的运动微分方程不可压缩粘性流体的运动微分方程yxzodydzdx第一个下标表示应力所在平面的法线方向第二个下标表示应力本身的方向。第3页/共89页第4页/共89页根据达朗伯原理，作用于微元平行六面体上的各力对通过中心M并与z轴相平行的轴的力矩之和应等于零。又由于质量力和惯性力对该轴的力矩是四阶无穷小量，可以略去不计，故有Mdydxyxo图8-2 分析切向应力之间的关系用图第5页/共89页再略去四阶无穷小量，又因再略去四阶无穷小量，又因 ，故得，故得同理可得同理可得 （8-2）这样，在式（这样，在式（8-1）中九个应力只有六个是独立的。）中九个应力只有六个是独立的。假若流体的粘度在个方向上都是相同的可得广义牛顿内摩擦定律第6页/共89页其意义为：切向应力等于动力粘度和角变切向应力等于动力粘度和角变形速度的乘积。形速度的乘积。在粘性流体中，由于粘性的影响，流体微团除发生角变形以外，同时也发生线变形。第7页/共89页 现将切向应力和法向应力的关系式代入式(8-1),化简可得不可压缩粘性流体的运动微分方程：纳维-斯托克斯（Navier-Stokes）方程第8页/共89页 如果是没有粘性的理想流体，则 为零，于是纳维-斯托克斯方程变成理想流体的欧拉运动微分方程。如果没有加速度，则 都为零，于是上述方程变成欧拉平衡微分方程。所以说，上述纳维-斯托克斯方程式不可压缩流体的最普遍的运动微分方程。第9页/共89页 以上三式加上不可压缩流体的连续方程 或 共有四个方程，原则上可以求解不可压缩粘性流体运动问题中的四个未知数 和p。但是，实际上由于流体流动现象很复杂，要利用这四个方程去求解一般可压缩粘性流体的运动问题，在数学上还是很困难的。所以，求解纳维-斯托克斯方程，仍然是流体力学的一项重要任务。第10页/共89页 现在来研究粘性流体在大雷诺数下平滑地绕流某静止物体（例如机翼的翼型，图8-9）的情况。在紧靠物体表面的薄层内，流速将由物体表面上的零值迅速地增加到与来流速度 同数量级的大小，这种在大雷诺数下紧靠物大雷诺数下紧靠物体表面流速从零急剧增加到与来流速度相同数量级的薄层体表面流速从零急剧增加到与来流速度相同数量级的薄层称为称为边界层边界层。在边界层内，流体在物体表面法线方向上的速度梯度很大，即使粘度很小的流体，表现出的粘滞力也较大，决不能忽略。所以，边界层内的流体有相当大的涡通量。当边界层内的有旋流离开物体而流入下游时，在物体后形成尾涡区域。在边界层外，速度梯度很小，即使粘度较大的流体，粘滞力也很小，可以忽略不计。所以可以认为，在边界层外的流动是无旋的势流。第11页/共89页第12页/共89页由此可见，当粘性流体绕过物体流动时，可以将物体外面的流场划分为两个区域：两个区域：在边界层和尾涡区域内在边界层和尾涡区域内，必须考虑流体的粘滞力，它应当被看作是粘性流体的有旋流动粘性流体的有旋流动；在边界在边界层和尾涡区以外的区域内层和尾涡区以外的区域内，粘滞力很小，可以看作是理想流理想流体的无旋流动体的无旋流动。实际上，边界层内、外区域并没有一个明显的分界面，一般在实际应用中规定从固体壁面沿外法线到速一般在实际应用中规定从固体壁面沿外法线到速度达到势流速度的度达到势流速度的99%处距离为边界层的厚度处距离为边界层的厚度，以 表示，见图8-9。解决大雷诺数下绕过物体流动的近似方法是以边界层理论为基础的。用微型测速管直接测量紧靠机翼表面附近的流速得知，实际上边界层很薄，通常边界层的厚度仅为弦长的几百分之一。例如在汽轮机叶片出汽边上，最大边界层厚度一般为零点几毫米。从图8-9中可以看出，流体在前驻点O处速度为零，所以边界层的厚度在前驻点处等于零，然后沿着流动方向逐渐增加。为了清晰起见，在图8-9上将边界层的尺寸放大了。另外，边界层的外边界和流线并不重合，流线伸入第13页/共89页边界层内，这是由于层外的流体质点不断地穿入到边界层里去的缘故。总结上面所述，边界层的基本特征有：总结上面所述，边界层的基本特征有：与物体的长度相比，边界层的厚度很小；与物体的长度相比，边界层的厚度很小；边界层内沿边界层厚度的速度变化非常急剧，即速度边界层内沿边界层厚度的速度变化非常急剧，即速度梯度很大；梯度很大；边界层沿着流体流动的方向逐渐增厚；边界层沿着流体流动的方向逐渐增厚；由于边界层很薄，因而可近似地认为，边界层中各截由于边界层很薄，因而可近似地认为，边界层中各截面上的压强等于同一截面上边界层外边界上的压强；面上的压强等于同一截面上边界层外边界上的压强；在边界层内粘滞力和惯性力是同一数量级的；在边界层内粘滞力和惯性力是同一数量级的；边界层内流体的流动与管内流动一样，也可以有层流边界层内流体的流动与管内流动一样，也可以有层流和紊流两种流动状态。和紊流两种流动状态。第14页/共89页第15页/共89页对平板而言，层流转变为紊流的临界雷诺数为 边界层从层流转变为紊流的临界雷诺数的大小决定于许多因素，如前方来流的紊流度、物体壁面的粗糙度等。实验证明，增加紊流度或增加粗糙度都会使临界雷诺数值降低，即提早使层流转变为紊流。如机翼前端的边界层很薄，不大的粗糙度凸出就会透过边界层，导致层流变为紊流。第16页/共89页 现在根据边界层的特征，利用不可压缩粘性流体的运动微分方程来研究边界层内流体的运动规律。为简单起见，只讨论流体沿平板作定常的平面流动,x轴与壁面相重合，如图8-11所示。假定边界层内的流动全是层流，忽略质量力，则不可压缩粘性流体平面定常流动的微分方程和连续方程为 （8-36）第17页/共89页图8-11 推导层流边界层的微分方程用图第18页/共89页可以利用边界层每一处的厚度都很小的特征，来比较方程组（8-36）中各项的数量级，权衡主次，忽略次要项，这样便可大大简化该方程组。边界层的厚度 与平板的长度 相比较是很小的，即 ，而y的数值限制在边界层内，并满足不等式为了把方程组（8-36）变换成无量纲的，引入坐标与平板长度 、分速度与来流速度 ，压强与 之比，即引入无量纲物理量：第19页/共89页将它们代入方程组（8-36），整理后得（8-37）式中 。很显然，在边界层内，以及y与 是同一数量级，于是可取（符号表示数量级相同），所以得到如下一些数量级：第20页/共89页然后，再来求出其它各量的数量级，由连续方程因此 ，于是又得到以下数量级：为了便于讨论，将各项的数量级记载方程组（8-37）相应项的下面。现在来分析方程组（8-37）各项的数量级，以达到简化方程的目的。第21页/共89页惯性项 和 具有相同的数量级1，而惯性项和 也具有另一个相同的数量级 ，比较这两个惯性项的数量级，方程组（8-37）中第二式中各惯性项可以忽略掉。另外，比较各粘性项的数量级，可知 与 比较，可以略去；又 与 比较，可以略去；最后，比较 和 的数量级，也可以略去。于是在方程组（8-37）的粘性项中只剩第一式中的一项 根据边界层的特征，在边界层内惯性项和粘性项具有同样的数量级，由方程组（8-37）可知，必须使 和 同数量级，所以 ，即 反比于 。这表明，雷诺数越大，边界层相对厚度越小。第22页/共89页 这样，将式（8-37）中的某些项略去，再变换成有量纲量，便得到了层流边界层的微分方程（称为普朗特边界层方程）：（8-38）其边界条件为 （8-39）式中 是边界层外边界上势流的速度分布，可由势流理论来决定。对于沿平板流动，从方程组（8-38）第二式得到一个很重要的结论：在边界层内压强p与y无关，即边界层横截面上各点的压强相等，第23页/共89页 。而在边界层外边界上，边界层内的流动与外部有势流动相合。所以压强 可以根据势流的速度 由伯努力方程来决定，即因为 ，即 ，这就是说，压强项和惯性项 具有同一个数量级。对于在壁面上的各点，由式（8-38）的第一式可得（8-40）第24页/共89页 方程组（8-38）是在物体壁面为平面的假设下得到的，但是，对于曲面物体，只要壁面上任何点的曲率半径与该处边界层厚度相比很大时（机翼翼型和叶片叶型即如此），该方程组仍然是适用的，并具有足够的精确度。这时，应用曲线坐标，x轴沿着物体的曲面，y轴垂直于曲面。虽然层流边界层的微分方程（8-38）比一般的粘性流体运动微分方程要简单些，但是，即使对最简单的物体外形，这方程的求解仍是很复杂的。由于这个缘故，解决边界层问题的近似法便具有很大的实际意义。边界层的动量积分关系式为近似揭发提供了基础。第25页/共89页 在定常流动的流体中，沿边界层划出一个单位宽度的微小控制体，它的投影面ABDC（图8-12）由作为x轴的物体壁面上的一微元距离BD、边界层的外边界AC和彼此相距dx的两直线AB和CD所围成。现在应用动量方程来研究该控制体内的流体在单位时间内沿x方向的动量变化和外力之间的关系。图8-12 推导边界层的动量积分关系式用图第26页/共89页 单位时间经过AB面流入的质量和带入的动量分别为单位时间经过CD面流出的质量和带出的动量分别为对于不可压缩流体，根据连续方程从边界层外边界AC面流入的质量和带入的动量必分别为第27页/共89页式中 为边界层外边界上的速度。这样，可得单位时间沿 方向经控制面的动量通量为现在求作用在该控制体上沿x方向的一切外力。作用在AB、CD和AC诸面上的总压力沿x方向的分量分别为式中 是A与C之间的平均压强。壁面BD作用在流体上的切向应力的合力为于是，作用在该控制体上沿x方向诸外力之和为第28页/共89页其中略去了二阶微量。根据动量方程，即单位时间经控制面流体动量的通量等于外力之和，就可得到定常流动条件下卡门的边界层动量积分关系式：（8-41）由本章第四节已知，在边界层内p=p(x)；从以后的计算中可知，在给定截面上 。所以，上式两个积分都只x的函数，因此式中的偏导数可改写为全导数，上式成为在推导中对壁面上的切向应力 未作任何本质的假设，所以，式（8-41）对层流和紊流边界层都适用。边界层外边界上的速度 可以用实验或解势流问题的第29页/共89页办法求得，并可根据伯努力方程求出 的数值。所以，在边界层的动量积分关系式（8-41）中，实际可以把 和 看作已知数，而未知数只有 和 三个。因此，要解这个关系式，还需要两个补充关系式。通常把沿边界层厚度的速度分布 以及切向应力与边界层厚度的关系式 作为两个补充关系式。一般在应用边界层的动量积分关系式（8-41）来求解边界层问题时，边界层内的速度分布是按已有的经验来假定的。假定的 愈接近实际，则所得到的结果愈正确。所以，选择边界层内的速度分布函数 是求解边界层问题的重要关键。第30页/共89页平板的层流边界层的近似计算平板的层流边界层的近似计算 在实际应用中，大都采用边界层的动量积分关系式（在实际应用中，大都采用边界层的动量积分关系式（8-8-4141）对边界层进行近似计算。这方法比较简单，所得的结果也）对边界层进行近似计算。这方法比较简单，所得的结果也有足够的精确度。现在以纵向流动中的平板层流边界层为例加有足够的精确度。现在以纵向流动中的平板层流边界层为例加以说明。均匀来流速度为以说明。均匀来流速度为 的不可压缩粘性流体纵向流过一的不可压缩粘性流体纵向流过一块极薄的平板，在平板上下形成边界层。取平板的前缘点块极薄的平板，在平板上下形成边界层。取平板的前缘点o o为为坐标原点，坐标原点，x x轴沿着平轴沿着平板（即平行于板（即平行于 ），），y y轴垂直于平板，如图轴垂直于平板，如图8-158-15所示。因为顺来所示。因为顺来流方向放置的是极薄流方向放置的是极薄的平板，可以认为不的平板，可以认为不引起流动的改变。所引起流动的改变。所以，在边界层外边界以，在边界层外边界图图8-15 8-15 平板上的层流边界层平板上的层流边界层第31页/共89页上上 根据伯努力方程根据伯努力方程可知，在边界层外边界上的压强也保持常数，所以，在整个可知，在边界层外边界上的压强也保持常数，所以，在整个边界层内每一点的压强都是相同的，即边界层内每一点的压强都是相同的，即 ，这样，边界层的动量积分关系式（这样，边界层的动量积分关系式（8-418-41）变成）变成由于在上式中由三个未知数由于在上式中由三个未知数 ，所以需要再补充两，所以需要再补充两个关系式。个关系式。第一个补充关系式：假定层流边界层内的速度分布以第一个补充关系式：假定层流边界层内的速度分布以 的幂级数表示为的幂级数表示为根据下列边界条件来确定待定系数根据下列边界条件来确定待定系数 第32页/共89页第33页/共89页于是，层流边界层中速度的分布规律为于是，层流边界层中速度的分布规律为 （8-478-47）第二个补充关系式：利用牛顿内摩擦定律和式（第二个补充关系式：利用牛顿内摩擦定律和式（8-478-47）得出：）得出：（c c）为了便于计算边界层厚度，先求下列两个积分式：为了便于计算边界层厚度，先求下列两个积分式：（d d）（e e）第34页/共89页将式（将式（c c）、）、式（式（d d）和式（和式（e)e)代入边界层动量积分关系式代入边界层动量积分关系式（a a），），得得因为在平板壁面前缘点处边界层厚度为零，即因为在平板壁面前缘点处边界层厚度为零，即 ，所以积分常数所以积分常数C=0C=0。于是得边界层厚度为于是得边界层厚度为将式（将式（d d）和式（和式（e e）代入式（代入式（8-438-43）和式（）和式（8-448-44），得），得（8-488-48）（8-498-49）第35页/共89页（8-508-50）将式（将式（8-488-48）代入式（）代入式（c c），），得切向应力为得切向应力为在平板一个壁面上由粘滞力引起的总摩擦阻力为在平板一个壁面上由粘滞力引起的总摩擦阻力为摩擦阻力系数为摩擦阻力系数为（8-518-51）（8-528-52）（8-538-53）第36页/共89页平板的紊流边界层的近似计算平板的紊流边界层的近似计算 现在研究不可压缩粘性流体纵向流过平板的紊流边界层的近似计算。人们对流体在圆管中作紊流流动的规律已完整地研究过，普朗特曾经作过这样的假设：沿平板边界层内的紊流流动与管沿平板边界层内的紊流流动与管沿平板边界层内的紊流流动与管沿平板边界层内的紊流流动与管内紊流流动相同。内紊流流动相同。内紊流流动相同。内紊流流动相同。于是，就借用管内紊流流动的理论结果去找积分关系式的两个补充关系式。这时，圆管中心线上的最大速度 相当于平板的来流速度 ，圆管的半径r相当于边界层的厚度 ，并且假定平板边界层从前缘开始就是紊流。第37页/共89页 与圆管内一样，紊流边界层内速度分布的规律也假定是1/7指数规律，这与实验测得的结果很符合，于是有（a）与式（a）相应的切向应力公式为式（5-14），即（b）其中沿程损失系数 在 的范围内用勃拉休斯公式计算：将此式代入式（b），得第38页/共89页在以上雷诺数范围内，平均流速 约等于 ，将 代入上式，得现在将圆管中心线上的 和r 用边界层外边界上的和 代替，则得（c）在边界层内沿平板壁面的压强是不变的，据此，将式（a）和式（c）代入边界层的动量积分关系式（8-41），并化简整理可得：（8-54）（8-56）（8-55）第39页/共89页将式（8-54）代入式（c），得切向应力为（8-57）在平板一个壁面上，粘滞力引起的总摩擦阻力为（8-58）摩擦阻力系数为（8-59）根据实验测量，式中比较精确的系数数值是0.074，即（8-60）在推导平板紊流边界层的公式时，借用了圆管中紊流速度分布的1/7指数规律和切向应力公式，所以，以上所得结果只适用于一定的范围。实验证明，式（8-60）适用于 ；当 时，这个公式就不准确了，因为这时紊流边界层内的速度分布符合对数规律：第40页/共89页式中 ，为切向应力速度。由式（8-61）计算所得的平板紊流边界层的摩擦阻力系数 与 之间的关系画成如图8-16中的曲线3，普朗特和施利希廷（H.Schlichting）根据这条曲线写成如下的经验公式：这公式的适用范围可以达到 。后来，舒尔兹-格鲁诺（Schultz-Grunow）对平板紊流边界进行了极其细致的测量，发现在边界层内靠外侧部分的速度分布有规则地偏离于圆管内对数规律的速度分布。他根据大量实测结果提出，平板紊流边界层的摩擦阻力系数的内插公式为其相应曲线如图8-16中曲线4所示，比曲线3的偏离小。第41页/共89页第42页/共89页边界层的基本特性边界层的基本特性边界层的基本特性边界层的基本特性边界层内的流态边界层内的流态边界层内的流态边界层内的流态层流层流层流层流紊流紊流紊流紊流速度分布规律边界层厚度位移厚度动量损失厚度切向应力总摩擦阻力摩擦阻力系数表表8 81 1 边界层的近似计算公式边界层的近似计算公式 现将平板层流边界层和紊流边界层的各近似计算公式列于表8-1。第43页/共89页 从表中可以看出，平板的层流边界层和紊流边界层的重大差别有：1）紊流边界层内沿平板壁面法向截面上的速度比层流边紊流边界层内沿平板壁面法向截面上的速度比层流边紊流边界层内沿平板壁面法向截面上的速度比层流边紊流边界层内沿平板壁面法向截面上的速度比层流边界层的速度增加得快，也就是说，紊流边界层的速度分布曲线界层的速度增加得快，也就是说，紊流边界层的速度分布曲线界层的速度增加得快，也就是说，紊流边界层的速度分布曲线界层的速度增加得快，也就是说，紊流边界层的速度分布曲线要饱满得多，这与圆管中的情况相似。要饱满得多，这与圆管中的情况相似。要饱满得多，这与圆管中的情况相似。要饱满得多，这与圆管中的情况相似。2）沿平板壁面紊流边界层的厚度比层流边界层的厚度增沿平板壁面紊流边界层的厚度比层流边界层的厚度增沿平板壁面紊流边界层的厚度比层流边界层的厚度增沿平板壁面紊流边界层的厚度比层流边界层的厚度增长得快，长得快，长得快，长得快，因为紊流的 与 成比例，而层流的 则与 成比例，在紊流边界层内流体微团发生横向运动，容易促使厚度迅速增长。3）在其他条件相同的情况下，平板壁面上紊流边界层中在其他条件相同的情况下，平板壁面上紊流边界层中在其他条件相同的情况下，平板壁面上紊流边界层中在其他条件相同的情况下，平板壁面上紊流边界层中的切向应力的切向应力的切向应力的切向应力 沿着壁面的减小要层流边界层中的减小慢些。沿着壁面的减小要层流边界层中的减小慢些。沿着壁面的减小要层流边界层中的减小慢些。沿着壁面的减小要层流边界层中的减小慢些。4）在同一在同一在同一在同一 下，紊流边界层的摩擦阻力系数比层流边下，紊流边界层的摩擦阻力系数比层流边下，紊流边界层的摩擦阻力系数比层流边下，紊流边界层的摩擦阻力系数比层流边界层的大得多，界层的大得多，界层的大得多，界层的大得多，这是因为层流中的摩擦组力只是由不同流层之间发生相对运动而引起的，紊流中还有流体微团得很强烈的横向掺混，因而产生更大的摩擦阻力。第44页/共89页平板的混合边界层的近似计算平板的混合边界层的近似计算 边界层内的流动状态主要由雷诺数决定。当雷诺数增大到某一临界值时，边界层由层流转变为紊流，成为混合边界层，即平板前端是层流边界层，后部是紊流边界层，在层流转变为紊流边界层之间有一个过渡区。在大雷诺数下可以看成在某一截面上层流突然转变为紊流。由于混合边界层内的流动情况十分复杂，所以在研究平板混合边界层的摩擦阻力时，为了简化计算，作了以下两个假设（图8-17）：1）在平板的A点层流边界层突然转变为紊流边界层；2）紊流边界层的厚度变化、层流速度和切向应力分布都从前缘点O开始计算。第45页/共89页第46页/共89页 根据这两个假设，用下列方法计算平板的混合边界层的总摩擦阻力。令 代表混合边界层的总摩擦阻力，代表层流边界层的总摩擦阻力，代表紊流边界层的总摩擦阻力，则式中，为转变点A至前缘点O的距离，和 各为层流边界层和紊流边界层的摩擦阻力系数。从上式可得到混合边界层的摩擦阻力系数为式中 ，取决于层流边界层第47页/共89页A1050170033008700表82 A值与临界雷诺数 的关系转变为紊流边界层的临界雷诺数 ，见表8-2。这样，平板的混合边界层的摩擦阻力系数可按以下二式计算：当 相当于A=1700时，式（8-65）在图8-16中以曲线5表示。综上所述，层流边界层的摩擦阻力系数比紊流边界层的摩擦阻力系数要小得多，所以层流边界层段越长，即层流边界层的转变点A离平板前缘越远，平板的摩擦阻力就越小。第48页/共89页曲面边界层的分离现象曲面边界层的分离现象 如前所述，当不可压缩粘性流体纵向流过平板时，如前所述，当不可压缩粘性流体纵向流过平板时，如前所述，当不可压缩粘性流体纵向流过平板时，如前所述，当不可压缩粘性流体纵向流过平板时，在边界层外边界上沿平板方向的速度是相同的，而且整在边界层外边界上沿平板方向的速度是相同的，而且整在边界层外边界上沿平板方向的速度是相同的，而且整在边界层外边界上沿平板方向的速度是相同的，而且整个流场和边界层内的压强都保持不变。当粘性流体流经个流场和边界层内的压强都保持不变。当粘性流体流经个流场和边界层内的压强都保持不变。当粘性流体流经个流场和边界层内的压强都保持不变。当粘性流体流经曲面物体时，边界层外边界上沿曲面方向的速度曲面物体时，边界层外边界上沿曲面方向的速度曲面物体时，边界层外边界上沿曲面方向的速度曲面物体时，边界层外边界上沿曲面方向的速度 是改是改是改是改变的，所以曲面边界层内的压强也将同样发生变化，对变的，所以曲面边界层内的压强也将同样发生变化，对变的，所以曲面边界层内的压强也将同样发生变化，对变的，所以曲面边界层内的压强也将同样发生变化，对边界层内的流动将产生影响。曲面边界层的计算是很复边界层内的流动将产生影响。曲面边界层的计算是很复边界层内的流动将产生影响。曲面边界层的计算是很复边界层内的流动将产生影响。曲面边界层的计算是很复杂的，这里不准备讨论它。这一节将着重说明曲面边界杂的，这里不准备讨论它。这一节将着重说明曲面边界杂的，这里不准备讨论它。这一节将着重说明曲面边界杂的，这里不准备讨论它。这一节将着重说明曲面边界层的分离现象。层的分离现象。层的分离现象。层的分离现象。第49页/共89页图818 机翼翼型上的边界层边界层外边界M流体经过机翼翼型（或叶片叶型）的流动如图流体经过机翼翼型（或叶片叶型）的流动如图流体经过机翼翼型（或叶片叶型）的流动如图流体经过机翼翼型（或叶片叶型）的流动如图8-188-18所示，所示，所示，所示，以以以以 和和和和 表示无穷远处流体流动所具有的速度和压强。表示无穷远处流体流动所具有的速度和压强。表示无穷远处流体流动所具有的速度和压强。表示无穷远处流体流动所具有的速度和压强。流体绕过翼型前驻点后，沿上表面的流速先增加，直增流体绕过翼型前驻点后，沿上表面的流速先增加，直增流体绕过翼型前驻点后，沿上表面的流速先增加，直增流体绕过翼型前驻点后，沿上表面的流速先增加，直增加到曲面上某一点加到曲面上某一点加到曲面上某一点加到曲面上某一点MM，然后降低。由伯努利方程可知，然后降低。由伯努利方程可知，然后降低。由伯努利方程可知，然后降低。由伯努利方程可知，相应的压强先降低（相应的压强先降低（相应的压强先降低（相应的压强先降低（dp/dx0dp/dx0)(dp/dx0)。MM点处边界层外边界上的速度最大，而压强最低。沿点处边界层外边界上的速度最大，而压强最低。沿点处边界层外边界上的速度最大，而压强最低。沿点处边界层外边界上的速度最大，而压强最低。沿曲面各点法向的速度剖面和压强变化曲线的示意如图曲面各点法向的速度剖面和压强变化曲线的示意如图曲面各点法向的速度剖面和压强变化曲线的示意如图曲面各点法向的速度剖面和压强变化曲线的示意如图8-8-1919所示。图中实线表示流线，虚线表示边界层的外边界。所示。图中实线表示流线，虚线表示边界层的外边界。所示。图中实线表示流线，虚线表示边界层的外边界。所示。图中实线表示流线，虚线表示边界层的外边界。第50页/共89页图819 曲面边界层分离的形成示意图第51页/共89页第52页/共89页第53页/共89页 首先从流体在边界层内流动的物理过程来说明曲面首先从流体在边界层内流动的物理过程来说明曲面首先从流体在边界层内流动的物理过程来说明曲面首先从流体在边界层内流动的物理过程来说明曲面边界层的分离现象。当粘性流体流经曲面时，边界层内边界层的分离现象。当粘性流体流经曲面时，边界层内边界层的分离现象。当粘性流体流经曲面时，边界层内边界层的分离现象。当粘性流体流经曲面时，边界层内的流体微团被粘滞力阻滞，动能损耗越大，减速也越甚。的流体微团被粘滞力阻滞，动能损耗越大，减速也越甚。的流体微团被粘滞力阻滞，动能损耗越大，减速也越甚。的流体微团被粘滞力阻滞，动能损耗越大，减速也越甚。在曲面的降压加速段中，由于流体的部分压强势能转变在曲面的降压加速段中，由于流体的部分压强势能转变在曲面的降压加速段中，由于流体的部分压强势能转变在曲面的降压加速段中，由于流体的部分压强势能转变为流体的动能，故流体微团虽然受到粘滞力的阻滞作用，为流体的动能，故流体微团虽然受到粘滞力的阻滞作用，为流体的动能，故流体微团虽然受到粘滞力的阻滞作用，为流体的动能，故流体微团虽然受到粘滞力的阻滞作用，但仍有足够的动能使它继续前进。但是，但仍有足够的动能使它继续前进。但是，但仍有足够的动能使它继续前进。但是，但仍有足够的动能使它继续前进。但是，在曲面在曲面的升压的升压的升压的升压减速段中，流体的部分动能不仅要转变为压强势能，而减速段中，流体的部分动能不仅要转变为压强势能，而减速段中，流体的部分动能不仅要转变为压强势能，而减速段中，流体的部分动能不仅要转变为压强势能，而且粘滞力的阻滞作用也要继续损耗它的动能，这就使流且粘滞力的阻滞作用也要继续损耗它的动能，这就使流且粘滞力的阻滞作用也要继续损耗它的动能，这就使流且粘滞力的阻滞作用也要继续损耗它的动能，这就使流体微团的动能损耗加大，流速迅速降低，其边界层不断体微团的动能损耗加大，流速迅速降低，其边界层不断体微团的动能损耗加大，流速迅速降低，其边界层不断体微团的动能损耗加大，流速迅速降低，其边界层不断增厚。当流体流到曲面的某一点增厚。当流体流到曲面的某一点增厚。当流体流到曲面的某一点增厚。当流体流到曲面的某一点S S时，靠近物体壁面的时，靠近物体壁面的时，靠近物体壁面的时，靠近物体壁面的流体微团的动能已被消耗尽，这部分流体微团便停滞不流体微团的动能已被消耗尽，这部分流体微团便停滞不流体微团的动能已被消耗尽，这部分流体微团便停滞不流体微团的动能已被消耗尽，这部分流体微团便停滞不前。跟着而来的流体微团也将同样停滞下来，以致越来前。跟着而来的流体微团也将同样停滞下来，以致越来前。跟着而来的流体微团也将同样停滞下来，以致越来前。跟着而来的流体微团也将同样停滞下来，以致越来越多的被停滞的流体微团在物体壁面和主流之间堆积起越多的被停滞的流体微团在物体壁面和主流之间堆积起越多的被停滞的流体微团在物体壁面和主流之间堆积起越多的被停滞的流体微团在物体壁面和主流之间堆积起来。与此同时，在来。与此同时，在来。与此同时，在来。与此同时，在S S点之后，压强的继续升高将使部分点之后，压强的继续升高将使部分点之后，压强的继续升高将使部分点之后，压强的继续升高将使部分流体微团被迫反向逆流，并迅速向外扩展，造成边界层流体微团被迫反向逆流，并迅速向外扩展，造成边界层流体微团被迫反向逆流，并迅速向外扩展，造成边界层流体微团被迫反向逆流，并迅速向外扩展，造成边界层的分离。在的分离。在的分离。在的分离。在STST线上一系列流体微团的切向速度等线上一系列流体微团的切向速度等线上一系列流体微团的切向速度等线上一系列流体微团的切向速度等第54页/共89页于零，于零，于零，于零，S S点称为边界层的分离点。分离时形成的漩涡，点称为边界层的分离点。分离时形成的漩涡，点称为边界层的分离点。分离时形成的漩涡，点称为边界层的分离点。分离时形成的漩涡，不断地被主流带走，在物体后部形成尾涡区。在渐缩渐不断地被主流带走，在物体后部形成尾涡区。在渐缩渐不断地被主流带走，在物体后部形成尾涡区。在渐缩渐不断地被主流带走，在物体后部形成尾涡区。在渐缩渐扩管的渐扩段中，或其他形式的渐扩管中，也同样有可扩管的渐扩段中，或其他形式的渐扩管中，也同样有可扩管的渐扩段中，或其他形式的渐扩管中，也同样有可扩管的渐扩段中，或其他形式的渐扩管中，也同样有可能出现边界层的分离现象。能出现边界层的分离现象。能出现边界层的分离现象。能出现边界层的分离现象。从以上的分析可得如下结论：粘性流体在压强降低从以上的分析可得如下结论：粘性流体在压强降低从以上的分析可得如下结论：粘性流体在压强降低从以上的分析可得如下结论：粘性流体在压强降低区内流动（加速流动）时，不会出现边界层分离，区内流动（加速流动）时，不会出现边界层分离，区内流动（加速流动）时，不会出现边界层分离，区内流动（加速流动）时，不会出现边界层分离，只有只有只有只有在压强升高区内流动（减速流动）时，才有可能出现分在压强升高区内流动（减速流动）时，才有可能出现分在压强升高区内流动（减速流动）时，才有可能出现分在压强升高区内流动（减速流动）时，才有可能出现分离，形成漩涡。离，形成漩涡。离，形成漩涡。离，形成漩涡。尤其在主流的减速足够大的情况下，边尤其在主流的减速足够大的情况下，边尤其在主流的减速足够大的情况下，边尤其在主流的减速足够大的情况下，边界层的分离就一定会发生。例如，在圆柱体和球体这样界层的分离就一定会发生。例如，在圆柱体和球体这样界层的分离就一定会发生。例如，在圆柱体和球体这样界层的分离就一定会发生。例如，在圆柱体和球体这样的钝头体的后半部分，当流速足够大时，便会发生边界的钝头体的后半部分，当流速足够大时，便会发生边界的钝头体的后半部分，当流速足够大时，便会发生边界的钝头体的后半部分，当流速足够大时，便会发生边界层的分离，这是由于在钝头体的后半部分有急剧的压强层的分离，这是由于在钝头体的后半部分有急剧的压强层的分离，这是由于在钝头体的后半部分有急剧的压强层的分离，这是由于在钝头体的后半部分有急剧的压强升高区，主流减速加剧的缘故。若将钝头体的后半部分升高区，主流减速加剧的缘故。若将钝头体的后半部分升高区，主流减速加剧的缘故。若将钝头体的后半部分升高区，主流减速加剧的缘故。若将钝头体的后半部分改为充分细长形的尾部，成为圆头尖尾的所谓流线型物改为充分细长形的尾部，成为圆头尖尾的所谓流线型物改为充分细长形的尾部，成为圆头尖尾的所谓流线型物改为充分细长形的尾部，成为圆头尖尾的所谓流线型物体（如叶片叶型和机翼翼型），就可使主流的减速大为体（如叶片叶型和机翼翼型），就可使主流的减速大为体（如叶片叶型和机翼翼型），就可使主流的减速大为体（如叶片叶型和机翼翼型），就可使主流的减速大为降低，防止边界层内逆流的发生，避免边界层的分离。降低，防止边界层内逆流的发生，避免边界层的分离。降低，防止边界层内逆流的发生，避免边界层的分离。降低，防止边界层内逆流的发生，避免边界层的分离。第55页/共89页第56页/共89页绕过圆柱体的流动绕过圆柱体的流动 卡门涡街卡门涡街 为了进一步说明边界层分离这一重要现象，可为了进一步说明边界层分离这一重要现象，可为了进一步说明边界层分离这一重要现象，可为了进一步说明边界层分离这一重要现象，可考察粘性流体绕过圆柱体的流动情况。把一个圆柱考察粘性流体绕过圆柱体的流动情况。把一个圆柱考察粘性流体绕过圆柱体的流动情况。把一个圆柱考察粘性流体绕过圆柱体的流动情况。把一个圆柱体放在静止的流体中，然后流体以相当于几个雷诺体放在静止的流体中，然后流体以相当于几个雷诺体放在静止的流体中，然后流体以相当于几个雷诺体放在静止的流体中，然后流体以相当于几个雷诺数（数（数（数（）的很低的速度）的很低的速度）的很低的速度）的很低的速度 绕流它。在开始绕流它。在开始绕流它。在开始绕流它。在开始瞬间与理想流体绕流圆柱体一样，流体在前驻点速瞬间与理想流体绕流圆柱体一样，流体在前驻点速瞬间与理想流体绕流圆柱体一样，流体在前驻点速瞬间与理想流体绕流圆柱体一样，流体在前驻点速度为零，而后沿圆柱体左右两侧流动，流动在圆柱度为零，而后沿圆柱体左右两侧流动，流动在圆柱度为零，而后沿圆柱体左右两侧流动，流动在圆柱度为零，而后沿圆柱体左右两侧流动，流动在圆柱体的前半部分是降压，速度逐渐增大到最大值，在体的前半部分是降压，速度逐渐增大到最大值，在体的前半部分是降压，速度逐渐增大到最大值，在体的前半部分是降压，速度逐渐增大到最大值，在后半部分是升压，速度逐渐下降，到后驻点重新等后半部分是升压，速度逐渐下降，到后驻点重新等后半部分是升压，速度逐渐下降，到后驻点重新等后半部分是升压，速度逐渐下降，到后驻点重新等于零（图于零（图于零（图于零（图8-218-21，a a）。）。）。）。以后逐渐增大来流速度，也即以后逐渐增大来流速度，也即以后逐渐增大来流速度，也即以后逐渐增大来流速度，也即增大雷诺数，使圆柱体后半部分的压强梯度增加，增大雷诺数，使圆柱体后半部分的压强梯度增加，增大雷诺数，使圆柱体后半部分的压强梯度增加，增大雷诺数，使圆柱体后半部分的压强梯度增加，以致引起边界层的分离（图以致引起边界层的分离（图以致引起边界层的分离（图以致引起边界层的分离（图8-218-21，b b）。）。）。）。随着来流雷随着来流雷随着来流雷随着来流雷诺数的不断增加，圆柱体后半部分边界层中的流体诺数的不断增加，圆柱体后半部分边界层中的流体诺数的不断增加，圆柱体后半部分边界层中的流体诺数的不断增加，圆柱体后半部分边界层中的流体微团受到更大的阻滞，分离点微团受到更大的阻滞，分离点微团受到更大的阻滞，分离点微团受到更大的阻滞，分离点S S一直向前移动。一直向前移动。一直向前移动。一直向前移动。第57页/共89页第58页/共89页当雷诺数增加到大约当雷诺数增加到大约当雷诺数增加到大约当雷诺数增加到大约4040时，在圆柱体的后面便产生时，在圆柱体的后面便产生时，在圆柱体的后面便产生时，在圆柱体的后面便产生一对旋转方向相反的对称漩涡（图一对旋转方向相反的对称漩涡（图一对旋转方向相反的对称漩涡（图一对旋转方向相反的对称漩涡（图8-218-21，c c）。）。）。）。雷诺雷诺雷诺雷诺数超过数超过数超过数超过4040后，对称漩涡不断增长并出现摆动，直到后，对称漩涡不断增长并出现摆动，直到后，对称漩涡不断增长并出现摆动，直到后，对称漩涡不断增长并出现摆动，直到 时，时，时，时，这对不稳定的对称漩涡分裂，最后这对不稳定的对称漩涡分裂，最后这对不稳定的对称漩涡分裂，最后这对不稳定的对称漩涡分裂，最后形成几乎稳形成几乎稳形成几乎稳形成几乎稳定的、非对称性的、多少有些规则的、旋转方向相定的、非对称性的、多少有些规则的、旋转方向相定的、非对称性的、多少有些规则的、旋转方向相定的、非对称性的、多少有些规则的、旋转方向相反的交替漩涡，称为卡门涡街反的交替漩涡，称为卡门涡街反的交替漩涡，称为卡门涡街反的交替漩涡，称为卡门涡街（图（图（图（图8-218-21，d d）。）。）。）。它以它以它以它以比来流速度比来流速度比来流速度比来流速度 小得多的速度小得多的速度小得多的速度小得多的速度 运动。图运动。图运动。图运动。图8-228-22所示为所示为所示为所示为流体以不同雷诺数绕过圆柱体的流动情况，从中可流体以不同雷诺数绕过圆柱体的流动情况，从中可流体以不同雷诺数绕过圆柱体的流动情况，从中可流体以不同雷诺数绕过圆柱体的流动情况，从中可以清楚地看出圆柱体后面的尾流中一对漩涡的形成以清楚地看出圆柱体后面的尾流中一对漩涡的形成以清楚地看出圆柱体后面的尾流中一对漩涡的形成以清楚地看出圆柱体后面的尾流中一对漩涡的形成过程。图过程。图过程。图过程。图8-238-23所示为在椭圆柱体后形成的卡门涡街。所示为在椭圆柱体后形成的卡门涡