25黄红梅_pdc钻头旋转流场-稿件

《25黄红梅_pdc钻头旋转流场-稿件》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、旋转PDC钻头流场数值模拟计算研究黄红梅1，翟应虎2（1.中国石油勘探开发研究院机械所，北京910信箱100083；2.中国石油大学石油工程学院，北京昌平102249）摘要：利用计算流体动力学技术对旋转PDC钻头的三维湍流流场进行了数值模拟。模拟中考虑了钻头的切削齿和钻头的射流喷嘴对流场的影响。流场的网格划分采用局部加密的混合网格形式，保证了计算精度和运算速度的要求。流场的控制方程在单一旋转参考坐标系中建立，流场的湍流模型采用的是重整化群模型(renormalizationgroup,RNG,　model)。旋转流场的三维模拟结果揭示了

2、旋转流场与静止流场的区别，为进一步的钻头水力结构优化分析奠定了理论基础。关键词：PDC钻头；流场；数值模拟；旋转目前，国内采用计算流体动力学技术对PDC钻头流场的研究，由二维、单喷嘴发展到最近的三维、旋转、非对称多个喷嘴的流动数值模拟[1]，取得了很大的进步。但在钻头物理模型的选取上，没有考虑真实的钻头模型，如钻头上的切削齿和射流喷嘴等。因此，模拟结果用来指导PDC钻头的水力结构优化设计还有一定的差距。1旋转PDC钻头的流场数值模拟1.1建立钻头流动的几何模型及网格划分选用一只BK432PDC钻头进行流场数值模拟。钻头的几何模型由CAD

3、软件完成，如图1所示。该钻头有四片刮刀，2对喷嘴,喷嘴出口直径均为11mm。靠近钻头中心部位的喷嘴倾斜角度分别为12度、18度；另两只喷嘴倾斜角度分别为22度、25度。利用已建立的PDC钻头几何模型，在CAD软件中建立了钻头的流动模型。对流动模型进行网格划分。在流场的网格划分中，将几何模型分为3部分，最上面的部分为逐渐稀疏的棱柱体网格；中间部分为过渡网格，其中中心液流为棱柱体网格，外圈为四面体网格；最下面的部分包含有4片刮刀和切削齿，此部分为四面体网格，局部加密。采用这种局部加密的网格形式，既满足了计算精度的要求，又兼顾了模拟的运算速度

4、。图1　钻头的几何模型1.2旋转流场的控制方程5旋转PDC钻头的流场中，只需考虑钻头的旋转，静止的井壁不需要考虑，因此这种流动属于单一旋转参考坐标系中的流动，需要在单一旋转参考坐标系中进行模型建立。在单一旋转参考坐标系中，流体的加速度作为一个附加项出现在旋转坐标系运动方程中的动量方程里，这时将绝对速度或相对速度各自作为独立参数，两种速度之间的关系如下：　　　　　　　　　　（1）其中，指角速度向量（即旋转坐标系的角速度）是旋转坐标系中的位置向量。惯性坐标系中的动量方程表达如下：　　　　　　　　（2）式中为静压力，为应力张量　　　　　　　　

5、　　（3）旋转坐标系中用绝对速度表达的动量方程左侧如下：　　　　　　　　　（4）根据与的关系代入上式，消去得　　　　　　（5）其中为哥氏力。以上方程未考虑与时间相关的项，因此不能用动量方程的相对速度表达式准确的模拟随时间变化的角速度。旋转坐标系中的连续性方程，不论是使用相对速度，还是使用绝对速度，都可表示如下：　　　　　　　　　　　　　（6）1.3旋转流场的湍流模型如果模拟的流场是具有大量漩涡的湍流流动（气旋流动，漩涡射流），就应该使用更高级的湍流模型，如：RNG模型，可实行的模型或者雷诺应力模型。选择哪种模型依赖于漩涡的强度，它可以用

6、漩涡的数量来度量。漩涡数量定义为角动量轴向流量与轴向流量动量的比值。对于较弱的中等涡流（S<0.5），RNG模型和可实行的模型比标准5模型要更好一些。对于强度较高的漩涡流动（S>0.5），应该使用雷诺应力（RSM）模型。PDC钻头的旋转流动属于中等强度的涡流，可采用RNG模型。RNG模型是将非稳态Navier-Stokes方程对一个平衡态作Gauss统计展开，并用对脉动频谱的波数段作滤波的方法，从理论上导出高Re数模型，所得出的方程形式上同标准模型完全一样，但不同的是5个系数之值不是根据实验数据而是由理论分析得出。这套系数的最新结果为［

7、2］：；；；；（7）；；高Re数方程采用上述系数时就构成了重整化群模型(renormalizationgroup,RNG,　model)。显然，与标准模型相比，RNG模型的最大特点在于在方程产生项的系数的计算中引入了主流的时均应变率。这样在RNG模型中之值不仅与流动情况有关，而且在同一问题中也还是空间系数的函数。由于RNG模型采用高Re数方程，因而在接近壁面处要采用壁面函数法来处理。1.4流动的边界条件及计算介质流动的入口条件根据钻头的工作排量来估算。设定钻头工作排量为30升每秒，则初始速度为11.7m/s，即，,,。流动的出口条件采用

8、的是无回流的出口边界条件，采用出口截面局部单向化假设。固体壁面边界条件为壁面上满足无滑移条件，近壁区采用壁面函数法处理。钻头的旋转速度为100转每分，速度的松弛因子设为0.4。计算介质采用清水介质。2模拟结