非牛顿流体cmc水溶液三维流动的k模型计算.docxVIP

非牛顿流体cmc水溶液三维流动的k模型计算 非牛顿流体的湍流流动模型 通过比较搅拌槽中非牛顿流的数值，模拟是计算流量的一个新领域，它是由流量性质的差异和水流本身的复杂性决定的。同时,湍流模型中模型常数的测定一般都是针对牛顿流体来测定的,很少有针对非牛顿流体的。比如对k-ε模型中的系数C2,它是通过测定格栅后湍流的衰减得到的。尽管有利用假塑性流体的实验研究,但是其并未给出C2与流动指数和湍流雷诺数的明确关系。目前尚没有适用于非牛顿流体的专用湍流模型,同时也难于确定现有湍流模型是否适用于非牛顿流体的湍流计算,因而对非牛顿流体湍流流动数值模拟的研究非常少。Takami用雷诺应力模型模拟了弯曲管中假塑性流体的湍流流动,流动指数n从1到0.9时,计算与实验值吻合较好,但在n=0.76时,时均速度的误差达到15%,湍流强度的误差达到30%。Venneker利用k-ε模型模拟了假塑性流体CMC水溶液(n=0.77)在搅拌槽内的速度场与湍流参数,计算采用“黑箱”模型,指定桨叶叶端处的边界条件。计算结果与实验结果定性符合,计算得到的轴向及径向速度均较实验结果偏低。 本文尝试利用k-ε模型计算假塑性流体CMC水溶液在搅拌槽内的三维流动场,并用PIV的测量结果进行验证,以检验计算结果的准确性,同时分析非牛顿流体湍流流动的一些特点。 1 k-模型中粘度的控制方程 对变量?的守恒方程可以用下式表示 式中(Ⅰ)为时间项,(Ⅱ)为对流项,(Ⅲ)为扩散项,(Ⅳ)为源项。上式对应于质量、动量、湍流动能和耗散率的守恒方程时,?分别为1、uj、k和ε。湍流模型采用标准的k-ε模型。圆柱坐标系下,控制方程中各项的含义参看文献。 在k-ε模型中有效粘度μe是分子粘度与湍流粘度的和。湍流粘度μT=Cμρk2/ε。分子粘度采用非牛顿流体的表观粘度。CMC水溶液属于假塑性流体,其表观粘度可以由下式表示 ηa=Kγ˙n?1=Kηa=Κγ˙n-1=Κ12tr(γ˙)2??????√12tr(γ˙)2n-1 12tr(γ˙)2=212tr(γ˙)2=2?v?r?v?r2+vr+1r?w?θvr+1r?w?θ2+?u?z?u?z2+r??rr??rwrwr+1r?v?θ+1r?v?θ2+1r?u?θ+w?z1r?u?θ+w?z2+?v?z+?u?r]2?v?z+?u?r]2 上式ηa是假塑性流体的表观粘度;12tr(γ˙)212tr(γ˙)2是应变速率张量的第二不变量。 2 计算策略 2.1 搅拌雷诺数的确定 计算所采用的搅拌槽槽体为圆柱形,均布四块挡板,搅拌桨为四斜叶涡轮。搅拌槽直径T=0.5 m,液位高H=0.4 m,挡板宽为T/10。搅拌桨直径D/T=0.6,桨叶离底距离C/T=0.2。搅拌采用下压式操作。工作介质为质量分数0.6%的CMC水溶液,其流变特性可用幂律模型来表示。利用RV-Ⅱ型同心圆筒粘度计测定了其流变特性,得到稠度系数K=0.156 kg/(s1-n·m),流变指数n=0.76。计算中搅拌转速为80 r/min,按照Metzner方法,即以非牛顿流体的表观粘度ηa代替牛顿流体的粘度η计算的搅拌雷诺数Re约为1 525。 根据流动的对称性,计算域选取了槽体的1/4。计算采用的网格是结构化的六面体网格40×46×36(r×θ×z),共有61 945个网格,这种网格的划分比较复杂,但在计算过程中的收敛性比较好。 2.2 计算方法 计算方法参见文献。 3 结果与讨论 3.1 叶下流场分布 图1为相同条件下水与CMC水溶液的宏观流动场。斜叶涡轮搅拌桨是介于轴流桨和径流桨之间的一类搅拌桨,它所产生的流动有时类似于轴流桨,有时又类似于径向流桨。对于宏观流动场,本文计算结果与文献的测量结果吻合。搅拌桨上方的主体循环与轴流搅拌桨所产生的循环类似,但桨叶排出流的速度方向却介于径流桨和轴流桨之间。虽然搅拌桨的离底距离已经很小,但在桨叶下部还是出现了二次循环。斜叶涡轮搅拌桨在两种介质中的总体流动形式一致,但在细节上却有较大差异。介质为水时,流体在槽体上部形成大的循环流动,流动状态较好;而介质为CMC水溶液时,自桨叶沿轴向向上,流动减弱较明显,桨叶排出的流体尚未经过上部的主体流动就在桨叶附近流回桨叶区,从而在叶端处形成一个明显的涡旋。造成这种流动状态的原因有不同的解释,湍流生长抑制理论认为高分子溶液分子的存在,改变了桨叶处近壁湍流脉动的能量平衡,从而抑制了小尺度涡旋的生成,减少了涡旋向外部的扩散,最终改变了外部的湍流结构。 图2为CMC水溶液在挡板前后的宏观速度场。从图中可以看出,在挡板以前区域,轴向运动较强,在整个r-z平面上形成一个整体循环。而在挡板后面区域,流动场较为复杂。图3给出了两个高度位置下,水与CMC水溶液在挡板后面的轴向速度分布。从图中可以看出,在桨叶高度位置水与CM