黏弹性流体在内管做轴向运动的偏心环空中的速度分布.pdfVIP

维普资讯 大 庆 石 油 学 院 学 报 第 29卷 第 l期 2005年 2月 JOURNAlOFDAQING PETR()IEUM INSTITUTE Vo1．29 No．1 Feb ． 2005 黏弹性流体在 内管做轴向运动的 偏心环空中的速度分布 杨树人 ，王春生，杨 英 (大庆石油学院 石油工程学院，黑龙江 大庆 163318) 摘 要：利用下随体 Maxwell本构模型，建立了双极坐标系下黏弹性流体在内管做轴 向运动的偏心环空中定常流的 控制方程 ，利用控制体积法对其离散化处理，并用ADI方法求解流场的速度分布。分析了偏心距等参数对速度分布和流 量的影响．结果表明．偏心距对流场的速度分布和流量影响较大． 关 键 词：黏弹性流体；本构方程；偏心环空；速度分布 中图分类号：TE933．2 文献标识码 ：A 文章编号：1000—1891(2005)01—011O—O2 l 控制方程 假设黏弹性流体为各向同性不可压流体，聚驱井井筒由2个轴线相互平行且距离为 的偏心圆管构 成，半径为尺。的外管静止，半径为尺，的内管 以匀速 7．22。做轴向运动，且又有一个平行于环形空间轴线 的 压力梯度 P作用于流体上 ，黏弹性流体在偏心环空中做定常等温层流流动．在此假设条件下，柱坐标系下 流动的运动方程和下随体 Maxwell本构方程 为 ÷，-C，，-+÷，-C，+户一0， ㈩ S．k+ 0Sik— q0Aik， (2) 式 (1～2)中：r为径向坐标 ；S，为应力张量的协变分量；0为周向向坐标 ；S，为S，下随体导数；。为松弛时 间； 为流体黏度；A，．为一阶Rivlin—Erickson变形张量的协变分量．由式 (1)和式(2)可得到控制方程 攀a-，+’攀a6I+’一… ㈦、 式中：u，为速度． 利用双极 坐标变换 ，见图 1．由于物理平 ／ 一 面左右对称，见图 1(a)。只需将右半圆经坐标 变换成计算平面即可，见图2(b)．柱坐标和 双极坐标变换关系式为_2一 、r一 (r—R，)／(r。一R )； (4a) -y一 0 ， (4b) ＼ 式中：rc1为流场外边界处的 r坐标．由式 (3) 和式(4)可得双极坐标系下的控制方程为 (a)物理平面 (b)计算平面 图 l 坐标变换 收稿 日期：2004—10—12；审稿人：韩洪升；编辑：王文礼 基金项 目：国家 自然科学基金资助项 目50374018) 作者简介：杨树人(1963一)，男，博士生，副教授 ，主要从事非牛顿流体力学方面的研究 · 110 · 维普资讯 第 1期 杨树人等：黏弹性流体在内管做轴向运动的偏心环空中的速度分布 { }+ {丽 )+[Rl一 。_R_)p