六斜叶式搅拌器流场数值模拟

大学Zhengzhou UniversityCae课程论文六斜叶式搅拌器流场数值模拟Numerical Simulation of Shell-side Fluid-flow in the Six pitched blade stirrer专业班级：过程装备与控制工程3班作者：郝苒杏作者学号：20090360310完成时间：2012年12月16日目录摘要 (1)Abstract (1)1、背景与意义 (1)2、研究现状 (2)3、数学物理模型 (2)3.1基本控制方程 (2)3.2湍流模型介绍 (3)4、六斜叶搅拌器fluent数值模拟 (3)4.1搅拌器结构 (3)4.2几何建模 (4)4.3网格划分 (4)4.4模型求解设置 (5)4.5边界条件设置 (6)4.6残差设置 (7)4.7初始化并且迭代求解 (8)5结果分析 (8)5.1网格独立性考核 (8)5.2搅拌器流场速度矢量分析 (9)5.3搅拌器压力场分析 (10)6结论 (11)7参考文献 (11)六斜叶式搅拌器流场数值模拟摘要本文以常规六斜叶搅拌器设备为研究对象，采用数值模拟的方法，研究了搅拌器搅拌釜的流场特性的分布规律。

研究结果表明:六斜叶搅拌器流动呈现为一个位于搅拌叶片外侧的大漩涡和一个位于叶片下方的小漩涡，两个漩涡之间存在流体和能量的交换，在六斜叶搅拌器中，桨叶区湍动能较大，能量耗散率高。

将CFD技术应用于搅拌器搅拌流场的分析，基于Naives-Stokes方程和标准k－e 紊流模型，求解搅拌器的湍流场，数值模拟的结果对搅拌器水力优化设计具有指导意义。

AbstractIn this paper, numerical simulation is eateries out to study the flow fields in three stirred tanks such as the general Pitched blade turbines(PDT)，the standard RUSHTON，and a stirred equipment with special usage. The results show that there is a large-scale vortex in the outer of the blade and a small vortex below the blade. The ruction stirred is vary little flow exchange between the vortices. The region of the stirred bale has a relative large turbulence and high turbulence dissipation rate. Stirrer CFD technology is applied to the analysis of the flow field, which is based on the Naives-Stokes equations and the standard k-e turbulence model and to solve agitator turbulence field. The numerical simulation results of the agitator is helpful to guide the design of its hydraulic optimization.1、背景与意义搅拌与混合是应用最广泛的过程操作之一，搅拌设备也大量应用于化工、轻工、医药、食品、造纸、冶金、生物、废水处理等行业中。

由于相际接触面积大、传热传质效率高、操作稳定、结构简单、制造方便等优点，使得搅拌设备既可以当做反应器应用于很多场合，例如在合成橡胶，合成纤维和合成塑料这三大合成材料的生产中，搅拌设备作为反应器的约占反应器总数的85%一90%。

同时也有大量的搅拌设备并不是仅用在化学反应中应用物料的混合、传热、传质以及制备乳液、悬浮液等。

在很多化工过程中，例如水煤浆和原油的输送是煤化工，石油化的重要特征，这种高浓度的液体输送前需要有相应的搅拌过程来防止进行前可能的沉淀。

在发酵工业中，搅拌操作同样占有非常重要的地位。

发酵工业涉及到很多有氧呼吸的微生物，同时氧气在发酵液中的溶解度一般都很低。

为了保证微生物基本代活动所需要的氧气，氧气的迅速有效的供给尤为重要。

有氧发酵过程中所涉及到的搅拌操作主要是气液传质和分散。

此外，(l)发酵过程中一般都伴随有中间补给，搅拌操作可以使补给原料和基料迅速混合，避免了局部的浓度过高。

(2)微生物的代活动和搅拌过程都能产生大量的热，这些可以通过搅拌来强化传热从而使搅拌釜的物料温度保持均匀。

(3)可以使发酵液中的菌体和固体基质均匀的悬浮。

在实现混合操作的过程中，转轮的搅拌推流形式起着很重要的作用。

不同的转轮造成的搅拌推流效果差别很大，而不同的生产过程有不同的搅拌推流目的。

本文将CFD软件应用于搅拌器的搅拌流场分析，对以后的设计和分析具有指导性的意义。

2、研究现状流体混合问题实验研究方法始于十九世纪八十年代，1883年Reynolds 在水平放置的玻璃管加入染色剂，然后观察试管流体的流动揭示了层流和湍流的存在。

经过一个世纪的实验研究和理论探索，现在的流体混合技术已进去了告诉发展的时期。

到这个世纪五，六十年代流场宏观特征的测量已经普遍得到应用。

近年来，随着石油，化工等领域的快速发展，通过流场部流动的了解来优化设计的要求越来越迫切。

在计算机技术、数字技术、激光技术、图形处理技术这些相关技术快速发展的带动下，测量技术也不断更新，对流场部细节的描述也变成了研究的重点。

搅拌混合技术在20世纪中期得到了迅速发展，研究的主要问题是常规搅拌器在不同粘度的牛顿和非牛顿单相体系，固体颗粒以及小气泡在液体中混合的非均相体系中搅拌器的功率、混合时间、均匀程度等宏观特征进行试验研究。

