工程流体力学第7章.pptVIP

§7.4 平板边界层计算 设有一极薄的静止光滑平板顺流放置在二维恒定均速流场中，如图。 平板首端设为坐标原点，取平板表面为x轴，来流速度为U0且平行于x轴。因为 平板极薄，可认为对流场没有影响，因此边界层外边界上的速度处处相等，且 等于来流速度U0。根据伯努利方程，由于流速不变，边界层外边界上的压强处处 相等，即 。对于不可压缩均质流体而言，密度ρ是常数，可以提到积 分符号外，（7.7）式可写为 (7.8) 一 平板层流边界层 首先要补充两个关系式。 第一个关系式为边界层内流速分布关系式ux= ux(y)。这里 假定层流边界内流速分布和管流中的层流速度分布相同， 即 将其应用于层流边界层，管流中的r0对应于边界层的厚度 δ，r对应于(δ-y) ，umax对应于U0，u对应于ux。这样， 上式可写为 或 (7.9) 第二个补充关系式为切应力与边界层厚度的关系式 τ0=τ0(δ)。因为是层流，符合牛顿内摩擦定律， 求平板上的切应力，令y=0，得 整理简化得到 将(7.9)、(7.10)式代入(7.8)式，得 (7.10) 简化为 积分得 积分常数C由边界条件确定，当x=0时，δ=0，得到C=0。代 入上式得 化简后得到 式（7.11)即为平板层流边界层厚度沿x轴方向的变化关系。 (7.11) 将(7.11)式代入（7.10）式，化简后可得 上式即为平板层流边界层的切应力沿x轴方向的变化关系。 作用在平板上一面的摩擦阻力Df为 式中，b为平板的宽度，L为平板的长度。将(7.12)式代入上 式，积分后得到 如求平板两面的总摩擦阻力时，将上式乘以2即可。 (7.12) (7.13) 通常将绕流摩擦阻力写成如下形式 式中，Cf为摩阻系数；ρ为流体密度；U0为流体的来流速 度；A为切应力作用的面积，这里指平板面积。 由(7.13)式和(7.14)式比较可得 式中， ，ReL是为以板长L为特征长度的雷诺数。 (7.15) (7.14) 二 平板湍流边界层 对于湍流边界层的计算，同样要补充两个关系式。 这里假定从平板首端开始就是湍流边界层，并且不考虑平板 壁面粗糙度的影响。 借用圆管湍流水力光滑区的流速分布公式 将其应用于平板湍流边界层，管流中的r0对应于边界层的厚 度δ，umax对应于U0，u对应于ux。这样，上式可写 (7.16) 第二个补充关系式为切应力与边界层厚度的关系式同样借用 管流的关系式 将(7.16)、(7.17)式代入(7.8)式，得 积分得到 积分常数C由边界条件确定，当x=0时，δ=0，得到C=0。 化简后得到 (7.17) (7.18) 上式即为平板湍流边界层厚度沿x轴方向的变化关系。 将式(7.18)代入（7.17）式，可得 上式即为平板湍流边界层的切应力沿x轴方向的变化 关系。 作用在平板上一面的摩擦阻力Df为 将(7.19)式代入上式，积分后得到 (7.19) 如果用绕流摩擦阻力的通用形式（7.14）表示，摩 阻系数为 实验表明，上式中的0.072改为0.074，则与实验的 结果符合得更好。 (7.20) 可以看出：与层流边界层比较，当ReL增加时，湍流的Cf要比层流的Cf减小得慢些；在同一ReL的情况下，湍流的Cf要比层流的Cf大得多，即在湍流边界层内，产生更大的摩擦阻力。 三 平板混合边界层 前面讨论的是假定整个平板上的边界层都处于湍流状态，但 实际上，当雷诺数增加到某一数值后，而且平板长度 时，平板的前部是层流边界层，后部是湍流边界层，在层流 和湍流边界层之间还有过渡段。这种边界层称为混合边界 层。 由于混合边界层内流动情况很复杂，在进行计算时，作了两 个假设： 1 层流边界层转变为湍流边界层是在xcr处突然发生的，没有过渡段； 2 混合边界层的湍流边界层可以看作是从平板首端开始的湍流边界层的一部分。 根据假设，整个平板混合边界层的摩擦阻力，由转 捩点xcr前层流边界层的摩擦阻力和转捩点xcr后湍流 边界层的摩擦阻力两部分组成，即 式中 Cf,m、Cf,t、Cf,l分别为混合边界层、湍流边 界层、层流边界层的摩阻系数，这里用Cf,