颗粒堆积床内两相流动阻力及相间摩擦力

为了探究颗粒堆积多孔介质通道内两相流动阻力特性,使用不同尺寸的球形颗粒构建实验床,在碎片床冷却性实验装置(DEBECO)上进行了气-水两相流动实验.基于实验数据,对比验证了多孔介质内两相流动阻力压降预测模型,重点分析了颗粒床内气液两相相间摩擦力对阻力压降的影响.结果表明:对小尺寸(如d=1.5mm)颗粒堆积床,气液两相相间摩擦力对其两相流动阻力压降影响微弱,两相流动阻力压降随气速增大呈逐渐上升的趋势;对大尺寸(如d=6mm)颗粒堆积床,在低速下气液两相相间摩擦力对阻力压降影响显著,随气速增大,其阻力压降出现先下降后上升的趋势.     

磁流化床气固两相流动的数值模拟及实验验证

建立了床体为圆筒的磁流化床筒化数学模型。采用SIMPLEC算法对此类磁流化床床内的气固两相流动进行数值模拟，并对多种运行条件下床内气固两相流动的数值模拟结果进行分析，获得了磁流化床内气固两相在偏离或接近临界稳定区域的流化特性。通过用相对应的实验装置进行实验验证。结果证明：该简化的双相模型和SIMPLEC算法对磁流化床床内气固两相流场进行数值模拟具较好的适用性和准确性。     

蜡沉积凝析气-液-固变相态复杂渗流数学模型

含蜡凝析气藏开采过程中伴有凝析液、蜡析出呈气-液-固变相态多相复杂流动. 在渗流机理的实验研究基础上,对气-液-固复杂渗流的数学模型进行描述,以期揭示伴有蜡沉积的凝析气-液-固变相态复杂渗流规律. 根据相的变化受析蜡线、露点线的相态变化规律控制,分析相变特性会导致储层内油气饱和度的变化,揭示蜡沉积在多孔介质表面引起孔隙性质改变和气、液、固在储集层中的微观空间分布会影响凝析气的流动特征. 在蜡沉积变相态渗流实验和渗流数学模型的研究基础上,引入蜡沉积数学模型、考虑毛管力和表面张力的相对渗透率模型等,建立了相应的多孔介质中气-液-固变相态多相流-固耦合渗流数学模型. 数值模拟结果表明,蜡沉积可使部分区域孔隙减少,流动阻力增大. 不考虑蜡沉积会对产量预测明显偏大,蜡沉积可使产量下降. 该模型较好地反映了凝析气-液-固变相态复杂渗流的物理实质.     

多喷嘴对置式煤气化炉内气固两相流动及壁面沉积

对我国具有自主知识产权的多喷嘴对置式煤气化炉内的气固两相流进行了数值模拟分析.在合理简化和假设基础上建立了基于Eulerian-Lagrangian模拟方法的炉内气固两相流动模型,采用Realizable k-ε模型描述炉内复杂气相湍流运动,应用颗粒轨道模型随机追踪煤粉颗粒在湍流气流中的运动.通过数值计算获得了炉内气固两相的速度矢量、颗粒分布、颗粒运动轨迹,以及颗粒碰撞并沉积于壁面的通量分布.揭示了该型气流床煤气化炉内气固两相流动特征,并分析了入口速度和炉体上部高度对气固两相流动和颗粒在壁面沉积的影响规律.结果表明:对撞流显著增强炉内气固流动湍动,强化气固相互作用,并使煤粉颗粒在炉内有效分散,有利于化学反应高效进行;在喷嘴入口及顶部的壁面处颗粒沉积率较大.     

旋转圆盘上液固两相流冲刷磨损数值模拟研究

