《2022年三维搅拌器数值模拟上课讲义 .pdf》由会员分享，可在线阅读，更多相关《2022年三维搅拌器数值模拟上课讲义 .pdf（12页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、三 维 搅 拌 器 数 值 模 拟名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 1 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除搅拌器数值模拟1 引 言搅拌混合是一种常规的单元操作，具有广泛的应用背景，搅拌可以使物料混合均匀、使气体在液相中很好地分散、使固体粒子(如催化剂 )在液相中均匀地悬浮、使不相溶的另一液相均匀悬浮或者充分乳化，并可以强化相间的传质、传热。作为工业生产中工艺过程的一部分，搅拌效果直接影响到其它后续生产过程。
2、在利用超临界流体对废旧橡胶进行脱硫的课题中，脱硫反应釜中应用四叶涡轮搅拌器加强脱硫剂对溶胀橡胶的渗透作用。本文即对搅拌器在反应釜中产生的流场进行数值计算，分析搅拌流场特性，通过模拟得到流场结构及搅拌桨的速度矢量分布。2 搅拌器流场数值模拟2.1 四叶涡轮搅拌器 solidworks建模四叶涡轮搅拌器桨叶直径mm106D，叶片宽mm20a，厚mm2b，轮毂直径 20mm。三维模型建好后，保存为jiaobanqi.IGS 文件。图 1四叶涡轮搅拌器2.2 四叶涡轮搅拌器 Gambit 建模（1）将生成的 jiaobanqi.IGS文件导入 Gambit 中，得到 volume1。（2）建立搅拌槽模
3、型名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 2 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除本文采用平底圆柱形槽体，内径mm210T，槽内液位高度TH； 搅拌器安装在轴径mm16d的搅拌轴上，桨叶中心线离槽底高度3TC。图 2 搅拌槽尺寸1）建立圆柱体模型，此模型作为搅拌器的动区域，圆柱体尺寸高为60mm，半径 60mm。之后需对圆柱体进行平移，由于圆柱体的基准面都是建立在坐标原点所处的面上，本模型需使圆柱体沿着Z 轴平移，
4、设定Z 轴的平移量为-20，得到 volume2。2）以同样的方法分别建立高为40mm，半径为 8mm，高为 210mm，半径为 105mm，高为110mm，半径为8mm 的 3 个圆柱体，分别为volume3，volume4，volume5，其中 volume3无需平移， volume4 沿 Z 轴平移 -60，volume5沿 Z 轴平移 40。最终得到搅拌槽的模型如图3所示。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 3 页，共 12 页 - - - - - - - - -
5、 精品文档收集于网络，如有侵权请联系管理员删除图 3 搅拌槽模型（3）布尔运算本次模拟采用多重参考系模型( Multi-Reference Frame, MRF )。即在计算时，将计算域分成两大部分：一部分包含运动的叶片，即转子区，另一部分包含静止的槽体，称为定子区；两个区域的计算分别采用两个参考坐标系来进行，叶片所在区域（转子区）采用以叶片速度旋转的参考系，另一部分区域（定子区）使用静止参考系，具体设置见图4，参考系边界可以直接进行数据交换。a:转子区 ( 采用旋转坐标系 ) b:定子区 ( 采用静止坐标系 ) Es:两部分的界面图 4 多重参考系模型示意图因此在建立 Gambit 模型时进
6、行布尔减操作，即将搅拌槽的整体分为转子区和定子区两部分，同时在转子区去除搅拌器。首先搅拌槽整体区域与转子区相减，需保留转子区。以同样的方法完成搅拌槽整体区域与搅拌轴上部相减，转子区与搅拌器相减，转子区与搅拌轴下部相减。（4）网格划分目前，在复杂区域内生成网格的方法总的来说可以分为两大类：结构化网格和非结构化网格。当计算区域的几何结构形状比较规则时，可以实现应用结构化网格进行划分区域，在结构化网格中，每一节点及控制容积的几何信息必须加以存储，但该节点与其相邻点关系则可依据网格编号规律自动得出，因而不必存储这类信息，这是结构化网格的一大优点。当计算区域比较复杂时，即使应用专门的网格生成技术也难以处
7、理所求解的不规则区域，这时采用非结构化网格进行区域划分。在非结构化网格中，由于一个节点与其邻点的关系不是固定不变的，因此这种联结信息必须对每一个节点都显式地确定下来并加以存储。非结构化网格的这一特点对于网格的自动名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 4 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除生成、自适应处理及平行计算的实施带来不少方便，因此这种网格被普遍使用。本次模拟采用结构化网格，即六面体网格，此种网格计算过程较
8、容易，便于控制。首先对搅拌槽体即定子区进行网格划分，网格间距为2。对于转子区网格的划分，由于靠近搅拌轴，搅拌速度较快，因此转子区网格需较定子区的网格密，因此网格间距改为1，其他条件与定子区网格划分的条件相同。最后如图 5 所示为完成网格划分的搅拌槽。图 5 搅拌槽网格图（5）设定边界条件1）搅拌器和搅拌轴都设为壁面WALL 。即搅拌器Name 设为 jbq-w，Type选择 WALL；搅拌轴分为上下两个区域，靠近搅拌器的区域的搅拌轴Name 设为 shaft-down，Type同样选择 WALL ；另一区域的搅拌轴Name设为 shaft-up，Type同样选择 WALL 。2）搅拌槽底部和侧
9、面都设为WALL ，Name 为 jbc-w；搅拌槽顶部为自由液面，因此 Type选择 SYMMETRY ，Name为 jbc-top。3）搅拌槽中定子区和转子区的交接面为转子区的外表面，此表面处为两层面重合，因此选用Type 为 INTERFACE，Name 分别为 in-1，in-2，in-3，out-1，out-2，out-3，分别代表转子区外侧，转子区上表面，转子区下表面，定子区与转子区交界面的侧面，定子区与转子区交界面的上表面，定子区与转子区交界面的下表面。（7）定义实体区域如上文所述搅拌槽整体区域分为转子区和定子区，因此需对其进行设置。名师资料总结 - - -精品资料欢迎下载 -
10、- - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 5 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除将静区域 Name定为 jing，动区域 Name定为 dong，Type都选为 FLUID 。（ 8 ） 保 存Mesh文 件 ， 将 网 格 输 出 为 为jbq.msh 。2.3 Fluent求解计算以三维单精度 ( 3d )方式启动 Fluent，将在 Gambit 中建好的物理模型及网格导入 Fluent 求解器中，进行数值求解。具体求解步骤如下：（1）网格操作1）导入
11、jbq.msh 文件；在将网格导入Fluent 后，必须对网格进行检查，检查最小网格体积是否小于0,以便确定是否可直接用于CFD 求解。选择GridCheck 命令， Fluent 会自动完成网格检查，同时报告计算域、体、面、节点的统计信息。若发现有错误存在， Fluent会给出相关提示。2）设置计算区域的尺寸和单位制。Fluent 默认的长度单位为m,改为Gambit默认的单位 mm；同时需把角速度angluar velocity 的单位设置为 rpm。3）依次点击GridSwmooth/Swap Grid，分别点击Smooth，Swap 循环，直至 Swap信息中出现 Number fac
12、es swaps: 0为止；这步操作是对网格光顺以及对等角倾斜度 (skewness) 高的地方交换网格以便于后面的运算。（2）模型设置1）求解器的选择。因为后面所用的MRF 为稳态处理法，假设流动是稳定的，转子定子的作用效果是近似的平均，这种模型可用于转子定子之间的只有微弱的相互作用，或只需要求系统的近似解的场合，因此选取默认的3 维稳态求解器即可。2）选择湍流模型。选则k-epsilon2eqn，k-e 模型是最简单的完整湍流模型，是两个方程的模型，要解两个变量，速度和长度尺度。在Fluent 中，标准k-e 模型自从被Launder and Spalding 提出之后，就变成工程流场计算
13、中主要的工具了。适用范围广、经济、精度合理，所以它在工业流场和热交换模拟中广泛应用。 k-e 是个半经验的公式，是从实验现象中总结出来的，主要是基于湍流动能和扩散率。 k 方程是个精确方程，e 方程是个由经验公式导出的方程。k-e 模型假定流场完全是湍流，分子间的粘性可以忽略。标准k-e 模型只对完全是湍流的流场有效。3）设置物料特性。为简化模拟过程，本文中选用水为搅拌介质。4）操作条件的设置。由于本计算的问题需要计及重力影响，故设定Z 方向上的重力加速度分量值为-9.812sm。（3）设置边界条件。由于应用MRF 方法进行模拟，所以将转子区内的流体设定与搅拌桨相同转速进行旋转，而定子区内的流
14、体则是静止的。槽内壁面定义为静止壁面条件，搅拌轴及桨叶设为运动壁面，模型的自由液面处采用自名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 6 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除由边界条件。这里需对搅拌器、搅拌轴、定子区和转子区分别进行设置。本文搅拌器转速为rpm300N，故需在rpm300N时模拟计算。1）搅拌槽区域内的流体设置为水，因此在选择过程中都选择水为介质。2）定义搅拌器、搅拌轴、定子区、转子区的边界条件时，由于
15、搅拌过程中都随搅拌器旋转，因此选择Momentum 下选择 Wall Motion 中的 Moving Wall ，而转子区中的搅拌器和搅拌轴运动条件的设置则是选择Motion 下方的相对速度Relative to Adjacent Cell Zone 和转动 Rotational，保持 Speed为 0，即定义四叶涡轮搅拌器叶轮及转子区的搅拌轴与邻近的区城一起以300rpm 的转速同步转动。定子区的搅拌轴选择的运动条件是绝对速度Absolute 和转动Rotational，Speed为 300rpm。3）在定义转子区时需选择Motion 中 Motion Type 为 Moving Refe
16、rence Frame，即启用 MRF 模型， 在 Rotational Velocity 的 Speed里输入 300，即转子区内流体为水，且以300rpm 的转速在旋转。在定义定子区时流动形式则选择Stationary。4）除定义以上边界之外，还要定义交界面。点击DefineGrid Interfaces，选择转动转子区与静止定子区相交面（即在Gambit 里定义的 Interface边界），设置 3 组即可，分别为in1-out1z1，in 2-out2z2，in 3-out3z3，这个设置是为了使两个区域的数据能进行数据交换。（4）求解1）设置求解参数，本文基于稳态隐式分离求解算法，动
17、量方程按一阶迎风格式进行离散求解，压力速度耦合项采用 SIMPLE 算法。2）对流场初始化，选择从所有区域开始计算（Compute From all zones）。由于在稳态条件下求解流动场，因此设定初始条件为 0 。3）设定收敛残差。将各变量的收敛残差设为410，并勾选 plot 选项打开绘制参数随迭代次数变化的监视窗口。4）完成上述步骤的设置后，对文件进行保存，保存文件为jbq.cas。5）开始求解。设定迭代步数为1000，进行迭代计算，直至收敛。6）设置观测面。本次模拟观测搅拌轴所在的Y-Z 面、搅拌器所在的X-Y面上的搅拌速度云图和速度矢量图，设置X-Coordinate25 及 Z
18、为观测面。7）最后结果如图所示。3 数值模拟结果分析名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 7 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除四叶涡轮搅拌桨属于比较简单的开启式涡轮搅拌器，排出性能较高，有利于液体的快速循环流动，而剪切性能高能使流体间的速度梯度大，这都是搅拌过程所需要的。本文以搅拌槽的轴向纵截面为研究对象，在rpm300N时模拟计算，得到四叶涡轮搅拌器速度云图和速度矢量分布图。图 6 速度云图从流场的模拟结果
19、（图6，以 X-Coordinate25 为观测面）可以看出，搅拌器叶片射流处流体速度最高，随着流动远离叶片，速度逐渐降低，叶片喷射出的流体进入周围大量低速运动的流体中，卷吸周围流体，并沿轴向和径向扩散。其中代表高速的为高亮度区域。高亮度区域的面积越大，速度云图的分布范围越广，即说明搅动范围越大，搅拌效果就越好。从对四叶涡轮搅拌器的云图可以看出，搅拌器叶片附近高亮区域集中，证明搅拌器周围流畅搅拌效果最好。对定子区域部分搅拌效果的模拟不是特别理想，而且转子区和定子区之间交界面明显，分析原因可能是本次模拟过程中网格划分不够连续，网格质量不够高导致。之后的研究可以进一步针对网格划分进行。名师资料总结
20、 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 8 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除图 7 轴向速度云图图 7 为搅拌器的轴向速度云图（以X-Coordinate25 为观测面），从图中可以看出搅拌器下部靠近搅拌轴的区域，以及搅拌槽壁面附近流体的速度较高，这是由于叶轮推动液体时，流体先沿径向运动，碰到搅拌槽避免后会分别向上向下流动，即沿轴向流动，从而导致搅拌槽附近的轴向速度较大，而向下流动的流体接触到搅拌槽底部时会反转卷吸到搅拌器
21、的下部，从而造成搅拌器下部轴向速度的高亮区。图 8 径向速度云图从搅拌器的径向速度云图（图8，以 Z 为观测面）中可以看出，在搅拌器的转动下叶轮的顶部及叶片所夹区域的径向流速较大，这说明在搅拌器的转动名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 9 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除下，搅拌器周围及两叶片间的物料混合效果较好。图 9 速度矢量图从速度矢量分布（图9，以 X-Coordinate25 为观测面）可以看出，
22、四叶涡轮搅拌器呈现出对称的双涡环流场结构，隔离区中心分别位于搅拌器的上下方，符合径向流搅拌器特征。流体在叶轮出口处产生强烈的径向运动，在槽的底部和上部各产生一个循环区。在叶轮区，即叶轮边界所确定的区域，叶轮旋转直接推动液体，形成主流，液流具有射流特征，卷吸周围流体，周围的流体称为次流。中心速度随径向距离的增加而衰减，三个速度分量中以切向和径向速度为主，轴向速度可以忽略。循环流区，由于主流的卷吸作用造成液体流动，总流量可达叶轮排出流量的几倍。径向速度很小，随着远离叶轮，愈以轴向运动为主。从图中可以看出桨叶下方靠近轴心区域流体出现反向流动，这一区域称为诱导锥形区。速度矢量图比较准确地反映了搅拌器流
23、场的速度分布状态，对以后的实验有一定的参考意义。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 10 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除图 10 轴向速度矢量图图 10 为轴向速度矢量图（以X-Coordinate25 为观测面），从图中可以看出搅拌器下部速度较大，这是由于流体的卷吸作用造成的；搅拌器中心处轴向速度较小，沿着搅拌器的径向轴向速度增大。从搅拌器的径向速度矢量图（图11，以 Z 为观测面）可以看出搅拌器周围的
24、速度较大，因此搅拌器周围的物料搅拌效果较好。图 11 径向速度矢量图4 结论名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 11 页，共 12 页 - - - - - - - - - 精品文档收集于网络，如有侵权请联系管理员删除机械搅拌反应器的搅拌性能直接关系到产品的质量、能耗和生产成本。而在利用超临界流体脱硫再生废旧轮胎橡胶的实验中，反应釜中添加了搅拌系统，使反应更容易进行，因此对搅拌器进行数值模拟，详细了解搅拌器对反应釜中流场的影响是必要的。本文通过对四叶涡轮搅拌器的模拟，详细介绍了搅拌器的模拟过程，并得到了搅拌槽中流体的速度分布云图和速度矢量图，同时对其进行了分析，对以后利用反应釜搅拌器进行废旧橡胶脱硫试验有一定的指导意义。名师资料总结 - - -精品资料欢迎下载 - - - - - - - - - - - - - - - - - - 名师精心整理 - - - - - - - 第 12 页，共 12 页 - - - - - - - - -