《SK静态混合器中液液两相流的标准欧拉模型和拉格朗日模.pdf》由会员分享，可在线阅读，更多相关《SK静态混合器中液液两相流的标准欧拉模型和拉格朗日模.pdf（4页珍藏版）》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、第 39 卷第 6 期 当 代 化 工 Vol.39，No.6 2010 年 12 月 Contemporary Chemical Industry December，2010 收稿日期收稿日期：2010-09-27 作者简介作者简介：张晓峰（1974-），男，辽宁沈阳人，工程师，1996 年毕业于中南工业大学有色金属冶炼专业，现主要从事制铝设备设计研发工作。E-mail：，电话：024-25505555-6161。SK 静态混合器中液液两相流的标准 欧拉模型和拉格朗日模型 张晓峰，邢玉华（中国有色（沈阳）冶金机械有限公司，辽宁 沈阳 110141）摘 要：利用计算流体动力学（CFD）模拟 S
2、K 静态混合器中液液两相流动。采用欧拉模型中的代数滑移混合模型（ASM），预测通过混合元件的压降、速度场和两相的体积分数分布。拉格朗日方法被用来追踪静态混合器中离散相的运动轨迹，用粒子的运动轨迹来分析混合的停留时间。关 键 词：石油化学；CFD；液液流动；SK 静态混合器；欧拉方法和拉格朗日方法 中图分类号：O 351.2 文献标识码：A 文章编号：1671-0460（2010）06-0714-04 Standard Euler Model and Lagrange Model of Liquid-Liquid Two-Phase Flow in SK Static Mixer ZHANG X
3、iao-feng，XING Yu-hua（NFC(shenyang)Metallurgical Machinery Co.,Ltd.Liaoning Shenyang 110141,China）Abstract:Two-phase flow in SK static state mixer was simulated by calculation flow mechanics(CFD).The algebraic slippage commixture model in Euler model(ASM)was adopted to forecast the pressure drop，velo
4、city field and two-phase volume fraction during passing the mixing element.Lagrange method was used to track the motion trajectory of discrete phase in the static state mixer，then mixing residence time was analyzed by particle motion trajectory.Key words:Petrochemistry；CFD；Liquid-Liquid Flow；SK Stat
5、ic Mixer；Euler equation and Lagrange equation 在石油化学以及相关行业中液液流动间的传质是一个很重要的设计问题。这些系统中的界面面积、液滴的大小和形成方法都是很重要的参数1。最初液液两相流的研究开始于 20 世纪 50 年代，然而仍然不能完全理解影响液滴形成的动力学因素2。由于 CFD 提供了一个对不同模型过程描述的详细分析，使得 CFD 成为模拟多项流的一个很有用的研究工具。一 个 广 泛 应 用 的 液 液 分 散 的 关 系 式 是Middleman 提出来的等式（1）。这个关系基于流体界面速度和管内径 D 提出的。作者认为，当分散相粘度低的时
6、候可以采用这两个特征量4。5/25/3032=tfWeKDd （1）这个关系式包括了韦伯数和雷诺数的影响，并涉及到静态混合内的平均平衡液滴直径。Chen 和Libby 为 SK 静态混合器提出了下面这个关系式：18.075.032)(14.1cdWeDd=（2）Manileve 和 Pacek 做出来的 SK 静态混合器内液滴大小的实验结果与等式（2）很好的吻合，这被用作为研究的根据。因此，用 CFD 模拟时采用混合器几何结构，两相流系统和它们的物理性质都与这个实验研究相似5。1 模拟对象 实验和模拟的静态混合器包括标准的 SK 型混合元件，其排列方式如图 1b。将顺时针扭转和逆时针扭转 2
7、种不同结构的混合元件交替放置在 SK 型混合器内。单一的元件如图 1a，它的直径是 25 mm，高是 37.5 mm，厚度是 3 mm。混合器结尾管长是 70 mm。入口管长度有 2 种不同结构，第 1 种是 70 mm，第 2 种是 90 mm。因此，第 1 种结构和第 2 种结构的总长度分别为 515 mm 和 535 mm。由于内管直径第 39 卷第 6 期 张晓峰，等：SK 静态混合器中液液两相流的标准欧拉模型和拉格朗日模型 715 增加了 1 mm 和长度增加了 20 mm，离散相通过第 2种结构的入口长度更长。这降低了以不同速度进入计算域的两相之间的过渡剪切力。混合管内放置了10
8、个 SK 型混合元件，第 1 个元件顺时针放置，与实验模型相似5。本文对两种情况下两相层流进行数值模拟。两种情况下的连续相都采用粘度为0.050 Pas的CMC水溶液，离散相采用动态粘度为 0.047 Pas（硅油50）和 0.479 Pas（硅油 500）的硅油。CMC 溶液，硅油 50 和硅油 500 对应的密度分别为 998，949，957 kg/m3。1a-单一元件；1b-SK 型混合元件；1c-块结构网格（77）例子；1d-非结构网格（1.5）例子。图 1 SK 型混合元件图 Fig.1 SK type mixer element 2 数值模型 利用 CFD 模拟多项流有 4 种基本
9、模型。其分别是：多孔介质模型，拉格朗日模型，欧拉模型和界面模型。本文采用了欧拉模型（ASM）和拉格朗日模型。2.1 欧拉模型 CFD 中经常用的一种欧拉模型是代数滑移混合模型（ASM）。模型假设被模拟相间之间可以相互渗透。ASM 模型是一个稳态模型。采用标准的连续方程和动量方程（N-S 方程）求解两相混合，作为一个统一体，另外也包含离散相（p）的连续方程（3）。()0,=ipppy （3）根据连续相（q）和离散相（p）所占有的空间分数，定义它们的体积分数为 p和 q，它们能取 0和 1 之间的任何值。该模型能计算两流动相的当地体积分数和速度。ASM 模型的一个基本限制就是只能用在没有传质的两相
10、模型中。另外，假设分散相的液滴大小是常数，并由 CFD 用户给出。2.1.1 模拟条件 CFD 模型的最初阶段是划分数值网格和定义边界条件。第一种结构的混合器采用 GambitTM1.2.2划分为 2 种不同类型的网格。静态混合器的 2 种网格是块结构网格和非结构网格，而且每个都有 3 种密度。块结构网格的密度分别为 55、77 和 99，它们对应的计算单元分别是 105 000、296 000和 616 000 个。非结构（四面体）网格的密度等级分别是 2.5、1.5 和 1.25，它们对应的计算单元分别是 85 000、308 000 和 614 000 个。3 种不同密度的 2 种网格形
11、成了相似的网格单元数。这既能分析网格的影响，也能分析密度的影响。SK 型混合元件的两种网格如图 1c 和 1d。采用商业软件包 CFD 的 FluentTM5.4.8 计算雷诺数分别为 100、200 和 400 时的情况。在以前已经有 36 个例子模拟了两相流，有 18 个例子也采用的是块结构网格和非结构网格，其中包括 3 种不同的雷诺数，2 种不同粘度的离散相和 3 种不同密度的数值网格。油滴直径通过给定的流体条件和方程（2）计算出来。离散相经常是以 1%的体积分数分散流动的，并且在管内的直径为 1 mm。连续相和离散相在3 种不同雷诺数下计算的速度是 q和 p如表 1。表 1 连续相和离
12、散相的速度以及油滴的大小 Table 1 Droplet size and speed of continuous phase and discrete phase under different Reynolds number 入口速度/(ms-1)油滴直径/m 雷诺数 Re 连续相速度(Vq)离散相速度(Vp)硅油50 硅油500 100 0.202 1.278 1 7532 610200 0.405 2.556 619.9922.9400 0.811 5.112 219.2326.3第二阶段是定义静态混合器的边界条件，如入口、出口和固流边界。两流体入口采用速度入口，连续相入口处油的体积分
13、数为 0，离散相入口处油的体积分数为 1。混合器出口采用压力出口，所有的固液边界采用无滑移条件。2.1.2 模拟结果 计算结果是在 Fluent 标准亚松弛因子下算得。同样，采用默认的第 1 次逼近和标准的压力速度耦合求解。将通过第 1 个混合元件压降的模拟结果与研究两相流的实验数据5进行比较。在块结构网格的模型下，硅油 50 和硅油 500 的 CMC 溶液的压降与实验值十分吻合。这些都是在雷诺数为 100 和 200 的时候 3 种网格密度下测得的，对于雷诺数等于 400时，采用最后一种网格（99），也就是含有 616 000716 当 代 化 工 2010年12月 个单元的时候测得的与实
14、验值最吻合。当采用非结构网格时，分析 SK 型混合器的压降可以得出同样的结论。因此，网格类型并没有很大的影响，只是在雷诺数超过 200 时，网格密度应该大点。2.1.3 讨论 代数滑移混合模型的一个重要缺陷就是它要求分散相单元大小是实验上的固定值。另外，ASM模型提供了多项流模型的信息，包括两相流的压降和速度，还有每相的当地体积分数。模拟表明，数值网格的类型对数值结果影响并不大。然而，在雷诺数等于 400 的时候网格密度对模拟质量的影响非常大。由于 ASM 模型的限制，作者也采用标准拉格朗日模型对两相液液流进行了模拟。2.2 拉格朗日模型 拉格朗日模型既可以用在单相流也可以用在多项流中，它是计
15、算单个粒子的流动轨迹的。Fluent 提供了一个标准的程序来分别求解连续相和离散相的流动方程。这个程序包括 3 个基本步骤，首先求解连续相的连续方程和动量传递方程，然后定义离散相粒子的初始（注入）条件，最后可视化注入粒子的轨迹。注射条件定义了初始粒子坐标，初始的粒子速度分量和初始的粒子大小 dp。注入的粒子数量仅由计算机资源所限制。在模拟的最后一个阶段，逐一求解粒子的轨迹方程和其他补充方程。Fluent 已经给出了各自的定量关系式，如等式（4）和（5）。)(ppuutu=dd （4）putx=dd （5）等式（4）和（5）整合的时间取决于在连续的时间间隔内粒子的坐标 x 和瞬时速度 up，这里
16、 1/代表粒子的松弛时间。粒子轨迹模型也能用于液液分散中追踪单个液滴，并能给出它们的轨迹、当地剪切特性和它们的停留时间。其它液体和固体的相互作用是没有标准的，对于所考虑的具体问题，需要用户自己添加CFD 代码。例如，标准代码不包括评价液滴大小演变的模型方程，它取决于剪裂和聚合。2.2.1 模拟条件 用拉格朗日模型模拟前面我们所提到的牛顿流体 CMC 溶液和硅油 50 在雷诺数等于 100、200 和400 时的情况。对第 2 种混合器结构，即入口段长度为 90 mm，在 3 个不同雷诺数和 2 种非结构网格密度下进行 6 次数值模拟。两相入口定义为轴向速度入口。速度入口分布根据文献方程计算，在
17、管中（离散相入口）和环形中（连续相入口）采用完全发展层流流动。在混合器出口，典型的条件是强加的，即常压和追踪粒子的出口。2.2.2 模拟结果 采用拉格朗日模型模拟完成，离散相注入的粒子数高达 201 个。尽管采用不同的影响系数和模拟很长的停留时间，没有那个算例能记下所有模拟粒子离开混合管的时间。模拟继续进行，直到停留时间达到 1 000 s，平均停留时间大约在 0.52.5 s之间。然而，一些粒子仍然在混合器内，仍然停留在混合器或混合元件壁面上。例如，在网格数为 164 000 个和雷诺数等于 100 的时候，201 个注射粒子中只有 20 个被采集。在雷诺数为 200 和 400 的时候，采
18、集的粒子数分别为 50 个和 72 个。这得出来的结论是在这方面我们需要进一步更细致的研究。在混合器内粒子的百分比仍然保持稳定，在网格数为 91 000 个（粗网格）时略有偏高。2.2.3 讨论 在拉格朗日模型中数值网格密度变得重要，这或许是由于所用的网格相对比较粗糙。这些主要是没有完全确定模拟方案和长的计算时间。选择拉格朗日模型模拟液液流，它给出了离散相粒子的轨迹和停留时间，并希望它能定义静态混合器中离散相的剪应力分布。SK 型静态混合器中离散相的停留时间分布结果定性上是近似正确的。这些分布情况为静态混合器的优化奠定了基础。然而，拉格朗日粒子追踪所模拟的结果似乎要进一步扩展，因为静态混合器中
19、粒子运动轨迹的不确定性。特别地，细网格的应用应该被测试，因为他们可能降低混合器内意外粒子的采集。3 结 论 这 2 种不同类型的标准模型在模拟 SK 型混合器中两相流流动时给出了不同的信息，并且有助于它们的设计。该模型可以想象流动模型、预测压降、计算当地相的速度和体积分数、评价停留时间分布和导出分散相的整个剪切历史。然而，没有一个标准模型可以预测液滴大小在静态混合器内流动的演变情况。这需要一个具体的程序来描述液滴破裂和聚合的影响，这就是预期的进一步模拟。参考文献：1 Hamad F.A.,Khan M.K.,Pierscionek （下转第 722 页）722 当 代 化 工 2010年12月
20、 证单的管理。4.4 工程项目竣工结算阶段 在这个阶段，造价人员主要任务就是依据竣工图、招标文件、施工合同、设计变更和现场签证等文件，对工程项目的造价进行审核。而审查的重点是核算实际完成工程量、审定清单未列项目单价、处理合同价格的调整问题。重中之重的是对未列清单项目定价，对此部分价格可采取这样的原则，有类似项目，可依据类似项目单价参照执行，对于其他类型项目，双方协商解决。5 结束语 工程量清单计价的推行，对石油化工行业造价管理工作提出了更高的要求，为适应市场需求，必须不断学习新的工程技术知识，学习国外造价管理中的先进经验和技术，理论和实践相结合，使石化行业工程造价管理水平上升到一个新的台阶。参
21、考文献：1 杨鲁豫，徐惠琴.建设工程工程量清单计价规范 中华人民共和国国家规范S.2003.（上接第 710 页）安全联锁的电磁阀正常时为励磁的、断电联锁电磁阀处于非励磁状态。输入旁路开关、手动联锁开关采用常开触点（开关连线断开时，不会导致连锁输入信号被旁路）。加氢装置关键部位的检测元件均采用三取二方式。例如，在反应进料泵联锁中原料油缓冲罐液位低低联锁、反应进料流量低低联锁、反应进料加热炉燃料气压力低低联锁、高压分离器液位低低联锁、循环氢入口分液罐液位高高联锁的现场一次检测元件均采用三取二方式，既保证检测信号的高可靠性，又不会因仪表本身故障造成联锁误发生。在反应进料泵出口、高压分离器至低压分离
22、器出口管线阀门采用切断阀带冗余电磁阀或控制阀、切断阀同时配置的方式。5 结束语 加氢装置联锁较多，设计完善的联锁配置来构成装置安全仪表系统，为保证装置的安全、平稳运行提供了可靠保证。对已设置的安全仪表系统要定期检查、定期更换，保证安全仪表系统的长期、有效运行。参考文献：1 李大东.加氢处理工艺与工程M.北京：中国石化出版社，2004 （上接第 716 页）B.K.,et al.Comparison of experimental results and numerical predictions of drop diameter from a single submerged nozzle i
23、n a liquid-liquid systemJ.The Canadian Journal of Chemical Engineering,2001,79(3):322-328.2 Al Taweel A.M.,Chen C.A novel mixer for the effective dispersion of immiscible liquidsJ.Chemical Engineering Research&Design,1996,74:445-450.3 Villermaux J.The role of energy dissipation in contacting and mix
24、 ing devicesJ.Chemical Engineering Technology,1988,11:276-287.4 Middelman S.Drop size distributions produced by turbulent pipe flow of immiscible fluids through a static mixerJ.I&EC Process Design and Developments,1974,13(1):78-83.5 Manileve C.,Pacek A.W.Droplet break-up within static mixers,Interna
25、l reportR.UK.Birmingham:The University of Birmingham,2001.6 Kays W.M.,Crawford M.E.Convective heat and mass transferM.New York：McGraw&Hill,1993.（上接第 719 页）比较成熟的方法；(3)重点介绍了一些适宜的利用途径，其中碳捕捉技术是一种非常有前途的 CO2处理方法。参考文献：1 颜家保，张浩，于庆满 CO2回收技术及应用前景J.应用化工，2005，34(2)：76-78 2 吴克明，黄松荣，F.Concha温室气体 CO2的分离回收及其资源化J.武汉
26、科技大学学报:自然科学版，2001，24(4)：365-369 3 费维扬，艾宁，陈健温室气体 CO2的捕集和分离分离技术面临的挑战与机遇J.化工进展，2005，24(1)：1-4 4 魏晓丹.国内外 CO2的利用现状及进展J.低温与特气，1997，4：1-7.5 宋师忠，焦艳霞 CO2用途综述及生产现状J.华工科技市场，2003，12：12-15 6 黄汉生 温室效应气体 CO2的回收与利用J.现代化工，2001，21(9)：53-57 7 江怀友，沈平平，宋新民，等世界气候变暖及 CO2埋存现状与展望J.古地理学报，2008，10（3）：323-328 8 徐俊，张军营，潘霞,等CO2存储技术的研究现状J.煤炭转化，2005，28(3)：80-85