对旋转圆盘内液固两相流冲刷腐损过程进行了数值模拟研究.数值模拟研究采用欧拉-拉格朗日方法来模拟液固两相流动,即通过在欧拉坐标系下求解流体相的雷诺时均方程组来模拟流体流场,通过拉格朗日坐标系下的随机轨道模型来获得固体颗粒相的运动轨迹.将最终数值模拟结果和实验测量数据进行了比较,结果表明文中构建的综合模型基本上是正确和可行的.     

滴流床加氢裂化反应器内流体流动的数值模拟

采用计算流体力学(CFD)方法,利用Fluent软件中多孔介质模型对某实验室规模滴流床反应器进行数值模拟;对比考察3种多孔区域三相动量相互作用模型对计算结果的影响;通过模拟计算得到具了反应器内速度场、压力场及持液量的详细信息,并对压降和持液量进行分析对比.结果表明,模拟计算结果与文献中实验报道结果相符,在3种动量相互作用模型中双流体界面力模型能描述床层的微观流动特征(壁流).     

CuO/Al2O3作为载氧剂的流化床化学链燃烧数值模拟

以欧拉-欧拉双流体模型和气固非均相化学反应动力学为基础,嵌入了气固化学反应速率方程和反应内热源项的UDF(自定义函数)程序,对化学链燃烧燃料反应器——鼓泡流化床内气固两相流动及化学反应过程进行了数值模拟,并分析了甲烷进气速度对床内气固两相流动、传热及化学反应速率的影响.结果表明:随着甲烷进气速度增加,床内气固混合更加剧烈,气泡的产生、碰撞和破碎使得气固分布不均,流化质量下降,导致反应器内化学反应速率以及温度分布不均,床内局部存在的高温区域将使颗粒温度过高而烧结,降低了甲烷燃烧效率.     

基于CFD-DEM耦合并行算法的锥形喷动床内离散颗粒数值模拟

基于FLUENT软件的MPI平台,在二次开发框架内通过用户自定义函数文件实现了CFD-DEM耦合并行算法.该并行算法具有随CPU核数增加的较佳扩展性能和良好加速性能.通过构建空隙率标量场、重组基于局部平均并考虑气-固相互作用的气相控制方程,提高了气相控制方程的求解稳定性与收敛性;通过软球模型确定颗粒-颗粒及颗粒-壁面间的相互作用,以曳力模型为基础建立了相间相互作用关系.将该CFD-DEM耦合并行模型应用于锥形喷动床气固两相流动的数值模拟中,结果表明:该CFD-DEM耦合并行模型能较好地模拟锥形喷动床内气固流动行为;锥形喷动床内的颗粒流动经历了鼓泡阶段和稳定流态化阶段;在稳定流态化阶段,物料流动具有明显的周期性.基于锥形喷动床内颗粒相的瞬时速度场,获得了稳定流态化阶段颗粒的时均速度分布规律.     

高炉三维气固湍流和煤粉燃烧过程数值模拟

基于欧拉气相方程组、欧拉颗粒连续方程和动量方程以及拉氏颗粒能量和质量变化方程,建立并发展了高炉风口回旋区湍流气固两相流动和煤粉燃烧的三维数学模型.用所建模型分别对冷态模型内气固两相流动和某企业750 m3高炉风口回旋区内的气固两相三维流动与煤粉燃烧进行了数值模拟.采用三维激光相位多普勒分析仪(PDA)对冷态模型内气固两相流场进行了测量,实验结果与冷态两相流动的模拟结果基本一致.热态模拟结果给出了气相温度和组分浓度分布,模拟结果与实验测量结果较吻合,揭示了风口回旋区内气固两相流动和煤粉燃烧的基本性质和特点.     

深海采矿提升泵的数值模拟分析

针对深海底矿物粗颗粒浆体在提升泵内的流动问题，提出了一种新型粗颗粒-均质浆体两相流动模型．采用流体动力学软件CFX对两级提升泵内的粗颗粒固-液两相湍流进行了数值模拟.计算获得了粗颗粒在提升泵内的流动特征及其泵工作特性的仿真结果，并与两级提升泵的性能试验结果进行对比，两者能较好地吻合，从而验证了该两相流动模型和数值模拟的有效性和准确性.     

