工程流体力学第8章边界层基础及绕流运动

本章主要内容本章主要内容§§8.18.1边界层的基本概念§§8.28.2边界层微分方程§§8.38.3边界层的动量积分方程§§8.48.4平板边界层的近似计算§§8.58.5边界层的分离§§8.68.6绕流阻力中国海洋大学海洋工程系王树青高等流体力学§§8.18.1边界层的基本概念III外部势流II尾部流区域I边界层边界层外边界边界层外边界中国海洋大学海洋工程系王树青高等流体力学一、边界层的提出一、边界层的提出™绕流运动‰绕流运动为高雷诺数运动一般物体的特征长度在l=0.01－10m范围，当物体在空气或水中以速度U=0.1－100m/s运动时，相应的雷诺数约在100－109之间。普通汽车和船舶以正常速度行驶时，空气和水的雷诺数均在106以上。飞机绕流的雷诺数则更高，因此大Re数流动是普遍存在的现象。中国海洋大学海洋工程系王树青高等流体力学一、边界层的提出一、边界层的提出™绕流运动‰高雷诺数运动：忽略粘性，简化为理想流体运动‰达朗贝尔佯谬：1752年在《试论流体阻力的新理论》中考虑没有粘性的不可压缩流体，结果得到运动物体受到的阻力为零的结论。他本人不满意这个结论，但又得不到正确的解释，成为一个所谓达朗贝尔佯谬。法国数学家、力学家、哲学家。中国海洋大学海洋工程系王树青高等流体力学一、边界层的提出一、边界层的提出™边界层的提出‰普朗特1904年，在德国举行的第三届国际数学家学会上，德国著名的力学家普朗特第一次提出了边界层的概念（论粘性很小的流体的运动）。他认为对于水和空气等黏度很小的流体，在大雷诺数下绕物体流动时，黏性对流动的影响仅限于紧贴物体壁面的薄层中，而在这一薄层外黏性影响很小，完全可以忽略不计，这一薄层称为边界层。普朗特的这一理论，在流体力学的发展史上有划时代的意义。中国海洋大学海洋工程系王树青高等流体力学一、边界层的提出一、边界层的提出™边界层的主要内容（1）固壁附近边界层内的流动，粘性力和惯性力同量级，必须考虑粘性的影响，为有旋运动；（2）边界层以外的流动区域,该区域内流体速度变化很小，可近似看成是理想流体.粘性剪切流无粘性流场中国海洋大学海洋工程系王树青高等流体力学一、边界层的提出一、边界层的提出™优点引入边界层后，可将流场的求解可分为两个区进行9边界层内流动必须计入流体的粘性影响可利用动量方程求得近似解。9边界层外流动视为理想流体流动，可按势流求解。它的提出为解决粘性流体绕流问题开辟了新途径，并使流体绕流运动中一些复杂现象得到解释中国海洋大学海洋工程系王树青高等流体力学二、边界层的形成二、边界层的形成u∞u∞u∞u∞uuux=0.99u∞边界层区主体区或外流区U层流边界层过渡区湍流边界层xy中国海洋大学海洋工程系王树青高等流体力学二、边界层的形成二、边界层的形成由于具有粘性，紧贴壁面的流体必然附着于物体表面上，其速度为零。近壁面的流体相继受阻而减速。u∞u∞u∞u∞uuux=0.99u∞边界层区主体区或外流区‰一方面，随着流体向前流动，速度受到影响的区域逐渐扩大‰另一方面，随着与板面法向距离的增大，板面对流体的减速作用逐渐减弱。中国海洋大学海洋工程系王树青高等流体力学二、边界层的形成二、边界层的形成u∞u∞u∞u∞uuux=0.99u∞边界层区主体区或外流区™在离板面一定距离之外的流体速度就基本上未受板面影响接近了的主流速度。减速作用发生在紧邻板面的很薄的流体层中，这一薄层称之为边界层;‰边界层内：沿板面法向的速度梯度很大，剪应力不可忽略。——粘性流体的流动‰边界层外：不存在速度梯度或速度梯度很小，剪应力可以忽略。——理想流体运动中国海洋大学海洋工程系王树青高等流体力学中国海洋大学海洋工程系王树青高等流体力学三、边界层的主要特征三、边界层的主要特征(1)与物体的特征长度相比，边界层的厚度很小δL。(2)边界层内沿厚度方向，速度梯度很大，为有旋运动。(3)边界层厚度沿流体流动方向是增加的。(4)由于边界层很薄，可以近似认为边界层中各截面上的压强等于同一截面上边界层外边界上的压强值。u∞u∞u∞u∞uuux=0.99u∞边界层区主体区或外流区中国海洋大学海洋工程系王树青高等流体力学层流边界层过渡区域紊流边界层层流底层三、边界层的主要特征三、边界层的主要特征(5)在边界层内，粘性力与惯性力同一数量级。(6)边界层内的流态，也有层流和紊流两种流态。中国海洋大学海洋工程系王树青高等流体力学三、边界层的主要特征三、边界层的主要特征层流边界层过渡区域紊流边界层层流底层™层流边界层向紊流边界层的转变ν=xURex0cxx=ν=ccxURe0cxx=5105×=中国海洋大学海洋工程系王树青高等流体力学四、边界层的厚度四、边界层的厚度1．名义厚度定义为速度达到外流速度99%时离壁面的垂直距离，称为名义厚度δ(x)。u∞u∞u∞u∞uuux=0.99u∞边界层区主体区或外流区中国海洋大学海洋工程系王树青高等流体力学四、边界层的厚度四、边界层的厚度2．位移厚度又称边界层流量排挤厚度，是指由于边界层的存在，使得外部流动按理想流体处理时，其流动的虚拟边界向壁面以外移动的距离。xy0u099.0uxuδ1δ∫∫δδρ−ρ=δρ00010dyudyuux∫δ−=δ001)1(dyuux中国海洋大学海洋工程系王树青高等流体力学四、边界层的厚度四、边界层的厚度2．位移厚度†表明：由于流体的粘性作用，存在着流动被阻滞了的边界层，为了满足连续性方程，流道就得扩张，才能让一定量的流体通过，因此流线向外偏斜，被排移了δ1的距离；也就是说，由于边界层的存在排移了厚度为δ1的非粘性流体的流量。流线Yy=Y+δ1U∞δ1如图，兰线为一条流线，由于边界层的存在使它向上偏移了排量厚度δ1的距离中国海洋大学海洋工程系王树青高等流体力学四、边界层的厚度四、边界层的厚度3．动量损失厚度是指由于边界层的存在，边界层内所损失的动量折合成按理想流体处理时，具有相同动量的等效厚度。xy0u099.0uxuδ1δ∫∫δδρ−ρ=δρ0200220dyudyuuuxx∫δ−=δ0002)1(dyuuuuxx中国海洋大学海洋工程系王树青高等流体力学四、边界层的厚度四、边界层的厚度4．能量损失厚度是指由于边界层的存在，边界层内所损失的能量折合成按理想流体处理时，具有相同能量的等效厚度。xy0u099.0uxuδ1δ∫∫δδρ−ρ=δρ020202030dyuudyuuuuxxx∫δ−=δ020203)1(dyuuuuxx中国海洋大学海洋工程系王树青高等流体力学§§8.28.2边界层微分方程™将利用边界层流动的特点（边界层的厚度与物体的特征长度相比为一小量）对N-S方程进行简化从而导出层流边界层微分方程。™在简化过程中，假定流动为二维恒定不可压定常流，不考虑质量力。⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫=∂∂+∂∂∂∂+∂∂+∂∂ρ−=∂∂+∂∂∂∂+∂∂+∂∂ρ−=∂∂+∂∂0)(1)(122222222yuxuyuxuvypyuuxuuyuxuvxpyuuxuuyxyyyyyxxxxyxx⎪⎪⎪⎪⎭⎪⎪⎪⎪⎬⎫=∂∂+∂∂∂∂+∂∂+∂∂ρ−=∂∂+∂∂∂∂+∂∂+∂∂ρ−=∂∂+∂∂0)(1)(122222222yuxuyuxuvypyuuxuuyuxuvxpyuuxuuyxyyyyyxxxxyxxz将上述方程组无量纲化。为此考虑如图所示的半无穷绕流平板，假定无穷远来流速度U0，流动绕过平板时在平板附近形成边界层，其厚度为δ，平板前缘至某点的距离为L。取δ和L为特征量，可定义如下的无量纲量：Lxx=0Lyy=000Uuuxx=00Uuuyy=200Uppρ=⎪⎪⎪⎪⎪⎪⎪⎭⎪⎪⎪⎪⎪⎪⎪⎬⎫=∂∂+∂∂δδδ⋅δδ⋅∂∂+∂∂+∂∂−=∂∂+∂∂δδδ⋅δ⋅∂∂+∂∂+∂∂−=∂∂+∂∂110111)(Re111)(111)(Re1000000200200220020/000000020200200220020/000000yuxuyuxuypyuuxuuyuxuxpyuuxuuyxyyyyxxxxxyxx⎪⎪⎭⎪⎪⎬⎫=∂∂+∂∂∂∂+∂∂−=∂∂+∂∂0Re10000200200000000yuxuyuxpyuuxuuyxxxyxxLxx=0Lyy=000Uuuxx=00Uuuyy=中国海洋大学海洋工程系王树青高等流体力学普朗特边界层方程⎪⎪⎭⎪⎪⎬⎫=∂∂+∂∂∂∂ρµ+∂∂ρ−=∂∂+∂∂0122yuxuyuxpyuuxuuyxxxyxx边界条件：（1）y＝0：0,0==yxuu（2）y＝∞：0uux=两个方程三个未知数：0=∂∂yp中国海洋大学海洋工程系王树青高等流体力学紊流边界层方程⎪⎪⎭⎪⎪⎬⎫=∂∂+∂∂′′−∂∂∂∂+∂∂−=∂∂+∂∂0)(11yuxuuuyuyxpyuuxuuyxyxxxyxxρµρρ™说明：（1）边界层方程虽然比N-S方程简化了，但仍是非线性：（2）布拉休斯应用该方程求解平板的层流边界层的解；（3）近似解法－动量积分方程；中国海洋大学海洋工程系王树青高等流体力学§§8.38.3边界层动量积分方程ƒ边界层的动量积分方程是对边界层内流动的再简化。ƒ其推导过程有两种方法：一种是沿边界层厚度方向积分边界层的方程组，一种是在边界层内直接应用动量守恒原理。ƒ下面的推导采用第二种方法。中国海洋大学海洋工程系王树青高等流体力学xyO0uABCDdxxdsθδδ+δd™假设：（1）边界层很薄，忽略质量力：（2）流动为恒定平面流动；（3）dx很小，BD和AC近似为直线；™动量方程：∑=−−xACABCDFKKK中国海洋大学海洋工程系王树青高等流体力学xyO0uABCDdxxdsθδδ+δd™动量变化：∑=−−xACABCDFKKK‰AB面：∫δρ=ρ0dyuqxAB∫δρ=02dyuKxAB‰CD面：dxdyuxdyudxxKKKxxABABCD)(0202∫∫δδρ∂∂+ρ=∂∂+=‰AC面：dxdyuxqqqxABCDAC)(0∫δρ∂∂=ρ−ρ=ρdxdyuxuuqKxACAC)(000∫δρ∂∂=ρ=dxdyuxudxdyuxKxx)()(0002∫∫δδρ∂∂−ρ∂∂=∆中国海洋大学海洋工程系王树青高等流体力学™外力分析：∑=−−xACABCDFKKKxyO0uABCDdxxdsθδδ+δd‰压强沿y向均匀分布：0=∂∂yp‰AB面压强：ppdxxpp∂∂+‰CD面压强：dxxpp∂∂+‰AC面压强：dxxpp∂∂+21dxxpp∂∂+21‰BD面摩擦：0τ0τ=∑xFδp))((δ+δ∂∂+−ddxxppθ∂∂++sin)21(dsdxxppdx0τ−中国海洋大学海洋工程系王树青高等流体力学™动量积分方程dxdsdxxppddxxpppFx0sin)21())((τ−θ∂∂++δ+δ∂∂+−δ=∑dxdxdxpdxxpFx021τ−δ∂∂−δ⋅∂∂−=∑dxdxdxdpFx0τ−δ⋅−=∑00002τ−δ−=ρ−ρ∫∫δδdxdpdyudxdudyudxdxx卡门动量积分方程dxdyuxudxdyuxKxx)()(0002∫∫δδρ∂∂−ρ∂∂=∆中国海洋大学海洋工程系王树青高等流体力学™动量积分方程00002τ−δ−=ρ−ρ∫∫δδdxdpdyudxdudyudxdxx‰未知数：δ，p，u0，ux，τ0‰势流理论求：u0‰能量方程求：p‰补充方程：（1）边界层内的流速分布ux=f(y)（2）切应力τ0随边界层厚度δ的关系式τ0=g(δ)中国海洋大学海洋工程系王树青高等流体力学§§8.48.4平板边界层的近似计算一、三种计算一、三种计算xyOUxδLLδ层流边界层过渡区湍流边界层ν=LLUxReν=ccUxRe‰层流边界层：cLReRe‰混合边界层：cLReRe‰紊流边界层：cLReRe中国海洋大学海洋工程系王树青高等流体力学二、平板边界层的计算公式二、平板边界层的计算公式™恒定均匀来流的平板边界层，其外边界满足‰外边界上的流速处处相等，且等于来流速度；0,00==dxduUu‰外边界外按理想流体处理，据能量方程，由于流速不变，则外边界上的压