非牛顿幂律钻井液在螺旋管中流动与传热的实验与数值研究

时间:2026年2月2日
来源:International Journal of Thermal Sciences

编辑推荐:

功率流变钻井液螺旋管内热力耦合特性研究通过实验与数值模拟结合的方法,系统分析了流体密度、流变参数(K值和n值)及螺旋管曲率比对传热与流动阻力的影响规律。研究发现:降低流体密度和提升入口温度可强化对流换热;减小K值和n值协同增强传热同时降低流动阻力;增大曲率比可显著提升热力性能并放大流变参数影响。基于实验与数值数据,创新性地建立了包含Dean数、曲率比和流变指数的摩擦因子关联式,并构建了努塞尔数预测模型,为幂律流体螺旋管换热器工程优化提供理论支撑。

广告
   X   

余文毅|周敬志|何宇哲|张凯|李勋峰
中国科学院工程热物理研究所,北京,100190,中国

摘要

通过结合实验和数值方法,研究了非牛顿幂律钻井流体在螺旋管中的流动和传热特性。系统分析了流体密度、流变参数(粘度指数 K 和幂律指数 n)以及螺旋管曲率比的影响。结果表明,较低的流体密度和较高的入口温度都能提高传热性能。降低 Kn 值可以协同增强对流传热,同时降低流动阻力。此外,增加曲率比不仅改善了热流性能,还放大了 Kn 的效应。基于实验和数值数据,建立了一个考虑了迪恩数(Dean number)、曲率比和幂律指数影响的摩擦系数关联式,并提出了相应的努塞尔数(Nusselt number)预测模型。这些发现为幂律流体在螺旋管换热器中的工程应用和优化设计提供了有价值的理论指导。

引言

由于螺旋管的紧凑几何形状以及由二次流引起的增强传热效果,它们在化学工程、石油和能源行业中得到广泛应用 [1,2]。许多工业流体,包括钻井流体和聚合物溶液,表现出剪切稀化现象和幂律行为,这显著影响了流动阻力和传热。
分析螺旋管中非牛顿流动和传热的一个主要挑战是缺乏统一的分析框架。Metzner 和 Reed [3] 提出的广义雷诺数(generalized Reynolds number)为复杂流体的工程计算和理论分析提供了基础。早期的实验研究主要集中在开发经验关联式上。例如,Rajasekharan 等人 [4] 研究了不同曲率比(0.021–0.097)、雷诺数(3–25000)和幂律指数(n = 0.47–2.16)下幂律流体在螺旋管中的传热和压降。他们基于修正的格拉茨数(modified Graetz number)提出了努塞尔数关联式,但没有提供通用的摩擦系数关联式。Mishra 和 Gupta [5] 后来使用羧甲基纤维素水溶液研究了动量传递。通过引入基于壁面剪切应力的表观粘度,他们定义了雷诺数和迪恩数,并为非牛顿流体开发了通用摩擦系数关联式。尽管这些研究建立了有用的关联式,但仍需要对其机制进行理论理解。Mashelkar 和 Devarajan [6,7] 应用小曲率边界层近似来阐明二次流的形成和剪切稀化效应,并为工程设计开发了摩擦系数关联式。他们进一步将工作扩展到粘弹性流体。然而,小曲率假设限制了其在工业流动中的适用性。
尽管有这些基础研究,但完全理解螺旋管中的非牛顿流动和传热仍面临重大挑战。特别是 Zhou 和 Shah [8] 指出,现有的摩擦系数关联式往往缺乏普遍性和准确性,且大多数研究主要关注摩擦而不是传热。为了解决这些问题,Pimenta 和 Campos [9,10] 研究了在恒定壁温下垂直螺旋管中的层流摩擦和传热。他们确定迪恩数为 80 是流动行为的关键临界值。对于迪恩数高于 80 的情况,剪切稀化行为会降低摩擦系数并通过迪恩涡流增强传热,而弹性效应会增加摩擦并阻碍传热。他们提出了一个考虑粘弹性的统一关联式。Pawar 和 Sunnapwar [11,12] 通过实验和模拟研究了三种曲率比下螺旋管中牛顿流体和非牛顿流体的流动和传热。他们将努塞尔数与 Mujawar 和 Rao [13] 提出的无量纲数 M 相关联,该数表示惯性力与粘性力的比率。这个参数通过指数调节了曲率的影响。他们进一步为牛顿流体和非牛顿流体开发了努塞尔数和摩擦系数关联式,与以往的研究结果吻合良好。然而,由于迪恩数及其扩展具有明确的物理意义和广泛的验证,它们仍然是最常用的标准。
此外,几何配置的影响也得到了广泛研究。Kurnia 等人 [14] 使用数值模拟比较了平面、圆锥形和环形螺旋管(所有管截面均为正方形)中幂律流体的层流传热。他们引入了一个性能因子来评估相对于泵送功率的传热效率。幂律流体相对于水提高了传热效率,但产生了更高的压降,其中环形螺旋管的性能最佳。Valizadeh 等人 [15] 进行了参数化数值研究,比较了牛顿模型、幂律模型和 Carreau 模型。他们发现,在低速流动下,幂律模型和 Carreau 模型的预测结果相似。螺旋结构引起的二次流增加了外侧的壁面剪切应力,显示出几何形状和流变性质之间的强耦合。同时,Kushwaha 等人 [16] 研究了螺旋线圈壳结构,发现安装挡板可以增强传热但增加流动阻力,其效果随普朗特数(Prandtl number)显著变化。
尽管有这些进展,系统研究曲率、流变性质和流动条件对传热的耦合效应仍然不足。现有关联式在实际操作条件下往往存在较大偏差,无法准确捕捉由剪切稀化行为和曲率变化引起的动量和能量传输的耦合。需要综合实验和数值研究来阐明潜在机制和参数耦合效应,并为该领域相关的典型流变参数范围内流动和传热开发更可靠的预测关联式。
本研究旨在阐明螺旋管中幂律钻井流体的流动和传热机制。建立了一个实验平台,并开发了相应的数值模型,系统研究了流体密度、流变参数和曲率比对热流性能的影响。评估了常用摩擦系数和努塞尔数关联式的适用性。基于这些分析,提出了新的摩擦系数和努塞尔数预测关联式,以更准确地描述螺旋管中剪切稀化流体的热流行为。

