耦合EMMS曳力与简化双流体模型的气固流动模拟

提出了一种耦合EMMS曳力的简化双流体模型,该模型忽略固相黏度,用简单的经验关联式来计算固相压力,并且耦合考虑了介尺度结构的EMMS曳力模型来计算气固相间作用力。采用简化双流体模型成功模拟一个三维实验室尺度鼓泡流化床,数值模拟结果与完整双流体模型以及实验测量结果进行了比较,结果表明耦合EMMS曳力的简化双流体模型模拟结果与完整双流体模型耦合EMMS曳力的模拟结果基本相当,并且都与实验结果吻合良好,然而简化双流体模型的计算速度是完整双流体模型的两倍以上。这表明曳力模型在气固模拟中起着主导作用,而固相应力的作用是其次的,耦合EMMS曳力的简化双流体模型在实现工业规模气固反应器快速模拟中具有巨大潜力。     

液固散式流态化CFD模拟中曳力模型的影响

基于无黏性双流体简化模型在商业化软件平台上,通过增加用户自定义子程序考察了Gidaspow、SyamlalO’Brien、Di Felice、Gibilaro、Dallavalle和BVK曳力模型对液固散式流态化CFD模拟结果的影响行为,探讨了相应的影响机制。经与文献中不同颗粒Reynolds数的代表性实验数据对比后发现:BVK和Dallavalle曳力模型对床层膨胀高度和整体固含率的预测精度较高;BVK、Syamlal-O’Brien以及Dallavalle曳力模型给出的床内固含率径向分布较为准确;BVK曳力模型较为准确地再现了颗粒轴向速度的径向分布特征。BVK曳力模型的影响机制与液固散式流态化中颗粒动力学特性相符合,在所考察范围内其预测性能最优;Dallavalle曳力模型在其余5个传统模型中预测性能较优且形式简洁在程序中易于实现。     

曳力模型和湍流模型对内环流反应器数值模拟的影响

研究了曳力模型和湍流模型对气升式内环流反应器流体力学参数的影响,进一步证实了DBS-Local曳力模型在气升式内环流反应器中的适用性。结果表明：曳力模型决定了是否可以模拟出下降管中的气体;而曳力模型和湍流模型共同作用,决定了气含率模拟结果的准确性。Schiller-Naumann、Tomiyama、Grace、Ishii-Zuber这4种曳力模型均无法预测出下降管含气这一现象,而DBS-Local曳力模型能够模拟出下降管中的气体。DBS-Local曳力模型与standard k-ε mixture湍流模型组合,对气含率的预测值与实验值较为接近,而与RNG k-ε dispersed湍流模型组合,对轴向液速的预测值与实验值更为接近。     

基于介尺度稳定性条件的多相流曳力与群体平衡模型

介尺度结构和介尺度机制是化工、冶金、能源等过程工程中的重要科学问题。尽管多相流数理模型在过去的几十年中已取得长足进展,但仍存在准确性依赖可调参数、模型适用性有限、计算量大等问题,难以适应当前快速发展的新工艺和新过程开发的需求。实际上,基于平均化方法的多流体方程组需要若干子模型封闭,如相间作用力、聚并/破碎核函数以及湍流模型等;这些子模型决定了多流体模型的模拟准确性。从介科学角度发展介尺度物理模型为解决这些问题提供了新的思路。模型可以解析多相流非均匀结构演化的控制机制,进而改进或重构子模型。总结了基于介尺度稳定条件的两类介尺度封闭模型：一类用于封闭相间动量传递,如介尺度曳力;另一类用于封闭离散相特征参数的演化,如介尺度群体平衡模型,计算气泡或液滴尺寸。进而综述了这些模型在流化床、鼓泡塔、气升式环流反应器、搅拌槽、转定子乳化器等多相流设备中的应用,并展望了未来发展方向和关键科学问题。     

液固循环流化床的研究(Ⅱ)表观滑落速度及曳力系数

