鼓泡床床层高度的细网格DEM模拟

气固流化床DEM的细网格模拟中,采用传统方法计算网格空隙率和局部空隙率会对模拟结果造成较大偏差,给出一个精确面积分数模型和一个完全依赖颗粒环境的局部空隙率模型,从而更加合理地计算网格空隙率和局部空隙率。采用二倍颗粒直径的细网格模拟了小规模鼓泡流化床,模拟的气泡形状和尺寸与实验结果接近。模拟结果表明:采用给出的面积分数模型和局部空隙率模型能较好地模拟鼓泡流化床床层高度随时间变化的波形。     

基于两相两区模型的气固流化床DEM模拟

气固流化床DEM模拟中,通常采用面积加权平均法计算局部空隙率,为了考虑网格中显著非均匀结构对局部空隙率的影响,提出一个计算局部空隙率的两相两区模型。该模型将网格中的非均匀结构虚拟划分为稀相和密相,将实际网格区域划分为稀区和密区,并采用时空关联性原理识别稀区和密区;模型还采用自适应方法计算网格中颗粒分布的非均匀度;将非均匀度作为非均匀结构的影响权重计算密区颗粒的局部空隙率。模拟了鼓泡流化床,模拟结果表明:与传统的DEM相比,基于两相两区模型的DEM能更好地模拟气泡形态,且能捕捉气泡冒出床层和气泡破裂的复杂现象。     

鼓泡流化床离散模拟中的一个局部空隙率模型

气固流化床离散颗粒模拟中通常采用面积加权平均法计算局部空隙率,不能较好地反映流化床中显著的非均匀结构特性.为了考虑非均匀结构对局部空隙率的影响,文中提出一个适用鼓泡流化床的局部空隙率模型.将流场划分为稀区(气泡区)和密区(乳相区);非均匀结构对密区局部空隙率的影响通过引入局部空隙率下限和非均匀影响系数描述;将提升管流动...     

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

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

基于气泡EMMS模型的鼓泡床两相流动CFD模拟

采用两步法将气泡EMMS模型拓展到计算网格层次,用于模拟鼓泡床两相流动,改进的EMMS曳力模型不仅与空隙率有关,且与网格内的速度直接关联.求解是在操作条件下利用宏观稳定性约束条件寻优得到乳化相内空隙率关系式,再用计算网格层次上的微观守恒方程结合乳化相空隙率关系式封闭求解其余结构参数及对应的曳力系数,最后借助用户自定义函数将气泡EMMS曳力嵌入双流体模型中,对含Geldart A类颗粒的鼓泡床内气固两相流动行为进行模拟与验证.结果表明,该模型能较好地模拟鼓泡床内的非均匀流动结构,准确捕捉到上稀下浓的轴向颗粒浓度分布特征,轴向颗粒浓度分布结果与实验值较接近,平均相对偏差为6.4%.床层径向上颗粒分布均呈明显的环-核结构.在不同床层高度(0.6,0.8和1.1 m)处的径向颗粒浓度分布与实验值基本吻合,仅在底部(0.4 m)略低,相对偏差小于14%.而Gidaspow曳力模型在4个床层高度与实验数据差距甚远,相对偏差可达91%.     

一个基于双流体理论预测三维流化床内流动特性的数学模型(Ⅰ)气固体系散式或鼓泡流态化特性

