偏心圆环管中非牛顿流体流动数值计算.pdfVIP

偏心圆环管中非牛顿流体流动的数值计算 许彬，张敏，房萍萍 南京理工大学动力工程学院，南京（210094 ） 摘 要：通过坐标变换，在偏心圆环域中生成正交结构化网格，用有限差分法以及二阶精度 的离散格式，研究计算非牛顿流体在充分发展层流情况下的流动特性，并将其计算结果同牛 顿流体和实验结果相比较。在不同偏心率(0 ≤ ε * ≤ 0.6) ，不同半径率(0.2 ≤ r* ≤ 0.8)和不同流 体指数(0.2 ≤ n ≤ 1)范围内寻找流动滞止规律，为石油化工工业生产过程提供有参考性的结 果。 关键词：偏心圆环管，层流，幂律流体 1. 引言 在石油、化学和食品工业等领域，诸如压缩机、传输和打井过程中，经常会遇到非牛顿 流体流动的问题，因此研究非牛顿流体在圆环管中的流动特性对工业生产过程的安全性和效 益都十分必要。 非牛顿流体的剪切应变力与剪切应变张量成非线性的幂函数关系。对于一切流体可用 Ostwald-de Waele 关系式来描述这种关系， ⎡ 1 n−1 ⎤ −⎢ ∆ ∆ ⎥ τ i ,j K ( i ,j : i ,j ) ∆ i ,j µ a ∆ ij (1) ⎢ 2 ⎥ ⎣ ⎦ [1] 其中，K 为物性系数，n 为流动特性系数，∆ i j 是应变张量，µ a 为平均粘性系数 。对 于剪切稀化流体或者伪塑性流体来说 n 小于 1，对于剪切稠化流体或者膨胀流体 n 大于 1， 对于牛顿流体 n 等于 1。关于这方面的描述可参考 Skelland[2], Bird et al.[1], Fenner [3] ，Cho 和 Hartnett [4]等其他研究者的文章。 对于偏心圆环管中牛顿流体流动情况的研究，Piercy et al.[5]是最早的研究者之一，他 给出了充分发展层流速度场和摩擦因子的分析解。Snyder 和 Goldstein[6]等人通过使用双 极坐标系也进行了类似的研究。以上研究者的工作主要局限于牛顿流体，而对非牛顿流体的 研究和分析还是很少，而且很不完善。本文采用有限差分法，对非牛顿流体在具有较大偏心 率(0 ≤ ε * ≤ 0.6)和半径率(0.2 ≤ r* ≤ 0.8)范围内的流动进行数值计算和研究分析，并得到有 参考价值的结果。 2. 基本数学表达式 轴向垂直偏心圆环管的几何形状和坐标系如图 1 所示。外圆环半径为 r ，内圆环半径 o 为 r ，坐标原点设在外圆环的圆心，偏心距 ε 定义为两个圆心之间的距离，内圆环中 i 心点到外圆环边界的距离用形状函数R (ψ ) 来表示。在本文中，考虑的是常物性充分发展剪 切稀化层流在偏心圆环管中的流动情况。通过下面的坐标变换，可以将图 1 的偏心圆环管转 换成无量纲的单位圆域， r − r r i Ψ Ψ , (2) B (Ψ)r0 其中，无量纲形状函数B(ψ ) 为， -1- [ ] B (Ψ) (R (Ψ) − r ) / r 0