幂律流体在环形通道中的流动规律.doc

幂律流体在环形通道中的流动规律 0 前言 在许多工程领域中经常会遇到非牛顿流体在环空中流动的情况，例如在石油工程中泥浆或钻井液在钻杆和套管间的流动，类似的例子在化学工程、生物食品工业和摩擦润滑中都会经常遇到。按照非牛顿流体的分类，许多情况下都可将其看成是幂律流体。幂律流体在这样的环空中的流动规律直接关系到具体工艺过程的效率、成本和质量。因此研究幂律流体在环空中的流动规律有着非常重要的工程实际意义。 1 运动方程及求解 假设不可压缩的幂律流体在如图1所示的同心环空中作轴向稳定等温的层流流动，Ri为环形空间内径，Ro为环形空间外径，Rλ为环形空间内最大速度所对应的半径。 图1 环空的几何结构 这样幂律流体在环形空间的速度为： (1) 同时其偏应力张量为: (2) 式中为剪切速率。 这样运动方程可以简化为： （3） 引入有效压力： （4） （3）式可以简化为： （5） 定解条件为： （6） （7） 将（5）式对r积分，得到： （8） 根据（7）式，在处，剪切速率，剪切应力也应为零，故由（8）式解得： （9） 将（9）式代到（8）式有： （10） （1）当时，，，幂律流体的本构方程为： （11） 由（10）、（11）式可得： （12） 将上式从到r积分并利用定解条件（6），可得时的速度分布： （13） （2）当时，，，幂律流体的本构方程为： （14） 由（10）、（11）式可得： （15） 将上式从r到积分并利用定解条件（6），可得时的速度分布： （16） （13）式和（16）式即为幂律流体在环空中的速度分布。 从（13）式和（16）式可以看出，为计算环空内的速度分布，必须求出。显然在处，由（13）式和（16）式求得的速度应该相等，即： （17） 解出后即可由（13）式和（16）式求出速度分布。根据速度分布，可由下式求出幂律流体在环空内的体积流量Q： （18） 将（13）式和（16）式代入（18）式，通过分部积分并利用（17）式可得： （19） 将上式进行分部积分，可得： （19a） 当n=1时为牛顿流体（），根据（13）、（16）、（17）式和（19）式可求得牛顿流体在环空内流动时的速度分布、和体积流量Q的表达式： （） （20） （） （21） （22） （23） 3 算例分析 求解的方程（17）是积分限和积分表达式中均含有未知数的复杂的积分方程，这里采用迭代和数值积分相结合的方法进行求解。解出后利用数值积分法求解方程（13）和（16）及方程（19），得到环空内的速度分布及相应的流量，同时分析得到幂律流体在环空内流动时的流动规律。 3.1 各因素对环空速度分布 （1）流变指数n对速度分布和流量的影响 流变指数n对环空内速度分布有显著影响，当n小于1时，幂律流体的速度分布位于牛顿流体的上方，且n越小速度分布曲线越粗长；当n大于1时，速度分布曲线位于牛顿流体上方，且n越大抛物线分布曲线越平缓。这是因为流变指数n反映了幂律流体的非牛顿程度，n越小幂律流体的剪切减稠性能越显著，因而在相同的压降下速度分布越粗长；n越大非牛顿性越强，因而流体的