将考虑拟平衡状态下颗粒与流体相互作用的附加力添加到基于双流体理论动量方程的数学模型中,用于Geldart A类物料散式流态化和B类物料鼓泡/床层塌落特性的三维数值模拟。该模型主要特点是将表征颗粒离散属性的特征长度视为颗粒直径的同一数量级且只需曳力系数一个关联式来封闭控制方程。在商业软件CFX4.4平台上,通过增加用户自定义子程序模拟了长0.2 m、宽0.2 m和高0.5 m流化床内瞬态流动特性。为了检验数学模型的实用性和数值模拟的可靠性,首先考察了两种A类物料在表观气速为u_mf和1.5u_mf下的散式流态化特性,结果展示出床层均匀膨胀的固有属性。随后,考察了扰动对A类物料在网格尺度上的局部空隙率和固体速度分布以及在设备尺度上床层压降的影响,探索了B类物料在网格尺度上鼓泡和床层塌落以及在设备尺度上鼓泡过程中床层压降和塌落过程中平均床层高度和相界面标准偏差的动态特性。上述模拟结果与经典的Geldart理论、前人的实验或模拟结果相吻合,说明该模型可以用来预报三维气固流化床内A类物料散式流态化和B类物料鼓泡及塌落的时空特性。     

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

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

网格尺寸对循环流化床内气固两相流动的影响

将基于能量最小多尺度方法(EMMS)的曳力模型耦合到双流体模型中,并针对循环流化床内的气固两流动进行了模拟研究。采用全滑移壁面边界条件处理颗粒相,考察了3种网格尺度对轴向空隙率和出口颗粒循环量等气固流动特性的影响。计算结果表明,应用EMMS曳力模型处理相间作用力,同时在采用全滑移壁面边界条件处理颗粒相时,双流体模型能够正确预测轴向空隙率分布。采用网格尺寸为2.325 mm×20 mm时,模拟结果和实测数据吻合较好,表明在循环流化床的数值模拟中选择恰当的网格尺度是极为重要的。     

三维单旋导向挡板鼓泡流化床内气固流动的CPFD模拟

单旋导向挡板是一种广泛应用于流化床反应器中的横向内构件,研究其对鼓泡流化床内气固流动的影响具有重要意义.采用计算颗粒流体力学(CPFD)方法耦合Igci曳力模型研究带有两层单旋导向挡板的鼓泡流化床内的三维气固流动,并与空筒鼓泡流化床作对比.首先将轴向和径向时均固含率分布的模拟结果和实验数据进行对比,验证了该计算方法的可...     

鼓泡床取热器流动特性的数值模拟

应用计算颗粒流体力学(CPFD)方法,采用点源注射的进气方式对气固鼓泡流化床取热器内流动特性进行数值模拟。考察不同气速下床层膨胀高度、轴径向时均固含率分布、颗粒轴向速度分布及床层颗粒内循环流率的分布。模拟结果与实验数据吻合较好,表明该模型可以用于描述鼓泡流化床取热器内气固两相的流动规律。     

A fast algorithm for mass transfer on an unstructured grid formed by DEM particles

The Discrete Element Method (DEM) was used to generate a particle packing and transient mass transfer in the particle layer was simulated using a novel algorithm for determining the mass fluxes between DEM particles. This method is intended for simulating diffusion phenomena occurring during the dissolution of swelling polymers as the DEM particles can change size and distort the grid, whilst maintaining sharp boundaries between mesh elements. In this work, the robustness and accuracy of the algorithm was tested by solving a diffusion problem over the DEM mesh and comparing it to an accurate solution obtained by a finite difference (FD) approximation, whose discretisation was 10 times finer than that of the unstructured grid. Parametric studies were conducted where the random packing and the size distribution of the DEM particles were altered and the differences in concentration profiles were compared with the FD reference solution through the use of the mean squared error. Both steady state and transient cases were compared. Two methods to improve the accuracy of the DEM unstructured grid were devised and tested using porosity and tortuosity data. In one case the porosity and tortuosity were obtained from the steady-state simulations and in the other case, local convex hulls were used to calculate porosity. Although both these methods decreased the mean squared error up to a factor of two, they also resulted into increased complexity of the simulation. It was concluded that the use of an effective packing algorithm and a narrow particle size distribution are key to maintaining the accuracy of the DEM-generated unstructured grid method.     

Searching for a mesh-independent sub-grid model for CFD simulation of gas-solid riser flows

