基于cfd的絮凝搅拌桨水动力学特性研究.docxVIP

基于cfd的絮凝搅拌桨水动力学特性研究 废水处理是水处理技术的重要组成部分。废水处理的效果直接影响后续水质处理的效果，废水处理设备的水动力特性对废水处理效果的影响极为重要。由于工程个体的差异性,在实际设计过程中,很难做到针对具体问题提出最适合的设计值,如设计往复隔板絮凝池时,起端和末端流速、段数、转弯处过水断面积与廊道过水断面积比例、机械搅拌池中搅拌转速、桨板总面积与水流截面积的比例、桨板长宽比等。G为速度梯度,G值作为絮凝效果的评价指标也存在诸多争议。 计算流体力学CFD(computational fluid dynamics)在水处理絮凝动力学研究中的应用,对絮凝设备开发与设计产生了重要影响。CFD利用数值计算和借助计算机求解流体力学的方法,对实际设备中流动场进行数值模拟,具有精度高、耗时短、费用低,无需实验测量,可提供所需的整个流场内的有关数据等特点。本文尝试通过杯罐试验和对其内水流流态的数值模拟,研究絮凝搅拌桨的几何尺寸及其置于容器中高度这两个水动力学因素对絮凝效果的影响,提出相应的合理工况。 1 混合试验 1.1 试验计划 1.1.1 试验水样的配制 主要试验设备:ZR4-6智能型混凝试验搅拌机(六联)、圆形杯罐(1 L)、HACH2100A型浊度计、FA-2004型天平、酸度计和温度计等。 原水取自长江南京段,碱度88 mg/L(CaCO3),pH值7.86,水温20°C,色度为10。试验水样浊度为45～450 NTU,即反映全年逐月水质变化。为保证试验水样的一致性,现场取进厂原水沉淀,用静沉后的泥水与原水进行混合,配置成所需浊度的试验水样。混凝剂采用是以聚合氯化铝为核心组分的复合型聚合氯化铝(PAC),投加量为6 mg/L。 1.1.2 试验快全面转速 根据混凝控制指标要求,综合考虑杯罐试验工况,试验快搅转速600 r/min、快搅时间90 s、慢搅转速为120 r/min、慢搅时间25 min、沉淀时间20 min。 1.2 混凝剂混凝效果模拟 模拟全年原水水质,考察2种不同几何尺寸的搅拌桨对絮凝过程的影响。桨底到容器底部的距离为12 mm(装置标准配置)。桨1: 50 mm×40 mm×0.5 mm;桨2: 55 mm×20 mm×2 mm。取静止沉淀后的上清液浊度作为混凝沉淀效果的评价指标。试验结果,如表1所示。 模拟全年原水水质,考察搅拌桨入水深度的改变对絮凝的影响。采用前述桨 1,改变其在容器轴向的位置,使其上部分别距离容器底部高度H0为12、10、8 mm,比较其絮凝效果,试验结果如表2所示。 2 计算流场数值的模型 2.1 基于有限体积法的控制方程 采用FLUENT6.0软件对杯罐内流场进行三维数值模拟,计算基于有限体积法,对控制方程中的源项作线性化处理,在控制体内积分各控制方程,可得到变量在结点处的差分方程。经过差分离散后得到复杂的方程组,采用TDMA法求解,压力与速度耦合采用SIMPLE算法。 2.2 桨叶为两直叶片 数值模拟与试验容器的结构、几何尺寸一致。搅拌桨桨叶为两直叶片。容器直径D=100 mm,液面高度H=127 mm,容器高度H′=175 mm,计算域取整个容器内的液体及自由液面以上的空气。 2.3 多相流模型模型 考虑到桨叶搅拌时气—液两相的作用,采用多相流模型VOF进行模拟,定义空气为基本相而水为第2个相;容器底部z=0处和容器侧壁处r=D/2处定义为壁面边界条件。 2.4 计算值的模型 计算中采用k-ε模型,该模型形式简单并已成功模拟了许多复杂的流动问题,在数值计算中应用广泛。对牛顿流体模拟的控制方程参照文献。 2.5 外部网格 采用滑移网格法将计算域分为:内部网格——包含运动桨叶和轴的液体部分,如图1中(a)区域。内部动区域采用非结构化四面体网格,该区域梯度变化较大故采用细网格。外部网格——除内部网格以外的其他部分,图1中(b)区域所示。外部静区域用非结构化六面体网格,该区域内梯度变化小而采用粗网格,两部分网格相互独立,网格单元数共为119 944个,网格节点数共为85 093个。内部网格随搅拌桨一起转动,外部网格静止,两部分网格之间通过滑移界面进行插值处理。 2.6 采用多模计算评估指标 2.6.1 容器内充放电液流运动规律 搅拌桨为径向流搅拌桨。当搅拌桨旋转时自由液面会产生一个较小的下凹空穴,该实验现象与计算流场云图(如图2所示,颜色反映水流流速的大小)中黑线下方的空穴相一致,N为转速。容器内液体流场,如图3所示。流体在桨叶边缘处产生强烈的径向运动,遇到容器壁分别向上和向下运动,在桨的上下方都呈现出一个循环状(循环区)。在桨叶区(桨旋转所划过的圆周区域),桨叶旋转推动液体,液流具有射流特征,卷吸周围流体。表明计算结果能较可靠地反映搅拌桨旋转时流体的流动状态。 图4为容