在本文(Ⅰ)报实验研究的基础上,进一步探讨了液固循环流化床中的两相流动特性.液固两相表观滑落速度的实测结果表明,实验结果与广义流态化方法的预测值间存在较大的差异.表观曳力系数的研究表明,实验回归的表观曳力系数与文献中的传统预测值间也存在较大的差异,进一步体现出由于状态参数的径向不均匀性带来的循环流态化与散式流态化间的差别.     

拟泡-乳曳力模型在大型高温费托流化床反应器中的CFD模拟研究

以带冷却盘管的大型高温费托流化床反应器为研究对象,开展三维计算流体力学模拟研究。传统双流体模型基于局部平均的假设,认为单位控制体内气固两相均匀分布,网格尺寸必须足够小才能正确揭示局部非均匀结构的所有细节。采用双流体模型模拟大型工业化流化床装置时,将导致网格数量过于庞大,远超现有计算能力。为提高计算效率的同时不损失模拟精度,提出了基于局部非均匀假设、适用于粗网格的拟泡-乳三相非均匀曳力(PBTD)模型。该模型将流化床分为乳化相气体、乳化相颗粒以及气泡三相,分别建立守恒方程,体现气泡的非均匀特性对气固曳力的影响。乳化相内气固曳力以及气泡相与乳化相内颗粒的曳力分开考虑。采用PBTD模型耦合传质和反应模型,建立基于局部非均匀假设的高温费托合成反应器三维流动-传递-反应模型,包括各相守恒控制方程、气泡尺寸模型、相间物质和动量交换模型、高温费托合成反应动力学模型以及初始和边界条件,预测反应器内的流场和组分浓度分布。研究结果表明:在粗网格条件下,非均匀曳力模型可以预测床层内相含率的分布情况,预测的床层膨胀高度与经验公式计算值接近,偏差为1.2%。反应器出口气体组分的质量分数与试验测量值相近,偏差在1.5%~16.0%。模拟结果证实,基于非均匀假设的PBTD模型适用于模拟工业规模的鼓泡流化床反应器,对其设计开发和工业运行具有指导价值。     

气固流化床介尺度结构形成机制及过滤曳力模型研究进展

气固流化床中,介于颗粒与宏观尺度间的复杂的时空多尺度结构（介尺度结构）将完全改变气固相间作用规律,加大了流态化系统调控及预测的难度。为此,需要构建考虑结构影响的相间本构关系。其中,曳力作为影响流态化动力学特征的主导因素,对其研究尤为重要。从结构产生演化的机制出发,概述结构影响曳力的机理,以模型构建流程的角度对结构和过滤两类模型进行总结,并重点综述过滤模型构建在提升准确性、有效性、通用性和考虑更多物理机制方面的最新进展。研究表明：提升模型通用性和考虑真实系统中更丰富的物理机制仍是建模中亟待解决的问题,结合结构演化机制理性建模和充分发挥机器学习数据分析处理优势或是曳力建模进一步发展的关键。     

基于EMMS曳力模型的提升管进料混合段模拟

提升管进料混合段是催化裂化提升管反应器最关键的区域。为找到一种合理的方法以改善提升管进料段内油剂两相流动及混合状况,采用了EMMS曳力模型,对提升管进料混合段气固两相流混合及流动进行三维计算流体力学(CFD)模拟,并与实验数据进行了对比;分析了气固两相流动及混合特性,还模拟分析了不同进料角度对提升管进料混合段内二次流的影响以及两相流动、混合状况。结果表明,EMMS/Matrix曳力模型能够较为准确地模拟进料混合段内气固两相流动、混合过程;当喷嘴斜向上喷射进料时,射流影响区颗粒流混合不均匀,颗粒流恢复稳定流型所需时间长,且边壁受二次流影响,出现"高浓度、高返混"区域,工业过程中易引起结焦。由此提出了一种新型提升管进料段结构的改进方案,能合理利用二次流,实现颗粒流均匀混合和流动。     

曳力模型对模拟鼓泡塔气含率的影响