并联管路气液两相流分流实验与数值模拟

为了研究气液两相流在并联立管中的分流状况,在实验室中建立了带有并联立管的气液两相流流动回路。建立了气液两相流分流数学模型,通过数值模拟方法对气液两相流分流进行研究。实验过程中发现了并联立管中三种不同的流动状态分别为弹状流-弹状流、搅混流-搅混流以及弹状流-搅混流。实验结果表明当并联立管中的流动状态为弹状流-弹状流时,两根立管中的压降相同,但是并联立管中的液塞排出并不具有同步性,可能产生偏流状况。在气液流速较高的搅混流-搅混流时两相流在并联立管中始终处于均匀分配状态。介于两者之间的弹状流-搅混流,并联立管中的压降产生较大的偏差,气液两相流流量更倾向于搅混流一侧流出,此时并联立管中产生严重的偏流状况。数值模拟过程中,通过对并联立管出口持液率进行监测,在较低的气液速下两相流并不能够均匀分配,但是在较高的气液速下,两相流始终处于均匀分配,可见所建立的数学模型具有良好的准确性。     

液力缓速器瞬态两相流动大涡模拟及性能预测

为了全面掌握液力缓速器各个流动单元内速度场和压力场的分布特性,提取了全流道几何模型作为计算区域,利用CFD平台的大涡模拟法和多流动区域耦合计算的滑动网格法对液力缓速器内部气液两相流动进行了三维瞬态数值模拟,将混合模型与欧拉模型交替运用在其多相流模型中,获得了不同充液率下液力缓速器内流场结构的变化及两相体积分数分布情况,分析了流场内二次流、脱流及涡旋的产生机理,并计算了缓速器的外特性.结果表明:数值计算结果与试验结果吻合很好,误差在8％以内;流场计算十分准确,运用的大涡模拟方法可以有效地模拟液力缓速器流场内的真实液流结构,其结果可以用来指导液力缓速器的设计及其结构优化.     

液力偶合器两相流动的研究

应用先进的CFD平台的滑动网格法与混合模型对液力偶合器部分充液的两相流场进行数值模拟计算和分析,得到偶合器部分充液时内部速度和压力分布以及不同充液率下流场结构的变化。基于速度、压力数值解计算了液力偶合器叶轮转矩和不同充液率下液力偶合器的原始特性。将性能预测结果与实验结果进行对比分析,结果表明,数值模拟计算结果与实验结果吻合较好,从而验证了数值模拟方法的正确性。     

预测气液两相分层流界面剪切应力的新方法

针对水平管气液两相分层流界面剪切应力的预测,建立了基于计算流体力学（CFD）方法的剪切滑移壁面模 型,即将气液界面处理为固定的具有均匀剪切应力的水平滑移壁面。利用FLUENT 软件模拟了气相流动,气液界面采 用剪切滑移壁面边界条件,通过对比实验测量的气相流场分布得到气速和界面剪切收敛值,并由收敛值拟合得到界 面摩擦因子与气相雷诺数的关联式。应用这一关联式的预测结果与实验间接测量值吻合较好,并且在气、液相雷诺 数9 4006ReG650 000 和21 0006ReL630 000 范围内优于Taitel 和Dukler 与Sidi-Ali 和Gatignol 模型。表明剪切滑移壁 面模型可以有效地控制关键参数,提高预测效率和精度,为CFD 方法预测气液两相分层流界面剪切应力提供了新的 思路。     

超重力旋转床中气液两相流动与传质过程的数值模拟