并积累了大量的设计经验和评估体系。

但是目前将试验中的搅拌装置直接用于工厂生产中的批量或者大规模工业生产还有难度。

对于试验设备的直接放大仍然是十分危险和毫无把握的事情。

很多情况下还是采用传统的一级逐渐放大的方法，但这种方法不仅设计周期相当长，另外也会耗费大量的人力和物力。

据有关部门统计，仅仅在美国一个国家，每年因为搅拌器设计的不合理而导致的经济损失就高达几十亿美元。

另外，随着新型材料和新的搅拌技术方法源源不断的出现，对传统的搅拌设备和研究方法也带来了新的挑战。

常规的搅拌器显然已经无法满足现代工业高效率，低消耗，环境无污染等要求，而传统的搅拌器的设计方法也面临着一系列的问题。

是否能准确的描述和模拟高粘度非牛顿流体、复杂的多相流、多相藕合、各种化学反应己经成为当今极为迫切需要解决的问题。

只有从根本上有效的解决这些问题才能使搅拌器的研究设计达到一个新的高度。

近年来，光学、图像信息处理、数值计算以及计算机等领域的飞速发展使上述问题的解决成为可能。

20世纪末期激光多普勒测速仪和粒子成像测速仪等的出现，以及计算流体力学的发展使得精确测量流场结构的细节成为可能。

这些细致分析使人们读搅拌设备部的流体特性了解更加深刻，这些使得过程装备的安全和优化设计成为现实，同时也提高了计算效率和降低了失败的风险，并达到提高反应产率的目的。

在这种工业上的需要和科学技术的发展双重背景下面，使得搅拌混合技术迈进了一个全新的发展阶段。

3、数学物理模型数值模拟结果的可靠性取决于物理模型及其依靠数学手段描述现实问题的数学模型的准确性。

对管壳式换热器进行有效的数值模拟，建立准确可信的物理和数学模型是其关键。

3.1基本控制方程体流动与传热要受到物理守恒定律的支配，基本的守恒定律包括:质量守恒定律、动量守恒定律、能量守恒定律。

如果流动处于湍流状态，系统亦要遵守附加的湍流方程，控制方程是这些守恒定律的数学描述。

质量守恒方程(连续性方程)：()()()0u v w t x y zρρρρ∂∂∂∂+++=∂∂∂∂ （1） 动量守恒方程：⎰⎰⎰⎰-=+∂∂τττρρτρA A n dA n p d f dA v V d t V )( （2）能量方程：T p S gradT c k div )T U (div t )T (+⎪⎪⎭⎫ ⎝⎛=+∂∂ρρ （3）在本文中由于不牵涉热量交换，所以不用能量方程。

3.2湍流模型介绍k-ε模型是fluent 提供的一种湍流模型，它把涡粘系数和湍动能及湍动能耗散联系在一起，湍动粘度模型为：ερμμ2k C t = （4） μC 是常量。

在模型中，表示湍动耗散率的ε被定义为：⎪⎪⎭⎫ ⎝⎛∂∂⎪⎪⎭⎫ ⎝⎛∂∂=k j k i x u x u ''ρμε （5） 在标准k-ε模型中，k 和ε是两个基本未知量，与之相对应的输运方程为：()()()()⎪⎪⎩⎪⎪⎨⎧+-++∂∂+∂∂=∂∂+∂∂+--++∂∂+∂∂=∂∂+∂∂εεεεεερεεσμμρερερεσμμρρS k C G C G k C x x x u t S Y G G x k x x ku k t B K j t j i i K M b K j k t j i i 2231)(])[(])[((6) 其中，K G 是由于平均速度梯度引起的湍动能k 的产生项，B G 是由于浮力引起的湍动能k 的产生项，M Y 代表可压湍流中脉动扩的贡献，ε1C 、ε2C 和ε3C 为经验常数，k σ和εσ是k 方程和ε方程的湍流Prandtl 数，k S 和εS 是用户自定义的源项。

在本文中，常数取以下值： 09.0=μC 、44.11=εC 、92.12=εC 、0.1=k σ、3.1=εσ。

4、六斜叶搅拌器fluent 数值模拟4.1搅拌器结构六叶型搅拌桨叶片真实形状如图4-1所示。

图4-1 搅拌器桨叶片模型图4-2 搅拌器整体结构4.2几何建模模拟所用搅拌槽体为平底圆柱体，槽体直径T=300mm，液面高度H=300mm，采用六斜叶搅拌桨，直径D=T/3，也叶片倾斜角为45度，其中心线距底面距离为T/3，工作介质为水，搅拌器转速为3.53rev/s，搅拌器主要尺寸如图4-3所示。

图4-3 搅拌器主要几何尺寸4.3网格划分通过GAMBIT生成计算网格。

GAMBIT是目前CFD分析中最好的前置处理器之一，它包含功能强大的几何建模能力以及先进的网格划分工具，可以划分出包含边界层等有特殊要求的高质量的网格。

GAMBIT可以生成并处理结构化网格或者非结构化网格，主要包括的二维网格有三角形和四边形网格，三维网格有四面体、六面体、楔形和金字塔形网格。