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

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

液固流态化动态过程中相间作用力的数值模拟及实验验证

选择恰当的相间作用力模型是液固流态化动态特性CFD建模的关键。首先采用Richardson-Zaki关联式验证了稳态操作条件下整体固含率的实验结果，然后在基于颗粒动理学理论的欧拉-欧拉双流体模型中比较了Wen-Yu、Gidaspow、Syamlal-O’Brien、Dallavalle和TGS 5个曳力计算公式对液固流化床收缩和膨胀特性的数值模拟结果，进而探讨了Moraga等提出的升力模型影响行为及主要相间作用力影响机制。与实验测量数据比较结果表明：收缩过程中Syamlal-O’Brien和TGS曳力模型对响应时间预测较为准确，TGS曳力模型对整体固含率的预测精度较高；膨胀过程中TGS曳力模型对响应时间和整体固含率的预测优于其他模型。整体而言，基于静止颗粒群绕流直接模拟得到的TGS曳力模型忽略了颗粒-颗粒相互作用，与液固散式体系中颗粒动力学特性相符合。升力模型对动态特性模拟结果影响较小，CFD模拟时根据选择体系可予以适当忽略。     

不同曳力模型及颗粒碰撞恢复系数对短接触旋流反应器内气固流场的影响

为考察曳力模型和颗粒碰撞恢复系数对短接触旋流反应器内流动特性的影响,基于双流体模型,结合颗粒动力学理论,对反应器内气固两相流场进行模拟研究。分别采用Gidaspow、WenYu和Syamlal-O’Brien 3种曳力模型,考察颗粒速度特性以及固含率径向分布。对比分析不同曳力模型的计算结果表明,Syamlal-O’Brien模型计算结果与实验结果误差较大,WenYu模型在反应器边壁附近区域的计算结果误差较大,Gidaspow模型计算结果与实验结果最为吻合。此外,颗粒碰撞恢复系数较小时,所得计算值小于实验测量值,当恢复系数为0.95时颗粒扩散效果最好,计算结果与实验数据吻合度最高。     

不同曳力模型及颗粒碰撞恢复系数对短接触旋流反应器内气固流场的影响

为考察曳力模型和颗粒碰撞恢复系数对短接触旋流反应器内流动特性的影响,基于双流体模型, 结合颗粒动力学理论,对反应器内气固两相流场进行模拟研究。分别采用Gidaspow、Wen & Yu和Syamlal-O’Brien 3种曳力模型, 考察颗粒速度特性以及固含率径向分布。对比分析不同曳力模型的计算结果表明,Syamlal-O’Brien模型计算结果与实验结果误差较大,Wen & Yu模型在反应器边壁附近区域的计算结果误差较大,Gidaspow模型计算结果与实验结果最为吻合。此外,颗粒碰撞恢复系数较小时,所得计算值小于实验测量值,当恢复系数为0.95时颗粒扩散效果最好,计算结果与实验数据吻合度最高。     

浆态床内固含率轴向分布的数值模拟

采用DBS曳力模型计算气液相间作用,分别采用Gidaspow曳力模型、经Brucato修正的Gidaspow曳力模型和Schiller?Naumann曳力模型计算液固相间作用,忽略气固间的直接作用,对比了浆态床内不同颗粒粒径体系轴向固含率的模拟和实验结果. 结果表明,不同液固相间曳力模型对气含率的预测影响不大;在颗粒粒径较大(140 ?m)的体系中,较低表观气速下气液DBS与液固Schiller?Naumann曳力模型组合模拟的固含率随床高度增加而减小,与实验结果吻合,而其它曳力模型组合的模拟结果较差,轴向分布较均匀;在颗粒粒径较小(35 ?m)的体系中,几种曳力模型组合的模拟结果均与实验结果吻合较好,轴向分布较均匀.     

B类颗粒在鼓泡流化床中流动特性的数值模拟

采用欧拉双流体模型对鼓泡流化床中气-固两相流动行为进行数值模拟。模拟结果表明,采用结构曳力模型能够较好地预测B类颗粒在鼓泡流化床中的流动行为。对比初始流态化颗粒浓度图和完全流态化颗粒浓度分布图,可以发现结构曳力模型能够较好地展现鼓泡流化床中气泡的运动特性。当比较不同曳力模型下的模拟结果时,结构曳力模型比传统曳力模型能够更好地预测颗粒的径向分布。     

低雷诺数下气固两相中浓稀相界面对微尺度曳力的影响及建模

采用具有二阶精度的浸没边界-格子玻尔兹曼法对低雷诺数下流体流经不同颗粒聚团结构的过程进行了解析到颗粒表面的直接数值模拟(PR-DNS)。结果显示在颗粒聚团表面的浓稀相界面处,传统的微观均匀BVK曳力模型[AIChE Journal, 2007, 53(2): 489-501]的预测结果与PR-DNS结果有明显的差别;同时文献中所构建的考虑界面影响的微观曳力模型(Int. J. Multiphase Flow, 2020, 128: 103266)也无法准确预测稀相固含率不为0的情况。因此,本工作提出了一种将界面附近网格分解求取曳力的方法。通过与不同颗粒聚团结构的PR-DNS结果及其他曳力模型预测结果对比发现,新模型不但在稀相固含率趋近于0时与文献中模型具有相近的预测能力,且在稀相固含率不为0时,具有明显优于文献模型的相关系数及拟合优度。     