文中采用基于颗粒轨道模型的欧拉-拉格朗日法对超重力旋转床中的气液两相流动与传质进行了数值模拟研究。在合理简化丝网填料结构和考虑液滴凝并与分散的基础上,分别利用SIMPLE算法和颗粒轨道模型计算了超重力旋转床中的气流场和液滴的运动轨迹,进而计算了液相的传质系数。数值模拟所得的液相传质系数与氮气解吸水中溶解氧实验结果符合良好,表明模型能够用于模拟旋转床中流体力学和分散相内的传质过程。计算分析表明,对超重力旋转床,在一定的转速下,液体和气体流量以及填料内径的变化对体积传质系数有重要影响。     

大型筛板塔上两相湍流的数值模拟

运用CFD数值模拟方法，建立了大型筛板塔塔板上气液两相流双流体模型，数值求解了三维速度场分布。文中应用κ-ε模型封闭雷诺方程，构成了二方程湍流模型，附加关系采用了修正GRACE相间曳力模型，CFD模拟结果与实验值进行了对比，真实地反映了塔板上三维混乱的流体力学特性，并揭示了液相的回流现象。     

气液两相状态下有衬里水带水动力特性研究

消防水带内工质为有压气液两相流状态时,其工作过程中会产生无规律振动,极易造成水带磨损和水枪的难以掌控,影响灭火救援工作的开展。对气液两相状态下消防水带水动力特性展开研究,通过数值模拟和实验的方法探索水带内气泡流动、合并、破碎、脱离的过程,分析气相比例、运行压力对水带稳定性的影响。数值计算中考虑重力、表面张力和壁面黏附作用。研究结果表明:气相比例是导致水带不稳定运行的重要因素,但存在一个极大值;运行压力的增加会直接导致水带运行不稳定性的增大,但两者之间未发现定量关系;气液两相的剧烈扰动,是造成气泡的流动、合并、破碎的主要诱因,在气液两相流状态下,应尽量避免水带高压运行。     

液力缓速器部分充液流场大涡模拟及特性预测

利用CFD平台的大涡模拟法和滑动网格法对液力缓速器内部气液两相流动进行三维瞬态数值模拟,将混合模型与欧拉模型交替运用在其多相流模型中,得到速度场和压力场的分布特性,总结不同充液率下液力缓速器内流场结构的变化及两相体积分布情况.基于流场数值解计算了液力缓速器不同工况下的制动转矩,进而预测其性能,性能预测结果与试验结果吻合较好,误差在8％以内,表明流场计算是较准确的,其结果可用于液力缓速器的设计及其结构优化.     

旋流条件下的气液环空流流动规律研究

针对涡流工具排液效果的问题,开展了旋流条件下气液两相流动模型的研究。考虑到旋流中角速度的存在,研究中采用气液流动在径向和周向上的动量和角动量平衡的方法,建立了气液流动控制方程,计算了液膜厚度,气液相旋流强度等参数以及压降梯度,并进行涡流工具实验验证模型。涡流工具降低压降损失的机理结果表明,安装涡流工具后流动压降可以降低5%~20%。根据实验及模型,在低速（气相速度小于13 m/s）时,小旋流角和高旋流强度更利于降低压降,而高速（气相速度大于16 m/s）时,高旋流强度会增加额外摩擦阻力。旋流强度的衰减速度会随着液相速度增大而减小,而随气相速度增大而增大。该研究结果可对涡流工具进行优化设计,以达到最佳排液效果。     

孔板通气气液两相流动特性研究

采用实验与数值模拟相结合的方法研究了孔板通气气液两相流动特性和气层形成机理.应用高速录像观察了通气气液两相流流动形态,采用Level Set方法和RNG k-ε湍流模型对通气平板流场进行了数值计算.结果表明气层覆盖区域压力是判定气层形成的重要参数:泡间高压消失表明当地的气泡开始形成气层,当气层基本形成时,壁面压力趋于一致.旋涡结构在气层形成不同阶段具有明显的特征,在气层初始形成位置前后,肾形涡由于气泡接触相互作用而逐渐消失;在气层完全形成位置前后,旋涡结构重新生成.随着孔间距的增加,气层形成位置由于横向距离增加而向下游移动,但相互作用区域变化不大.