SK型静态混合器的数值模拟_方志林.docVIP

第 23卷 第 4期 广东石油化工学院学报 Vol. 23No. 4 2013年 8月 Journal of Guangdong University of Petrochemical Technology Aug. 2013 SK 型静态混合器的数值模拟 * 方志林 (巢湖学院 化生学院 , 安徽 巢湖 238000 摘要 :利用专业的多物理场耦合软件 COS MOL Multiphysics 对 SK 型 静态混 合器内 的不可 压缩层 流流场进 行了仿 真模拟 , 计 算过程中采用软件中的层流流动和稀物质的传递内置模块 , 得到了静态混合物中流体的浓度 场、 速度 场及压力 场分布。由 于螺旋混合元件的作用 , 流体做螺旋运动 , 通过对流体的剪 切、 旋 转 , 使流体 充分混合 , 结 果在出口 处几乎 混合完 全。利用 后处理软件对计算结果进行分析 , 为其性能测 试提供了理论依据。 关键词 :SK 型静态混合器 ; COS MOL Multiphysics; 数值模拟 中图分类号 :TQ051.7文献标识码 :A 文章编号 :2095-2562(2013 04-0084-03 SK 型静态混合器又称单螺旋形静态混合器 , 它的单元是扭转 180b 或 270b 的螺旋板 , 安装入管壳时相 邻的螺旋板分别为左旋和右旋。当流体逐次流过每个元件时 , 即被分割成越来越薄的薄片 , 其数量按元件 数的幂次方增加。由于混合单元的使用 , 使流体时而左旋 , 时而右旋 , 不断改变流动方向 , 不仅将中心液流 推向周边 , 而且将周边流体推向中心 , 从而造成良好的径向混合效果。与此同时 , 流体自身的旋转运动在 相邻元件连接处亦会发生 [1]。随着计算机技术的发展和计算流体力学方法的日趋成熟 , 越来越多的学者 采用计算流体动力学 (C FD 方法对 SK 型静态混合器的流体力学行为进行研究 , 大多数报道的结果显示 SK 型静态混合器的流场由雷诺数混合元件的结构及个数等条件决定。但在模拟计算时由于流体的流动 会影响到流体混合的效果 , 也就是说流体的速度场与浓度场会有耦合 , 常规的 CFD 模拟不能同时得到几 个物理场的同时模拟。 COMSOL Multiphysics 是以有限元法为基础 , 通过求解偏微分方程 (单场 或偏微分 方程组 (多场 来实现真实物理现象的仿真 , 被当今世界科学家称为 /第一款真正的任意多物理场直接耦合分 析软件 0。该软件的内置应用模块功能为这一物理耦合的研究提供了可能 [2]。 1SK 型静态混合器装置 COSMOL Multiphysics 仿真要经历建立选择空间 维度、 增加物理量、 几何模型、 定义物理参数、 划分有 限元网格、 求解、 可视化后处理等过程 [3]。 SK 型静态混合器的几何模型如图 1所示 , 它由 一个管道组成 , 管道具有三个交替旋转的扭曲叶片 , 这种混合技术特别适合层流混合 , 因为在流动区域 它只产生小的压力损耗。 图 1三个交替旋转叶片的静态混合器 2SK 型静态混合器仿真模拟 2. 1边界条件 进口边界条件为速度入口 ; 出口边界条件为自由流出口 ; 面默认为固壁无滑移光滑界面。层流时的速 *收稿日期 :2013-05-06; 修回日期 :2013-05-28 :( , 男 , , : 度沿管径按抛物线的规律分布 , 在管壁处速度为零 , 离开管壁以后速度渐增 , 到管中心处最大 , 其速度分布 方程为 : u =u max (1-r 2/r 2 (1 其中 , u 为管内任意处的速度 , m P s; u max 为流体在圆管中的最大流速 , m P s; r 为管内任意质点处到轴心距 离 , m; R 为管内半径 , m 。以溶解在水中的物质 (CaSO 4 的混合过程进行模拟。已知进入静态混合器的平 均流速 u -为 0. 01m P s, 入口的浓度为 5mol P L, 粘度 L 为 1@10-3pa #s, 密度 Q 为 1000kg P m 3, 扩散系数 D 为 5@10-8m 2P s 。雷诺数按 Re =d Q u P L 计算为 120, 流动状态为层流。 2. 2 数学模型 求解流体流动控制方程可以表示为 [4]: 9(?Q +di M (Q u ? =di M (#grad ? +S (2 其展开形式为 : +++=1#+1#+1#2+S (3 其中 :?为通用变量 , #为广义扩散系数 , S 为广义源项。 稀物质的传递方程可以表示为 : ¨#(-D ¨c+cu =0 (4 2. 3 物理模型 SK 混合叶片将所在的管道分为两部分 , 流体被迫沿叶片流动。其中管道的半径 R =6mm, 管道长度 84mm, 混合器旋转