引言鼓泡塔由于其良好的传热、传质特性而被广泛用于化工、生物制药、冶金等领域.近年来,计算流体力学(CFD)越来越多地被应用于研究鼓泡塔内部复杂的流体力学状态.然而,如何合理地描述气液相间作用及湍流模型是CFD模拟能够准确复现鼓泡塔内复杂流动状态的关键和难点.     

固体表面活性对液固流态化的影响

用硅烷化试剂对改性聚苯醚（PPO）颗粒进行表面脱活性处理，发现其流化性能得到改善，同时水在颗粒间流动的阻力略有减小。     

A gas pressure gradient-dependent subgrid drift velocity model for drag prediction in fluidized gas–particle flows

Due to the linear correlation between the subgrid drift velocity and the filtered drag force, modeling the drift velocity would be an alternative way to obtain the filtered drag force for coarse-grid simulations. This work aims to improve the predictability of models for the drift velocity using a new effective marker, the filtered gas pressure gradient, which is identified by momentum balance analysis. New models are constructed based on conditional averaging of the results obtained from fine-grid two-fluid model simulations of three-dimensional unbounded fluidized systems. A priori assessment is presented with the comparison between the proposed models and the best available Smagorinsky-type model with dynamic adjustment technique proposed in the literature. Results show that the proposed models give satisfactory performance. More important, the proposed models are demonstrated to have a better adaptability for cases under various physical conditions than the Smagorinsky-type model.     

Three‐component solids velocity measurements in the outlet section of a riser

Coincident (simultaneous) three‐component particle velocity measurements performed using two laser Doppler anemometry probes at the outlet section of a 9 m high cylindrical riser are for the first time presented for dilute flow conditions. Near the blinded extension of the T‐outlet a solids vortex is formed. Particle downflow along the riser wall opposite the outlet tube is observed, which is restricted to higher riser heights at higher gas flow rates. Increased velocity fluctuations are observed in the solids vortex and downflow region as well as at heights corresponding to the outlet tube. Contrary to the rest of the riser, in the downflow region time and ensemble velocity averages are not equal. Given the local bending of the streamlines, axial momentum transforms to radial and azimuthal momentum giving rise to the corresponding shear stresses. Turbulence intensity values indicate the edges of the downflow region. © 2016 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 62: 3575–3584, 2016     

Fluidization of micronic particles in a conical fluidized bed: Experimental and numerical study of static bed height effect

The numerical simulations and experimental data of bed hydrodynamics in a conical fluidized bed unit are compared. Experimental studies have been carried out in a bed containing TiO₂ particles belonging to A/C boundary of Geldart's classification with a wide particle‐size distribution. Thus, pressure measurements and an optical fiber technique allowed determining the effect of static bed height on the fluidization characteristics of micronic particles. Numerical simulations have then been performed to evaluate the sensitivity of gas‐solids drag models. The Eulerian multiphase model has been used with different drag models and three boundary conditions (BC) consisting of no‐slip, partial‐slip, and free‐slip. The numerical predictions using the Gidaspow drag model and partial‐slip BC agreed reasonably well with the experimental bed pressure drop measurements. The simulation results obtained for bed expansion ratio show that the Gidaspow model with the free‐slip BC best fit with the experimental data. © 2011 American Institute of Chemical Engineers AIChE J, 2012     

Effect of anisotropic microstructures on fluid–particle drag in low-Reynolds-number monodisperse gas–solid suspensions

Local structural anisotropy prevails in gas–solid suspensions. It causes strong fluctuations in the drag on individual particles. In this work, the anisotropy of microstructures is quantified by a second-order structure tensor, which is determined with a directionally dependent mean free path length. Direct numerical simulations of low-Reynolds-number flows past anisotropic and isotropic BCC, FCC, and random arrays of monodisperse spheres in sufficiently large domains are performed. The results show that, at the same solid volume fraction, the differences between the mean drag in principal directions of anisotropic arrays and that in isotropic arrays correlate well with functions of eigenvalues of the structure tensor for the anisotropic arrays. Anisotropic drag models for different arrays are proposed. Assessment of the model for random arrays shows that it well captures fluctuations in the mean drag at microscales of several sphere diameters, where the traditional model fails to give satisfactory predictions.