统一的泥沙沉降阻力系数及沉速公式.pdfVIP

．296． 第六届全国泥沙基本理论研究学术讨论会论文集 统一的泥沙沉降阻力系数及沉速公式 宋志尧 王小平 孔 俊 章卫胜 (河海大学水文水资源与水利工程科学国家重点实验室 南京 210098) 摘要：本文基于球体沉降的实验数据，结合大量前人的研究成果，从函数构造论的角度获得了球体泥 沙阻力系数与其粒径雷诺数关系曲线的近似数学表达式，并得出了相应球体泥沙的沉速公式。通过与许 多已有公式的比较研究表明，除超高雷诺数(2．6×105)的情况外，本文公式统一了层流区、过渡区及紊 流区球体阻力系数以及沉速与其粒径雷诺数之间的函数关系。本文公式结构较为简洁，计算结果与已有 的统一公式较为一致。最后，考虑天然泥沙颗粒的几何特性，修正球体泥沙沉速公式中的系数，得出了 实用的泥沙沉速公式，并与张瑞瑾、Zanke、Julien和Soulsby等公式进行了比较，结果均优于其他此类 公式，可用于天然泥沙的沉速预测确定。 关键词：阻力系数沙粒雷诺数泥沙沉速 1前言 沉速是泥沙的基本特性之一，是研究泥沙运动尤其是悬移质运动时的重要参数，如 水流挟沙能力(张瑞瑾、窦国仁)、泥沙起动流速(杨志达)、河口泥沙的絮凝及再悬浮(Metha) 等的研究。同时，在实际工程中，有关泥沙问题的研究如沉沙池设计、保滩促淤工程及 航道规划、设计和整治等，不可避免地要涉及泥沙沉速的确定。 天然泥沙颗粒的形态各异，其沉降过程复杂、影响其沉速的因素较多，不易从中认 识其规律性。为此，人们将泥沙颗粒概化为球体，通过理论分析和物理实验，对球体的 沉降过程及其沉速进行了全面的研究和深入的分析，得出了许多规律性的结论。其中关 于沉速研究的主要是球体的沉降机理可视为对称绕流运动，其绕流阻力F可表达为： 9 F：Cd孚d29竺 (1) 式中：o为阻力系数；d为球体直径；，，为水的容重；CO为泥沙沉速。并通过大量实验研 究得到了阻力系数与沙粒雷诺数Re=驯u(D为水的运动黏性系数)之间的关系曲线。 另一方面，泥沙沉速是指单颗粒泥沙在无限大静止清水体中匀速下沉时的速度(王昌 杰)。因为球体颗粒在水中受到的有效重力职可表达为 作者简介：宋志尧，男，1964年出生，浙江余姚人，教授。 三、泥沙运动基本理论及泥沙输移 ·297· (2) 暇=(Y。-Y)6兀d3 其中儿为泥沙容重。故沉速可表达为 舻1／西叫彳∥舻怯。√竽∥ 由上式可知，沉速计算的关键是。的确定，而。是尺P的函数，因此白又是∞、d 的函数，故沉速表达式(3)实际上是一个关于∞的隐函数。当然，对于特定的绕流状态， 沉速表达式可简化为实用的计算公式，如层流时的Stokes公式和紊流时的沉速公式，而 对介于层流与紊流状态之间的绕流过渡区，由于层流时的黏性阻力和紊流时的形状阻力 同时存在且关系复杂，尽管许多学者对此做了大量深入的研究，如岗恰洛夫、沙玉清、 张瑞瑾、窦国仁Zanke、Julien和Soulsby等，但目前尚无一致的泥沙沉速公式。 本文根据球体白一RP关系实验曲线，并借鉴前人的研究成果，通过函数构造获得 了。关于RP的函数表达式，统一了层流区、过渡区和紊流区的球体沉降阻力系数公式， 从而得到了相应一致的沉速公式，并讨论了天然泥沙的沉速计算公式。 2沉速阻力系数 2．1 研究回顾 早在1851年，Stokes在水流惯性力远小于黏性力的条件下，将Ⅳ一S方程线性化， 得到了a反比于RP的重要结论，即 Cd=詈 (4) 后来的研究表明，Stokes得出的式(4)只适用于Re0．4的情况。为此，不少学者做 了进一步的研究。1927年，Oseen在Stokes分析的基础上，在运动方程中保留了一部分 惯性力作用项，