实验装置

如图 1 所示,实验装置由四个主要回路组成:钻井流体回路、海水回路、电加热油回路和制冷回路。钻井流体由多级离心泵(容量为 0–5 m3/h)循环,其流量通过涡轮流量计和节流阀调节。由于钻井流体的热敏感性,直接加热可能导致壁温过高和流体性能下降

数值建模

在 SolidWorks 中创建了一个基于实验装置的三维模型,并在 ANSYS Fluent 中进行了模拟,该模型表示高温钻井流体通过螺旋管的向上流动,并伴有逆流的海水冷却。忽略了向环境的散热,这与测试段的隔热外墙一致。
钻井流体的热物理性质(见表 4),包括温度依赖的热导率和比热容

实验分析的结果与讨论

分别研究了三个关键参数——入口温度 Tin、广义雷诺数 Regen 和流体密度 ρ——在五个不同水平上的影响。由于高密度条件下的流动限制,某些高速度情况下的实验数据未能获得,因此共有 115 个有效数据集,如图 10 和图 11 所示。广义迪恩数 Degen 量化了离心效应和螺旋管中的二次流强度

结论

本研究采用实验和数值相结合的方法,系统研究了螺旋管中幂律钻井流体的流动和传热行为,重点关注流变参数和几何特性对 Nufc 的影响。主要结论如下:
  • (1)
    Tinρ 虽然在统计上不显著,但显示出一致的趋势:较低的 ρ 促进对流传热同时降低流动阻力,较高的 Tin 提高了整体
  • CRediT 作者贡献声明

    余文毅:撰写初稿、软件开发、方法论设计、概念构建。周敬志:监督工作、资源调配、方法论指导、实验研究。何宇哲:数据验证、实验研究、数据整理。张凯:数据验证、监督工作、实验研究。李勋峰:撰写审查与编辑、项目管理、方法论指导。

    利益冲突声明

    作者声明他们没有已知的财务利益或个人关系可能影响本文报告的工作。

    生物通微信公众号
    微信
    新浪微博


    生物通 版权所有