气-固环流反应器内瞬态流体力学特性的数值模拟

采用双流体模型和颗粒动力学理论,并考虑颗粒团聚现象对气固相间曳力的影响,对气-固环流反应器内的流体力学特性进行了数值模拟。模拟的时均固含率和颗粒速度与实验数据具有较好一致性,验证了模拟方法的可靠性。模拟结果表明：气-固环流反应器内瞬态固含率的分布具有典型聚式流态化的非均匀特征;压力脉动沿床层轴向的分布在一定程度上定性反映了气泡运动的信息;颗粒速度的时间序列和概率密度分布函数表明,床层各径向位置均存在颗粒的向上、向下运动,颗粒主体在床层内向上运动的同时还存在微观的内循环运动,模拟值为颗粒时均速度的径向非均一性宏观分布提供了合理的微观解释。     

散料颗粒密相气力输送动力参数优化及压降特性分析

为了降低管内压降,减小散料颗粒群非均匀悬浮输送所产生的湍流强度,以欧拉双流体模型为基础,应用流体力学Fluent仿真软件,建立粘度系数μ_s,k_(in)、曳力系数β在Syamlal-O’Brien和Gidaspow两种模型下形成的4种组合模型,并与实验结果进行对比,选出Gidaspow与Gidaspow组合模型作为气力输送计算模型,进而分析影响输送过程中压降变化的主要动力参数:粘度系数、曳力系数和碰撞恢复系数e_(ss),探讨前两者在该模型组合时管内压力的分布规律,得出碰撞恢复系数和入口风速v_g对管内压降变化特性的影响。     

气固下行床中壁面对介尺度曳力的影响规律

传统介尺度曳力模型的构建通常基于全周期域的细网格模拟数据集，未考虑壁面的影响，而真实流化床反应器系统中的壁面条件会影响壁面附近的非均匀结构，进而影响到壁面附近颗粒所受的曳力，因此研究壁面效应对介尺度曳力模型的影响很有必要。本文通过对拟二维全周期系统和不同床宽的周期下行床系统的介尺度曳力修正系数相对误差径向分布进行比较，探究了周期下行床中的介尺度曳力回归全周期系统曳力的现象，分析了曳力修正系数径向分布的影响因素。发现不同床宽下受壁面影响区域占整个床的比例改变较小，且曳力系数相对误差的值受该位置的固含率、固含率梯度、固相剪切率及粒化温度分布等因素影响，最后拟合了曳力系数相对误差关于径向位置的公式，并对过滤尺度对曳力系数相对误差的影响和一些变量关联曳力时是否已经隐性包含了壁面的影响进行了探讨。     

Experimental and numerical research for fluidization behaviors in a gas–solid acoustic fluidized bed

The effects of sound assistance on fluidization behaviors were systematically investigated in a gas–solid acoustic fluidized bed. A model modified from Syamlal–O'Brien drag model was established. The original solid momentum equation was developed and an acoustic model was also proposed. The radial particle volume fraction, axial root‐mean‐square of bed pressure drop, granular temperature, and particle velocity in gas–solid acoustic fluidized bed were simulated using computational fluid dynamics (CFD) code Fluent 6.2. The results showed that radial particle volume fraction increased using modified drag model compared with that using the original one. Radial particle volume fraction was revealed as a parabolic concentration profile. Axial particle volume fraction decreased with the increasing bed height. The granular temperature increased with increasing sound pressure level. It showed that simulation values using CFD code Fluent 6.2 were in agreement with the experimental data. © 2009 American Institute of Chemical Engineers AIChE J, 2010     

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

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

CFD Simulation of Mixing and Segregation of Binary Solid Mixtures in a Dense Fluidized Bed

The mixing and segregation behaviour of binary solid mixtures has been extensively studied through various experiments, while accurate CFD simulations are difficult to achieve due to process complexity and a lack of reliable constitutive relations. In this study, CFD simulations of a dense fluidized bed with glass and polystyrene particles were performed in order to identify a universal set of simulation parameters and models for simulating binary mixtures with different mixed and segregation behaviour. Through a comparison to experimental data, it was found that the EMMS drag model coupled with the Ma-Ahmadi solid pressure and radial distribution models predicted more a reasonable axial distribution of solid phases than the Syamlal O'Brien drag model coupled with the Lun et al. solid pressure and radial distribution models. The increase in the solid-solid drag further improved the simulation results.     

基于EMMS模型的搅拌釜内气液两相流数值模拟