This work aims to examine the effects of grid size in applying the two-fluid model (TFM), and thereby attempts to search for a mesh-independent sub-grid model for simulating gas-solid riser flows. To this end, we performed a series of TFM simulations over a periodic domain with various grid resolutions and drag closures. Of these drag models, EMMS/matrix model in its simplified version was chosen to be the focus of discussion. It was found that TFM simulation with a homogeneous drag model reaches its numerically asymptotic solution when the grid scale is as small as 10 times the particle diameter, but it still fails to capture the characteristic S-shaped axial voidage profile and highly over-predicts the solids flux. By comparison, EMMS/matrix model seems to reach a mesh-independent solution of the effect of sub-grid structures on the drag force, and predict successfully the axial voidage profile and the solids flux with even coarse grid. Therefore, the fine-grid TFM simulation is inadequate for gas-solid riser flows. We need sub-grid modeling of the heterogeneous structure.     

网格和数值格式对TBPBM-CFD模型模拟鼓泡塔的影响

基于气泡动力学属性的现有认识,把气泡分成大、小气泡,首次建立了完整的双气泡相-群平衡模型(TBPBM)以预测气泡尺寸.通过编写用户自定义程序实现了TBPBM模型、Luo破碎模型以及Prince 聚并模型,并耦合TBPBM与CFD双流体模型对直径440 mm鼓泡塔进行数值模拟,详细考察了网格与数值格式对TBPBM-CFD模型模拟结果的影响.结果表明,网格与数值格式对各物理变量的模拟结果影响非常大,特别是网格和体积分数方程对流项离散格式的影响最为显著.随着计算精度的提高,湍流耗散率和整体气含率分布梯度增大,气泡平均直径减少,大气泡所占气相比率降低,液相速度及气含率径向分布与实验值更趋吻合.     

Thermal conductivity of pressed powder compacts: tin oxide and alumina

Three aspects, which significantly reduce heat transfer through a polycrystalline material, are considered in this paper: porosity, grain boundary thermal resistance and the state of the grain–grain contacts. Tin oxide and alumina were chosen as model systems. Tin oxide, without a sintering additive, does not densify during thermal treatment but grain growth is not inhibited and consequently the microstructure can be varied. In alumina, variation of the thermal treatment conditions varies both grain size and porosity. Thermal conductivity measurements, using the laser-flash technique, reveal that the thermal resistance of a pressed powder compact is almost independent of temperature and at least a factor of 2.5 greater than a consolidated material with similar pore volume fraction and grain size. The reduced contact area of the grain–grain interfaces in the green body can explain this as demonstrated by numerical simulation. We also show that larger grain size increases the thermal conductivity of the porous ceramic.     

Experimental Study and DEM Simulation of Micro–Macro Behavior of TATB Granules During Compaction Using X‐ray Tomography

In this study, a noninvasive experimental method and a discrete element method (DEM) model were used to investigate the micro–macro behavior of TATB granules including global deformations, contact forces, stress distributions, local coordination number, local pore ratio, particle kinematics, and granule displacement during compacting process. The TATB granule deformation and structural evolution during compaction process were studied by X‐ray tomography and the deformation characteristics of deformable material. The size and number of granules with different sizes were determined from the CT images. The digital microstructure can be directly used for the numerical simulation setup. The DEM method was used for the numerical simulation of TATB granules during compaction. The contact forces and displacement vectors in compaction process were obtained. The behaviors and stress distributions at different pressures were also derived from the simulation. This study evaluated the possibility to better describe the behaviors of TATB granules.     

Heat transfer in aggregative and particulate liquid-fluidized beds

