1、Folie1Folie2边界层理论及其近似边界层理论及其近似5.15.1、边界层近似及其特征、边界层近似及其特征5.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程5.35.3、平板层流边界层的相似解、平板层流边界层的相似解5.45.4、边界层动量积分方程、边界层动量积分方程5.55.5、边界层的分离现象、边界层的分离现象Folie3 5.15.1、边界层近似及其特征、边界层近似及其特征 1、边界层概念的提出 业已知道，流动Re数（O.Reynolds，1883年，英国流体力学家）是用以表征流体质点的惯性力与粘性力对比关系的。根据量级分析，作用于流体上的惯性力和粘性力可
2、表示为:惯性力：粘性力：惯性力/粘性力：20203ULtULdtdVmFJLUAdydVF0Re00202LULUULFFJFolie4 5.15.1、边界层近似及其特征、边界层近似及其特征 因此，在高Re数下，流体运动的惯性力远远大于粘性力。这样研究忽略粘性力的流动问题是有实际意义的。这也是早期发展理想流体力学的重要依据，而且确实较成功地解决了与粘性关系不大的一系列流动问题，诸如绕流物体的升力、波动等问题，但对绕流物体阻力、涡的扩散等问题，理想流体力学的解与实际相差甚远，且甚至得出完全相反的结论，圆柱绕流无阻力的DAlembert疑题就是一个典型的例子。（DAlembert，法国力学家，17
3、17-1783）Folie5 5.15.1、边界层近似及其特征、边界层近似及其特征 那么，如何考虑流体的粘性，怎样解决绕流物体的阻力问题，这在当时确实是一个阻碍流体力学发展的难题，直到1904年国际流体力学大师德国学者L.Prandtl通过大量实验发现，虽然整体流动的Re数很大，但在靠近物面的薄层流体内，流场的特征与理想流动相差甚远，沿着法向存在很大的速度梯度，粘性力无法忽略。Prandtl把这一物面近区粘性力起重要作用的薄层称为边界层（Boundary layer）。Folie6Folie7 5.15.1、边界层近似及其特征、边界层近似及其特征Prandtl的边界层概念，为人们如何计入粘性的
4、作用开辟了划时代的途径，因此称其为粘性流体力学之父。对整个流场提出的基本分区是：（1 1）整个流动区域可分成理想流体的流动区域（）整个流动区域可分成理想流体的流动区域（势流区）和粘性流体的流动区域（粘流区）。势流区）和粘性流体的流动区域（粘流区）。（2 2）在远离物体的理想流体流动区域，可忽略粘）在远离物体的理想流体流动区域，可忽略粘性的影响，按势流理论处理。性的影响，按势流理论处理。（3 3）粘性流动区域仅限于物面近区的薄层内，称）粘性流动区域仅限于物面近区的薄层内，称为边界层。既然是粘流区，粘性力的作用不能忽为边界层。既然是粘流区，粘性力的作用不能忽略，与惯性力同量级，流体质点作有旋运动。
5、略，与惯性力同量级，流体质点作有旋运动。Folie8 5.15.1、边界层近似及其特征、边界层近似及其特征 2、边界层的特征（1）边界层定义 严格而言，边界层区与主流区之间无明显界线，通常以速度达到主流区速度的0.99U作为边界层的外缘。由边界层外缘到物面的垂直距离称为边界层名义厚度。（2）边界层的有涡性 粘性流体运动总伴随涡量的产生、扩散、衰减。边界层就是涡层，当流体绕过物面时，无滑移边界条件相当于使物面成为具有一定强度的连续分布的涡源。以二维流动为例说明之。Folie9 5.15.1、边界层近似及其特征、边界层近似及其特征此时，物面上的涡源强度为对于不可压缩流体，二维流动的涡量输运方程为
6、上式表明，由于粘性的影响，物面上的涡量一方面沿垂直流线方向扩散，另一方面，涡量沿主流方向迁移，并随之而逐渐衰减。涡量的扩散速度与粘性有关，涡量的迁移速度取决于流动速度。oxxyzyuyuxu2222yxdtdzzzzFolie10 边界层边界层 拟序结构产生及发展的过程拟序结构产生及发展的过程T-S 波L L涡涡发卡涡发卡涡变形、变形、拉伸、拉伸、诱发诱发湍流斑湍流斑动画演示：拟序涡的生成与发展动画演示：拟序涡的生成与发展 (Ma=0.7,(Ma=0.7,平板）平板）Folie11 边界层边界层 拟序结构产生及发展的过程拟序结构产生及发展的过程动画演示:发卡涡的演化过程（Ma=0.7,平板）M
7、ovie:Isosurface of Q=10Folie12 边界层边界层 拟序结构产生及发展的过程拟序结构产生及发展的过程不同时刻展向涡量的分布不同时刻展向涡量的分布 z=0.044z=0.044 发卡涡导致高剪切发卡涡导致高剪切诱发转捩诱发转捩Folie13 5.15.1、边界层近似及其特征、边界层近似及其特征（3）边界层厚度的量级估计 根据边界层内粘性力与惯性力同量级的条件，可估算边界层的厚度。以平板绕流为例说明。设来流的速度为U，在x方向的长度为L，边界层厚度为。惯性力：粘性力：由惯性力与粘性力同量级得到 22ULtULdtdVmFJ2LUAdydVFRe1 UL 22LLUFFJFo
8、lie14 5.15.1、边界层近似及其特征、边界层近似及其特征由此可见，在高Re数下，边界层的厚度远小于被绕流物体的特征长度。（4）边界层各种厚度定义（a）边界层排移厚度 在边界层内，理想流体的质量流量为 其中，ue为边界层外缘速度。0dyumeeiFolie15 5.15.1、边界层近似及其特征、边界层近似及其特征由于粘性的存在，实际流体通过的质量流量为 上述两项之差表示粘性存在而损失的流量，这部分流量被排挤到主流场中，相当于主流区增加了一层流体。0dyumeei0udymeFolie16 5.15.1、边界层近似及其特征、边界层近似及其特征 主流区所增加的厚度为 这部分主流区增加的流体厚
9、度是由边界层流体排挤入主流区造成的。因此，称其为排移厚度。（b）边界层动量损失厚度 在边界层内，在质量流量不变的条件下，理想流体通过的动量为 由于粘性的存在，实际流体通过的动量为01dyuuueeee011dyuuee0udyuKei02dyuKeFolie17 5.15.1、边界层近似及其特征、边界层近似及其特征 上述两项之差表示粘性存在而损失的动量，这部分动量损失用外流流速ue（理想流体）折算的动量损失厚度为（c）边界层能量损失厚度 在边界层内，在质量流量不变的条件下，以外流速度（理想流体）通过的动能为 由于粘性的存在，实际流体通过的动能为022dyuuuuueee021dyuuuueee
10、0221udyuEei0221udyuEeFolie18 5.15.1、边界层近似及其特征、边界层近似及其特征 上述两项之差表示粘性存在而损失的动能，上述两项之差表示粘性存在而损失的动能，这部分动能损失用主流流速这部分动能损失用主流流速u ue e（理想流体）折算（理想流体）折算的动能损失厚度为的动能损失厚度为:上述各种厚度的计算公式，对于不可压缩流上述各种厚度的计算公式，对于不可压缩流体而言，变为体而言，变为:032322121dyuuuuueeee02231dyuuuueee011dyuue021dyuuuuee02231dyuuuueeFolie19 5.15.1、边界层近似及其特征、边
11、界层近似及其特征（5）几点说明（a）实际流动中，边界层流动与理想流动是渐近过渡的，边界层的外边界线实际上是不存在的，因此边界层的外边界线不是流线，而是被流体所通过的，允许流体穿过边界线流动。在边界层内流线是向外偏的。Folie20 5.15.1、边界层近似及其特征、边界层近似及其特征（b）边界层各种厚度的定义式，即适用于层流，也适用于湍流。（c）边界层各种厚度的大小与边界层内流速分布有关。但各厚度的大小依次是:边界层厚度边界层排移厚度边界层动量损失厚度Folie21Ludwig PrandtlLudwig Prandtl介绍介绍 普朗特重视观察和分析力学现象，养成了非凡的直观洞察能力，善于抓住
12、物理本质，概括出数学方程。他曾说：“我只是在相信自己对物理本质已经有深入了解以后，才想到数学方程。方程的用处是说出量的大小，这是直观得不到的，同时它也证明结论是否正确。”普朗特指导过81名博士生，著名学者Blasius、Von Karman是其学生之一。我国著名的空气动力学专家、北航流体力学教授陆士嘉先生（女，19111986）是普朗特正式接受的唯一中国学生，唯一的女学生。郭永怀（1909-1968），冯卡门学生，空气动力学家，“两弹一星功勋奖章”。Folie225.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程Folie235.25.2、平面不可压缩流体层流边界层方
13、程、平面不可压缩流体层流边界层方程1、平壁面上边界层方程 根据Prandtl边界层概念，通过量级比较，可对N-S方程组进行简化，得到边界层近似方程。对于二维不可压缩流动，N-S方程为 选取长度特征L，速度尺度ue，时间尺度t=L/ue，边界层近似假定：0yvxu22221yuxuxpfyuvxuutux22221yvxvypfyvvxvutvyFolie245.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程（1）根据边界层定义，纵向偏导数远远小于横向偏导数。（2）法向速度远远小于纵向速度。（3）边界层内的压强与外流速度的平方成正比。将这些量级关系式代入到N-S方程中，
14、得到yxyLxLL,1,1,Re1uvutLuuvuLuLtveeee,Re1,/2eupFolie255.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程N-S方程组与各项量级比较:LuLLuyvxueeReReu1Reu 0eeLLuLuLuLuyuxuxpfyuvxuutueeeex2e2e2e2e2e222222uu LReuLu Re/Reu 1ReLuReLu ReLuReLu Re Re Re12e2e2/32e2e2222222LuLuLuyvxvypfyvvxvutveeeyFolie265.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层
15、流边界层方程在高Re数情况下，忽略小量得到忽略质量力，由第三个方程得到这说明，在高Re数情况下，在边界层内压力沿法向是不变的。0yvxu221yuxpfyuvxuutuxypfy10yp10Folie275.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程 边界层内的压力分布与边界层外边界线上的压力分布相等。也就是，p与y无关，仅是x和t的函数。即 忽略质量力，Prandtl边界层方程变为 边界条件：),(txppe0yvxu221yuxpyuvxuutuyp10euu y0 v0u 0yFolie285.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边
16、界层方程在边界层外边界线上，可按照理想流体势流方程在边界层外边界线上，可按照理想流体势流方程确定压强。即确定压强。即在定常流动情况下，有在定常流动情况下，有xpxuutueeee10yvxu22yuxuutuyuvxuutueee0yvxu22yuxuuyuvxuueeFolie295.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程综上所述，边界层基本特性可归纳为),(pp 0 Re1uv Re1eetxypyxLFolie305.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程2、曲壁面上的边界层方程 在实际流动中所遇到的物面常是弯曲的，因此
17、推导曲壁面上的边界层方程具有重要意义。在推导中，使用曲壁面上的边界层坐标系。其中，x轴贴着壁面，y轴垂直于壁面。在边界层内任取一点M，其坐标 x=ON y=NM M为M的邻点，MM的弧长为ds22)()(MMMMdsFolie315.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程 在x处，设曲壁的曲率半径为R(x)，有 则 仍以u和v分别表示边界层坐标系中的x和y方向的速度分量，则由正交曲线坐标系方程，得到连续方程dyMNNMdxRyRdyRMMMM ,)(1,RyR)()(RyR )()(212221222222HHdyHdxHdydxMMMMds0yRvyvxuR
18、yRFolie325.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程运动方程为:xudxdRyRRydxdRyRRvxvyRRyRuyuyRyuxuyRRxpyRRyRuvyuvxuuyRRtu332222222)()()(2)(11xvdxdRyRRydxdRyRRuxuyRRyRvyvyRyvxvyRRypyRuyvvxvuyRRtv3322222222)()()(2)(11Folie335.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程 假定物面的曲率半径R(x)与x向的特征长度L同量级，y的量级与边界层厚度同量级，故有 量级比较，简化
19、的边界层方程为Re1,1,1,11,1,1,1RyyxyLxRyRRyHdxddxdRRLuvutLuuvuLuLtveeee,Re1,/0yvxu221yuxpyuvxuutuypRu2Folie345.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程 这就是曲壁面上的边界层方程，与平壁面的方程相比，只是y方向的方程有所不同。为了和流动弯曲所产生的离心力相平衡，必须有y方向的压力梯度。以下估计这个压力梯度的量级大小。初步假定边界层内速度分布为线性分布。从y=0到y=s积分，有 2222,yRuyuRypyuueeeRRuppe31up ,31)0()(2e2Folie
20、355.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程 在Rs的情况下，此压差是个小量，可忽略不计。由此仍得出在曲壁面的边界层内，法向压力不变是个常数。这说明，在曲率半径不太小且变化不太大的情况下，曲壁面上的边界层方程与平壁面上的边界层方程完全相同。Folie365.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程Folie375.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程Folie385.25.2、平面不可压缩流体层流边界层方程、平面不可压缩流体层流边界层方程Folie395.35.3、平板层流边界层的相似解、平
21、板层流边界层的相似解 1908年，Prandtl学生Blasius利用边界层速度分布的相似性求解了平板层流边界层方程。对于零压梯度、定常、不可压缩流体平板层流绕流，边界层方程为相应的边界条件为Blasius假设，在平板上边界层内的速度分布具有相似性特征。即 0yvxu22yuyuvxuu0 0 0 ;yuvyuU()euyguFolie405.35.3、平板层流边界层的相似解、平板层流边界层的相似解根据量级比较，边界层厚度的量级为:引入流函数：根据相似假设，流函数的形式表示为：RexxxxUU xUyyx()xUf,uvyx Folie415.35.3、平板层流边界层的相似解、平板层流边界层的
22、相似解Folie425.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie435.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie445.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie45Folie465.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie475.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie485.35.3、平板层流边界层的相似解、平板层流边界层的相似解Folie495.45.4、边界层动量积分方程、边界层动量积分方程 边界层动量积分关系式是由Karman1921导出的，对近似求解边界层
23、特性具有重要作用。适应于层流边界层和湍流边界层。今在边界层内任取一控制体，控制体长度为dx，控制面为Aab、Abc、Acd、Ada。现对控制体应用动量定律，可得 由Aab面流入控制体的质量为 由Acd面流出控制体的质量为)(0 xabudymdxudyxmmxabcd)(0Folie505.45.4、边界层动量积分方程、边界层动量积分方程根据质量守恒定律，通过Abc流入控制体的质量为由Aab面流入控制体的动量为由Acd面流出控制体的动量为通过Abc流入控制体的动量在x方向的分量为 dxudyxmmmxabcdbc)(0)(02xabdyuKdxdyuxKKxabcd)(02dxudyxuKxe
24、bc)(0Folie515.45.4、边界层动量积分方程、边界层动量积分方程在Aab面上的作用力为在Acd面上的作用力为在Abc面上的力为在Aad面上的切应力为)(xpFeab)(dxdxdxdppFeecdddxdxdppFeebc2dxFad0Folie525.45.4、边界层动量积分方程、边界层动量积分方程现对控制体建立x方向的动量方程为整理后，得dxudyxuKdxdyuxKdxddxdxdppdxdxdxdppxpxeabxabeeeee)(0)(0202)()()(02)(00)(xxeedyudxdudydxduxdxdpFolie535.45.4、边界层动量积分方程、边界层动量
25、积分方程Folie545.45.4、边界层动量积分方程、边界层动量积分方程Folie555.45.4、边界层动量积分方程、边界层动量积分方程Folie565.45.4、边界层动量积分方程、边界层动量积分方程Folie575.45.4、边界层动量积分方程、边界层动量积分方程Folie585.45.4、边界层动量积分方程、边界层动量积分方程Folie595.45.4、边界层动量积分方程、边界层动量积分方程Folie605.45.4、边界层动量积分方程、边界层动量积分方程Folie615.45.4、边界层动量积分方程、边界层动量积分方程Folie625.45.4、边界层动量积分方程、边界层动量积分方
26、程Folie635.45.4、边界层动量积分方程、边界层动量积分方程Folie645.45.4、边界层动量积分方程、边界层动量积分方程Folie655.55.5、边界层的分离现象、边界层的分离现象Folie665.55.5、边界层的分离现象、边界层的分离现象Folie675.55.5、边界层的分离现象、边界层的分离现象从静止开始边界从静止开始边界层发展情况层发展情况扩张管扩张管（上壁有抽吸）（上壁有抽吸）Folie685.55.5、边界层的分离现象、边界层的分离现象Folie695.55.5、边界层的分离现象、边界层的分离现象Folie705.55.5、边界层的分离现象、边界层的分离现象1、边
27、界层分离现象 边界层中的流体质点受惯性力、粘性力和压力的作用。其中，粘性力的作用始终是阻滞流体质点运动，使流体质点减速，失去动能；压力的作用取决于绕流物体的形状和流道形状，顺压梯度有助于流体加速前进，而逆压梯度阻碍流体运动。以圆柱绕流为例说明边界层的分离现象。对于理想流体，流体微团绕过圆柱时，在OM段为加速减压区，压能转化为动能。在MF段为减速增压区，动能减小压能增加Folie715.55.5、边界层的分离现象、边界层的分离现象 对于粘性流体，在上述能量的转化过程中，由于粘性的作用，边界层内的流体质点将要克服粘性力作功而消耗机械能。因此微团在逆压区，不可能到达F点，而是在MF段中的某点处微团速
28、度降为零，以后来的质点将改道进入主流中，使来流边界层与壁面分离。在分离点下游的区域，受逆压梯度的作用而发生倒流。分离点定义为紧邻壁面顺流区与倒流区的分界点。在分离点附近和分离区，由于边在分离点附近和分离区，由于边界层厚度大界层厚度大 大增加，边界层假大增加，边界层假设不在成立。边界层分离的必要设不在成立。边界层分离的必要条件是：条件是：逆压梯度逆压梯度和和物面粘性物面粘性的的阻滞作用结果。阻滞作用结果。Folie725.55.5、边界层的分离现象、边界层的分离现象 仅有粘性的阻滞作用而无逆压梯度，不会发生边界层的分离，因为无反推力使边界层流体进入到外流区。这说明，顺压梯度的流动不可能发生边界层
29、分离。只有逆压梯度而无粘性的阻滞作用，同样也不会发生分离现象，因为无阻滞作用，运动流体不可能消耗动能而滞止下来。Folie735.55.5、边界层的分离现象、边界层的分离现象气流绕翼型的流动与边界层分离现象。需要指出的是：逆压梯度和壁面粘性阻滞作用是边界层分离的必要条件，但不是充分的，也就是说只有在一定的逆压梯度下，才有可能发生分离。Folie74Folie75Folie765.55.5、边界层的分离现象、边界层的分离现象2、在不同压力梯度区边界层的速度分布特征根据边界层动量方程，在壁面上压力梯度对边界层内流动速度分布产生一定的影响。对于顺压梯度的情况，有对于逆压梯度的情况，有xuuxpyue
30、e1220,0,022xuyuxpe0,0,022xuyuxpeFolie775.55.5、边界层的分离现象、边界层的分离现象对于零压梯度的情况，有由此可见，随着压力梯度的变号，边界层速度分布的曲率将改变符号。0,0,022xuyuxpeFolie785.55.5、边界层的分离现象、边界层的分离现象 对于顺压梯度区，压力沿程减小，速度沿程增加。在壁面处，另一方面，在边界层的外边界上，有 由此说明，在顺压梯度区，边界层内的速度沿y方向是单调增加的，分布曲线无拐点，是一条向外凸的光滑曲线，流动是稳定的。0,00220yyyuyu0,0022yyyuyuFolie795.55.5、边界层的分离现象、
31、边界层的分离现象 对于逆压梯度区（分离点前），压力沿程增加，速度沿程减小。在壁面处，有另一方面，在边界层的外边界上0,00220yyyuyu0,022yyyuyuFolie805.55.5、边界层的分离现象、边界层的分离现象 于是在边界层内 ，速度分布曲率从正变为负，在某点处必然有 这一点是速度分布的拐点。拐点的出现改变了速度分布的形状，在拐点以上为外凸型，在拐点以下为外凹型，存在拐点的速度分布型是不稳定的。0022yyyu y0Folie815.55.5、边界层的分离现象、边界层的分离现象 在最小压力点处，有 说明拐点在物面上，随着流体质点向下游流动，拐点向外边界移动，物面近区的速度分布愈来
32、愈瘦小，但当拐点移动到某点时，物面处出现 该点称为分离点。0,0022yyuxdp0,00220yyyuyuFolie825.55.5、边界层的分离现象、边界层的分离现象 在分离点下游区，有 0,0,00022yyyuyuxdpFolie835.55.5、边界层的分离现象、边界层的分离现象 发生了回流，回流把主流推离壁面，边界层假设失效。由上分析可见，逆压梯度愈大，边界层分离愈靠前。边界层分离后，流动特征发生了变化。如：（1）从分离点不断脱离出旋涡，在分离点下游形成不稳定的旋涡区，从而使得主流区由原来的无涡区变成有涡。（2）物面上压力分布由原来的几乎对称分布变成不对称分布，在分离点后出现低压区（或负压区），从而大大增加了绕流物体的阻力。Folie845.55.5、边界层的分离现象、边界层的分离现象Folie855.55.5、边界层的分离现象、边界层的分离现象