采用欧拉-欧拉模型对搅拌釜内气液两相流进行了三维CFD模拟,重点研究了采用不同曳力模型时CFD模拟对搅拌桨附近排出流区两相流动的预测能力。模拟结果表明CFD能准确地预测排出流区的液相速度分布,但采用传统的Schiller-Naumann曳力一定程度上低估了排出流区的气液相间曳力,导致在完全扩散区CFD预测的分布器和桨叶下方区域气含率偏小,而基于气液非均匀结构和能量最小多尺度（EMMS）方法得到的DBS-Global曳力模型能更准确地描述完全扩散区气液搅拌釜内流动情况。与传统曳力模型相比,采用DBS-Global曳力模型能显著提高对气含率的预测。     

Drag models coupled with CFD–PBM method for simulation in bubble columns

In this study, three-dimensional numerical simulation of gas–liquid flow in bubble columns was realized by using the computation fluid dynamics (CFD)–population balance model (PBM). The new drag model improves the stability-constrained multi-fluid (SCMF-C) model because of the consideration of the wake accelerating and the hindering effects for calculating the drag correction factor. The gas holdup, axial liquid velocity, and bubble size distribution (BSD) predicted by four drag models at 0.02 and 0.1 m/s were compared. The results revealed that the proposed drag model can provide excellent predictions for both bubbly and heterogeneous flows. Because the wake accelerating and the hindering effects were considered, reliable predictions were achieved for the gas holdup, and the problem of uniform gas holdup distribution was mitigated. Therefore, the SCMF-C model can be extended for nonuniform BSD. The gas holdup and liquid velocity increased, and the nonuniformity of radial results became pronounced at 0.1 m/s. The profiles of four drag models were similar at a low height, whereas the difference between the simulations of the four models became obvious with the variation of heights. The results of the four models were accurate, and the BSD was wide at 0.1 m/s. Subsequently, the feasibility of the four drag models was evaluated at 0.2 and 0.4 m/s. The results of the comparison revealed that the proposed drag model exhibited excellent feasibility at higher gas velocities and was powerful for the simulation of bubble columns.     

Macroscopic properties of liquid-solid circulating fluidized bed with viscous liquid medium

Experiments were conducted in a liquid-solid circulating fluidized bed to study the effect of liquid viscosity and solids inventory on pressure gradient, critical transitional liquid velocity, onset average solids holdup, axial solids holdup distribution, average solids holdup and solids circulation rate in circulating fluidization regime with riser operated in fixed inventory mode. The results indicate that critical transitional liquid velocity decreases with increase in liquid viscosity. The onset average solids holdup, on the other hand, increases with increase in either auxiliary liquid velocity or solids inventory. The variation of axial solids holdup distribution, average solids holdup and solids circulation rate with liquid viscosity when solid inventory was 0.15 m was dissimilar with either 0.25 m or 0.35 m solid inventory. Correlations were proposed for estimating the average solids holdup and are satisfactorily compared with experimental values.     

Profiling solid volume fraction in a conical bed of dry micrometric particles: Measurements and numerical implementations

In order to develop predictive process models in conical fluidized beds, there is an ongoing search for experimental methods and simulation tools to measure and model hydrodynamics parameters. Accordingly, experiments carried out in a conical fluidized bed containing micrometric TiO₂ particles with a wide particle size distribution. An optical fiber technique was employed to determine the effect of particles loading on the solid volume fraction. The axial and radial profiles of solid volume fraction have then been determined to evaluate the sensitivity of different models, including Syamlal-O?Brien (1988), Arastoopour et al., (1990) and Gidaspow (1994) drag models. The Eulerian-Eulerian model with frictional stress models and three different boundary conditions (BCs), consisting of no-slip, partial-slip and free-slip have been used in the numerical simulations. The Gidaspow model with the partial-slip BC gives the best agreement with experimental data for different particle loadings.     

Magnetohydrodynamics of trickle bed reactors: Mechanistic model, experimental validation and simulations

An isothermal one-dimensional two-fluid magnetohydrodynamic (MHD) model based on the volume average mass, charge, and momentum balance equations and the Maxwell's equations coupled via the Lorentz force and Ohm's law was developed for the prediction of the two-phase pressure drop and the total liquid holdup in trickle bed reactors experiencing a homogeneous transverse magnetic field. The slit model approximation and the drift flux Kozeny-Carman approach were extended for the derivation of appropriate drag force closures required in the conservation equations, respectively, in the trickle flow regime and in the dispersed bubble flow regime. The expression of liquid-solid drag was adapted to take into account the influence of the magnetic field on the laminar term and the damping of turbulent/inertial term via the Hartman number and the liquid-to-bed electrical resistance ratio. Associating these drag forces with the proposed model resulted in a fully predictive MHD approach for trickle beds. Several model limiting formulations were derived for an electrically conducting fluid flowing downwards with a stagnant gas (pure trickle flow) to yield liquid holdup, as well as for single-phase upward conditions to yield the single-phase pressure drops.