Wall-to-bed heat transfer in liquid fluidized beds, particulately and aggregatively fluidized, was studied. Glass particles fluidized with water gave particulate fluidization and lead particles with water gave aggregative fluidization. Local heat transfer coefficients and bed temperature profiles were measured. Heat transfer coefficients were found to be strongly dependent on particle size and porosity and increased with increasing particle size, but were independent of the height of the heater surface from the grid. Any variations in local bed properties, such as porosity do not affect wall-to-bed heat transfer. The heat transfer coefficients show a characteristic, maximum at porosities near 0.7 for both systems. Bed temperature profiles deviate considerably from open-pipe values.A two-resistance model for the heat transfer resistance agrees well with the data. Bed resistance is modeled by a radial eddy diffusivity, which indicates the mixing effectiveness in the core of the bed. Glass beds (particulate) show a maximum mixing effectiveness at porosities near 0.7 and the mixing effectiveness increases with particle diameter. Lead beds (aggregative) show two maxima in mixing effectiveness, the first between porosities of 0.5 and 0.6, and the second between porosities of 0.7 and 0.8. Mixing is greatest at an intermediate particle size in the case of lead beds. In both systems the fraction of the total resistance in the bed core increases as porosity decreases towards packed bed conditions.     

Influence of void fraction calculation on fidelity of CFD‐DEM simulation of gas‐solid bubbling fluidized beds

The correct calculation of cell void fraction is pivotal in accurate simulation of two‐phase flows using a computational fluid dynamics‐discrete element method (CFD‐DEM) approach. Two classical approaches for void fraction calculations (i.e., particle centroid method or PCM and analytical approach) were examined, and the accuracy of these methodologies in predicting the particle‐fluid flow characteristics of bubbling fluidized beds was investigated. It was found that there is a critical cell size (3.82 particle diameters) beyond which the PCM can achieve the same numerical stability and prediction accuracy as those of the analytical approach. There is also a critical cell size (1/19.3 domain size) below which meso‐scale flow structures are resolved. Moreover, a lower limit of cell size (1.63 particle diameters) was identified to satisfy the assumptions of CFD‐DEM governing equations. A reference map for selecting the ideal computational cell size and the suitable approach for void fraction calculation was subsequently developed. © 2014 American Institute of Chemical Engineers AIChE J, 60: 2000–2018, 2014     

Study of cold powder compaction by using the discrete element method

The discrete element method (DEM), based on a soft-sphere approach, is commonly used to simulate powder compaction. With these simulations a new macroscopic constitutive relation can be formulated. It is able to de-scribe accurately the constitutive material of powders during the cold compaction process. However, the force-law used in the classical DEM formulation does not reproduce correctly the stress evolution during the high density compaction of powder. To overcome this limitation at a relative density of about 0.85, the high density model is used. This contact model can reproduce incompressibility effects in granular media by implementing the local solid fraction into the DEM software, using Voronoi cells. The first DEM simulations using the open-source YADE software show a fairly good agreement with the multi-particle finite element simulations and experimental results.     

Computational study of bubble coalescence/break-up behaviors and bubble size distribution in a 3-D pressurized bubbling gas-solid fluidized bed of Geldart A particles

A computational study was carried out on bubble dynamic behaviors and bubble size distributions in a pressurized lab-scale gas-solid fluidized bed of Geldart A particles. High-resolution 3-D numerical simulations were performed using the two-fluid model based on the kinetic theory of granular flow. A fine-grid, which is in the range of 3–4 particle diameters, was utilized in order to capture bubble structures explicitly without breaking down the continuum assumption for the solid phase. A novel bubble tracking scheme was developed in combination with a 3-D detection and tracking algorithm (MS3DATA) and applied to detect the bubble statistics, such as bubble size, location in each time frame and relative position between two adjacent time frames, from numerical simulations. The spatial coordinates and corresponding void fraction data were sampled at 100 Hz for data analyzing. The bubble coalescence/break-up frequencies and the daughter bubble size distribution were evaluated by using the new bubble tracking algorithm. The results showed that the bubble size distributed non-uniformly over cross-sections in the bed. The equilibrium bubble diameter due to bubble break-up and coalescence dynamics can be obtained, and the bubble rise velocity follows Davidson’s correlation closely. Good agreements were obtained between the computed results and that predicted by using the bubble break-up model proposed in our previous work. The computational bubble tracking method showed the potential of analyzing bubble motions and the coalescence and break-up characteristics based on time series data sets of void fraction maps obtained numerically and experimentally.