论文部分内容阅读
搅拌槽是广泛应用于化工、环境工程及溶剂萃取等工业领域的反应器及混合设备之一，对搅拌槽内流体力学特性的研究对于反应器的放大和设计具有非常重要的意义。关于搅拌槽湍流流场的数值模拟研究很多，但是大多采用k-ε湍流模型，只能给出时均流场，而且速度分量，尤其是湍流动能的预测值远低于实验值。
大涡模拟从理论上克服了k-ε模型各向同性假设的缺陷，适用于复杂三维湍流流动的模拟计算。本文采用改进的内外迭代法对有挡板的Rushton桨搅拌槽进行了大涡模拟计算。湍流模型分别采用Smagrinsky和dynamic大涡模拟亚格子模型。数值模拟了搅拌槽中雷诺数Re=29,000的液相流动，以水为介质。计算得到的各速度矢量和湍流动能与Wu&Patterson的实验数据比较吻合得较好，尤其是轴向速度和湍流动能比k-ε模型的结果有很大的改善，在桨区附近得到湍流动能分布的双峰曲线，而这是一般k-ε模型不能模拟出来的。文中还给出了瞬时的流场矢量图。可以发现搅拌槽内的流动非常复杂，包括大大小小的涡团，并且槽内的流场是不对称的。计算结果表明：内外迭代法和大涡模拟结合完全可以得到很好的结果。而dynamic亚格子模型的效果和Smagrinsky亚格子模型的结果差不多。
考虑搅拌槽内的流体流动具有强烈的旋流，所以在气液两相k-ε模型方程中引入强旋流修正。本文在孙海燕提出的强旋流修正的基础上从理论和模型两个方面进行了完善以后应用于标准Rushton桨气液两相搅拌槽的数值计算。结果表明：旋流修正的效果非常好，旋流修正对气含率的影响尤为突出；无论是从理论分析还是计算结果的角度，本文所采用的方法都优于从前被广泛应用的Launder的修正方法。
尽管大涡模拟模拟单相体系已经非常成功，然而涉及到两相体系的研究仍然很少。已有的一些文献基本上都是采用Lagrange观点描述分散相。本文采用Euler-Euler观点对三维搅拌槽内气液两相流动建立控制方程，然后用大涡模拟湍流模型封闭模型方程。给出了不同时刻气液两相的瞬时及时均流场矢量图，从而可以更细致地了解搅拌槽内流场的流动细节。将大涡模拟计算得到的速度值、气含率分布值和文献中的实验值、旋流修正k-ε模型计算的结果进行了比较发现：大涡模拟得到的结果和实验值吻合得最好。