sk静态混合器中液液两相流动的cfd模拟.docxVIP

sk静态混合器中液液两相流动的cfd模拟 原油长距离输送和相关行业的传质是一个非常重要的设计问题。这些系统中的界面面积、液滴的大小和形成方法都是很重要的参数。最初液液两相流的研究开始于20世纪50年代,然而仍然不能完全理解影响液滴形成的动力学因素。由于CFD提供了一个对不同模型过程描述的详细分析,使得CFD成为模拟多项流的一个很有用的研究工具。 一个广泛应用的液液分散的关系式是Middleman提出来的等式(1)。这个关系基于流体界面速度和管内径D提出的。作者认为,当分散相粘度低的时候可以采用这两个特征量。 这个关系式包括了韦伯数和雷诺数的影响,并涉及到静态混合内的平均平衡液滴直径。Chen和Libby为SK静态混合器提出了下面这个关系式: Manileve和Pacek做出来的SK静态混合器内液滴大小的实验结果与等式(2)很好的吻合,这被用作为研究的根据。因此,用CFD模拟时采用混合器几何结构,两相流系统和它们的物理性质都与这个实验研究相似。 1 混合管路的制作 实验和模拟的静态混合器包括标准的SK型混合元件,其排列方式如图1b。将顺时针扭转和逆时针扭转2种不同结构的混合元件交替放置在SK型混合器内。单一的元件如图1a,它的直径是25 mm,高是37.5 mm,厚度是3 mm。混合器结尾管长是70mm。入口管长度有2种不同结构,第1种是70 mm,第2种是90 mm。因此,第1种结构和第2种结构的总长度分别为515 mm和535 mm。由于内管直径增加了1 mm和长度增加了20 mm,离散相通过第2种结构的入口长度更长。这降低了以不同速度进入计算域的两相之间的过渡剪切力。混合管内放置了10个SK型混合元件,第1个元件顺时针放置,与实验模型相似。 本文对两种情况下两相层流进行数值模拟。两种情况下的连续相都采用粘度为0.050 Pa·s的CMC水溶液,离散相采用动态粘度为0.047 Pa·s(硅油50)和0.479 Pa·s(硅油500)的硅油。CMC溶液,硅油50和硅油500对应的密度分别为998,949,957 kg/m3。 2 界面模型的研究 利用CFD模拟多项流有4种基本模型。其分别是:多孔介质模型,拉格朗日模型,欧拉模型和界面模型。本文采用了欧拉模型(ASM)和拉格朗日模型。 2.1 两相混合模型 CFD中经常用的一种欧拉模型是代数滑移混合模型(ASM)。模型假设被模拟相间之间可以相互渗透。ASM模型是一个稳态模型。采用标准的连续方程和动量方程(N-S方程)求解两相混合,作为一个统一体,另外也包含离散相(p)的连续方程(3)。 根据连续相(q)和离散相(p)所占有的空间分数,定义它们的体积分数为αp和αq,它们能取0和1之间的任何值。该模型能计算两流动相的当地体积分数和速度。ASM模型的一个基本限制就是只能用在没有传质的两相模型中。另外,假设分散相的液滴大小是常数,并由CFD用户给出。 2.1.1 网格单元数计算 CFD模型的最初阶段是划分数值网格和定义边界条件。第一种结构的混合器采用GambitTM1.2.2划分为2种不同类型的网格。静态混合器的2种网格是块结构网格和非结构网格,而且每个都有3种密度。块结构网格的密度分别为5×5、7×7和9×9,它们对应的计算单元分别是105 000、296 000和616 000个。非结构(四面体)网格的密度等级分别是2.5、1.5和1.25,它们对应的计算单元分别是85 000、308 000和614 000个。3种不同密度的2种网格形成了相似的网格单元数。这既能分析网格的影响,也能分析密度的影响。SK型混合元件的两种网格如图1c和1d。 采用商业软件包CFD的FluentTM5.4.8计算雷诺数分别为100、200和400时的情况。在以前已经有36个例子模拟了两相流,有18个例子也采用的是块结构网格和非结构网格,其中包括3种不同的雷诺数,2种不同粘度的离散相和3种不同密度的数值网格。 油滴直径通过给定的流体条件和方程(2)计算出来。离散相经常是以1%的体积分数分散流动的,并且在管内的直径为1 mm。连续相和离散相在3种不同雷诺数下计算的速度是υq和υp如表1。 第二阶段是定义静态混合器的边界条件,如入口、出口和固流边界。两流体入口采用速度入口,连续相入口处油的体积分数为0,离散相入口处油的体积分数为1。混合器出口采用压力出口,所有的固液边界采用无滑移条件。 2.1.2 双相流实验结果 计算结果是在Fluent标准亚松弛因子下算得。同样,采用默认的第1次逼近和标准的压力速度耦合求解。 将通过第1个混合元件压降的模拟结果与研究两相流的实验数据进行比较。在块结构网格的模型下,硅油50和硅油500的CMC溶液的压降与实验值十分吻合。这些都是在雷诺数为100和200的时候3种网格密度下测