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

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

Numerical investigation for the suitable choice of bubble diameter correlation for EMMS/bubbling drag model

Mesoscale bubbles exist inherently in bubbling fluidized beds and hence should be considered in the constitutive modeling of the drag force. The energy minimization multiscale bubbling(EMMS/bubbling) drag model takes the effects of mesoscale structures(i.e., bubbles) into the modeling of drag coefficient and thus improves the coarse-grid simulation of bubbling and turbulent fluidized beds. However, its dependence on the bubble diameter correlation has not been thoroughly investigated. The hydrod...     

A spatially‐averaged two‐fluid model for dense large‐scale gas‐solid flows

We present a spatially‐averaged two‐fluid model (SA‐TFM), which is derived from ensemble averaging the kinetic‐theory based TFM equations. The residual correlation for the gas‐solid drag, which appears due to averaging, is derived by employing a series expansion to the microscopic drag coefficient, while the Reynolds‐stress‐like contributions are closed similar to the Boussinesq‐approximation. The subsequent averaging of the linearized drag force reveals that averaged interphase momentum exchange is a function of the turbulent kinetic energies of both, the gas and solid phase, and the variance of the solids volume fraction. Closure models for these quantities are derived from first principles. The results show that these new constitutive relations show fairly good agreement with the fine grid data obtained for a wide range of particle properties. Finally, the SA‐TFM model is applied to the coarse grid simulation of a bubbling fluidized bed revealing excellent agreement with the reference fine grid solution. © 2017 American Institute of Chemical Engineers AIChE J, 63: 3544–3562, 2017     

Filtered and heterogeneity‐based subgrid modifications for gas–solid drag and solid stresses in bubbling fluidized beds

Two different approaches to constitutive relations for filtered two‐fluid models (TFM) of gas–solid flows are deduced. The first model (Model A) is derived using systematically filtered results obtained from a highly resolved simulation of a bubbling fluidized bed. The second model (Model B) stems from the assumption of the formation of subgrid heterogeneities inside the suspension phase of fluidized beds. These approaches for the unresolved terms appearing in the filtered TFM are, then, substantiated by the corresponding filtered data. Furthermore, the presented models are verified in the case of the bubbling fluidized bed used to generate the fine grid data. The numerical results obtained on coarse grids demonstrate that the computed bed hydrodynamics is in fairly good agreement with the highly resolved simulation. The results further show that the contribution from the unresolved frictional stresses is required to correctly predict the bubble rise velocity using coarse grids. © 2013 American Institute of Chemical Engineers AIChE J, 60: 839–854, 2014     

Heterogeneity analysis of gas-solid flow hydrodynamics in a pilot-scale fluidized bed reactor

Mesoscale drag model is of crucial significance for the reliability and accuracy in coarse-grid EulerianEulerian two-fluid model(TFM) simulations of gas-solid flow hydrodynamics in fluidized bed reactors.Although numerous mesoscale drag models have been reported in the literature,a systematic comparison of their prediction capability from the perspective of heterogeneity analysis is still lacking.In this study,in order to investigate the effect of several typical drag models on the hydrodynamic ...     

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

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

甲醇制烯烃流化床内流化特性的多尺度CFD模拟

近些年，我国成功开发了以煤为原料的甲醇制烯烃(Methanol to Olefins, MTO)生产工艺和技术，带动煤制烯烃产业的快速发展，保障了国家能源安全。流化床式反应器是MTO工业生产的核心反应装置，通过计算流体力学(Computational Fluid Dynamics, CFD)方法深入认知MTO流化床内的流化特性规律具有重要的意义，它可以从理论上更加准确地指导MTO流化床的优化与放大。本工作采用基于宏观?亚网格层次的气泡EMMS曳力和传统TFM耦合计算的多尺度CFD方法，对工业尺度MTO流化床内的多相流化行为进行了三维数值模拟。模拟结果表明，该多尺度CFD方法考虑了气泡结构对气?固相间曳力的影响，能较准确地预测MTO流化床内轴向颗粒浓度的“S-型”分布规律，且得到实验数据的验证；所预测的径向颗粒浓度分布呈现出经典的“环?核”分布规律，气体/颗粒的轴向时均速度在径向上的分布也与实际情况相互佐证，表明该多尺度CFD方法显著改善了基于均匀曳力的传统TFM对于宏观流场的预测能力。下一步工作将多尺度CFD方法拓展应用于MTO流化床优化放大及反应特性的研究。     

Gas backmixing in the dense region of a circulating fluidized bed

The gas backmixing characteristics in a circulating fluidized bed (0.1 m-IDx5.3-m high) have been determined. The gas backmixing coefficient (D_ba) from the axial dispersion model in a low velocity fluidization region increases with increasing gas velocity. The effect of gas velocity onD _ba in the bubbling bed is more pronounced compared to that in the Circulating Fluidized Bed (CFB). In the dense region of a CFB, the two-phase model is proposed to calculate D_bc from the two-phase model and mass transfer coefficient (k) between the crowd phase and dispersed phase. The gas backmixing coefficient and the mass transfer coefficient between the two phases increase with increasing the ratio of average particle to gas velocities (U_p/U_g).     

基于EMMS模型的气固鼓泡床的模拟及气泡特性的分析

基于EMMS曳力模型,采用双流体的方法对气固鼓泡床内的气固流动特性进行模拟,建立基于图像处理气泡特性的分析方法,重点研究了不同表观气速下气泡在床层内分布特性,包括气泡平均当量直径、气泡速度和气泡球形度的轴向分布,以及气泡的生命周期。研究结果表明,小气泡多集中在床层底部和壁面区域,而大气泡多集中在床层中间区域。随着表观气速的增加,床层高度不断增加,气泡的球形度降低,气泡的大小、出现频率、上升速度以及生命周期均增加;然而,当表观气速增大到一定程度,继续增加气速对气泡的上升速度影响不大。     

Acceleration of CFD simulation of gas-solid flow by coupling macro-/meso-scale EMMS model

Gas-solid flow features significant dynamic multi-scale structure; multi-scale modeling is therefore in order. In this article, the macro-scale EMMS model was coupled with a two-fluid method (TFM) elaborated by the meso-scale EMMS model resolving sub-grid scale heterogeneity to simulate the hydrodynamics of circulating fluidized bed (CFB) risers. The overall flow distribution under the steady state was approximately predicted by the macro-scale EMMS model, which serves as the initial condition for meso-scale TFM simulations reproducing the dynamic behavior of heterogeneous gas-solid flows. Using the solid circulation flux as criterion, it was shown that this coupling approach can significantly reduce the time required to reach the statistically steady state, as compared to the packed bed or homogeneously dispersed initial condition. It also suggests a general approach to speedup dynamic simulation in the multi-scale paradigm of computation.     

Kinetic simulation of a compact reactor system for hydrogen production by steam reforming of higher hydrocarbons

A bubbling fluidized bed membrane reactor for steam reforming of higher hydrocarbons is modelled, using n‐heptane as a model component to represent steam reforming of naphtha. The reformer is modelled as a bubbling fluidized bed reactor, consisting of two pseudo phases, a dense phase and a bubble phase, both in plug flow. In situ H₂ permselective membranes remove H₂ continuously as a pure product, greatly enhancing the H₂ yield per mole of heptane fed. A fluidized bed membrane reformer for higher hydrocarbons could give a very compact reactor system combining all the units from the pre‐reformer to the hydrogen purification system in a traditional steam reforming plant into a single unit.     

Effective particle diameters for simulating fluidization of non‐spherical particles: CFD‐DEM models vs. MRI measurements

Computational fluid dynamics—discrete element method (CFD‐DEM) simulations were conducted and compared with magnetic resonance imaging (MRI) measurements (Boyce, Rice, and Ozel et al., Phys Rev Fluids. 2016;1(7):074201) of gas and particle motion in a three‐dimensional cylindrical bubbling fluidized bed. Experimental particles had a kidney‐bean‐like shape, while particles were simulated as being spherical; to account for non‐sphericity, “effective” diameters were introduced to calculate drag and void fraction, such that the void fraction at minimum fluidization (ε_mf) and the minimum fluidization velocity (U_mf) in the simulations matched experimental values. With the use of effective diameters, similar bubbling patterns were seen in experiments and simulations, and the simulation predictions matched measurements of average gas and particle velocity in bubbling and emulsion regions low in the bed. Simulations which did not employ effective diameters were found to produce vastly different bubbling patterns when different drag laws were used. Both MRI results and CFD‐DEM simulations agreed with classic analytical theory for gas flow and bubble motion in bubbling fluidized beds. © 2017 American Institute of Chemical Engineers AIChE J, 63: 2555–2568, 2017     

Validation study on spatially averaged two‐fluid model for gas‐solid flows: II. Application to risers and fluidized beds

In our prior study (Schneiderbauer, AIChE J. 2017;63(8):3544–3562), a spatially averaged two‐fluid model (SA‐TFM) was presented, where closure models for the unresolved terms were derived. These closures require constitutive relations for the turbulent kinetic energies of the gas and solids phase as well as for the subfilter variance of the solids volume fraction. We had ascertained that the filtered model do yield nearly the same time‐averaged macroscale flow behavior in bubbling fluidized beds as the underlying kinetic‐theory‐based two‐fluid model, thus verifying the SA‐TFM model approach. In the present study, a set of 3D computational simulations for validation of the SA‐TFM against the experimental data on riser flow and bubbling fluidized beds is performed. Finally, the SA‐TFM predictions are in fairly good agreement with experimental data in the case of Geldart A and B particles even though using very coarse grids. © 2018 The Authors AIChE Journal published by Wiley Periodicals, Inc. on behalf of American Institute of Chemical Engineers AIChE J, 64: 1606–1617, 2018     

Particle‐resolved direct numerical simulation of gas–solid dynamics in experimental fluidized beds

Particle‐resolved direct numerical simulations (PR‐DNS) of a simplified experimental shallow fluidized bed and a laboratory bubbling fluidized bed are performed by using immersed boundary method coupled with a soft‐sphere model. Detailed information on gas flow and individual particles’ motion are obtained and analyzed to study the gas–solid dynamics. For the shallow bed, the successful predictions of particle coherent oscillation and bed expansion and contraction indicate all scales of motion in the flow are well captured by the PD‐DNS. For the bubbling bed, the PR‐DNS predicted time averaged particle velocities show a better agreement with experimental measurements than those of the computational fluid dynamics coupled with discrete element models (CFD‐DEM), which further validates the predictive capability of the developed PR‐DNS. Analysis of the PR‐DNS drag force shows that the prevailing CFD‐DEM drag correlations underestimate the particle drag force in fluidized beds. The particle mobility effect on drag correlation needs further investigation. © 2016 American Institute of Chemical Engineers AIChE J, 62: 1917–1932, 2016     

Comparative analysis of subgrid drag modifications for dense gas‐particle flows in bubbling fluidized beds

Many subgrid drag modifications have been put forth to account for the effect of small unresolved scales on the resolved mesoscales in dense gas‐particle flows. These subgrid drag modifications significantly differ in terms of their dependencies on the void fraction and the particle slip velocity. We, therefore, compare the hydrodynamics of a three‐dimensional bubbling fluidized bed computed on a coarse grid using the drag correlations of the groups of (i) EMMS, (ii) Kuipers, (iii) Sundaresan, (iv) Simonin, and the homogenous drag law of (v) Wen and Yu with fine grid simulations for two different superficial gas velocities. Furthermore, we present an (vi) alternative approach, which distinguishes between resolved and unresolved particle clusters revealing a grid and slip velocity dependent heterogeneity index. Numerical results are analyzed with respect to the time‐averaged solids volume fraction and its standard deviation, gas and solid flow patterns, bubble size, number density, and rise velocities. © 2013 American Institute of Chemical Engineers AIChE J, 59: 4077–4099, 2013     

EMMS-based Eulerian simulation on the hydrodynamics of a bubbling fluidized bed with FCC particles

Although great progress has been made in modeling the bubbling fluidization of Geldart B and D particles using standard Eulerian approach, recent studies have shown that suitable sub-grid scale models should be introduced to improve the simulation on the hydrodynamics of Geldart A particles. In this study, the flow structures inside a bubbling fluidized bed of FCC particles are simulated in an Eulerian approach employing the energy minimization multi-scale (EMMS) model (Chemical Engineering Science, 2008, 63: 1553-1571) as the sub-grid scale model for effective inter-phase drag force, using an implicit cluster diameter expression. It was shown that the experimentally found axial and radial solid concentration profiles and radial particle velocity profiles can be well reproduced.     

Particle fluctuations and dispersion in three-phase fluidized beds with viscous and low surface tension media

Particle fluctuations and dispersion were investigated in a three-phase (gas–liquid–solid) fluidized bed with an inside diameter of 0.102 m and height of 2.5 m. Effects of gas and liquid velocities, particle size (0.5–3.0 mm), viscosity (1.0–38×10⁻³ Pa s) and surface tension (52–72×10⁻³ N/m) of continuous liquid media on the fluctuating frequency and dispersion coefficient of fluidized particles were examined, by adopting the relaxation method base on the stochastic model. The fluctuations and dispersion of fluidized solid particles were successfully analyzed by means of the pressure drop variation with time, which was chosen as a state variable, based on the stochastic model. The fluctuating frequency and dispersion coefficient of particles increased with increasing gas velocity, due to the increase of bubbling phenomena and bed porosity in which particles could move, fluctuate and travel. The frequency and dispersion coefficient of particles showed local maximum values with a variation of liquid velocity. The two values of fluctuating frequency and dispersion coefficient of particles increased with increase in particle size, but decreased with increase in liquid viscosity due to the restricted movement and motion of particles in the viscous liquid medium. Both fluctuating frequency and dispersion coefficient of particles increased with decrease in surface tension of liquid phase, due to the increase of bubbling phenomena with decrease in σ_L. The values of obtained particle dispersion coefficient were well correlated in terms of dimensionless groups as well as operating variables.     

基于EMMS的循环流化床流域研究

流化床的设计、放大和优化需要对流域有基础的认识，然而气固系统的流域划分至今仍存在诸多争议。总结了气固流化系统流域划分的研究现状，并分析了流域划分的主要争议，发现文献中对快速床的界定存在分歧。通过耦合基于稳态EMMS的曳力模型开展双流体模拟，对不同气速和颗粒浓度下的循环流化床进行了数值研究。模拟结果捕捉到了颗粒回流、节涌等现象，据此确定了快速床的操作边界并绘制了流域图，该流域图能够展示循环床中的各流域形态。     

Effect of annular fins on heat transfer of a horizontal immersed tube in bubbling fluidized beds

Experiments were conducted in a bubbling air-fluidized bed to investigate the effect of annular fins of constant thickness on heat transfer. Steady state time averaged local heat transfer coefficient measurements were made by the local thermal simulation technique in a cold bubbling fluidized bed (90 mm ID, 260 mm tall) with horizontally immersed tube initially with no fin and then with three fixed annular fins of constant thickness. Silica sand of mean particle diameter 307 μm and 200 μm were used as the bed materials. The superficial velocity of air was from minimum fluidization conditions, u_mf, to approximately 3 × u_mf. The results indicate that, although the heat transfer coefficient falls with the use of fins, the total heat transfer rises as a result of the greater surface area. Increasing the particle diameter reduces the heat transfer coefficient not only for unfinned horizontal tube but also for annular finned horizontal tube at the same conditions of fluidized bed. Based on the experimental data, correlations are proposed for predicting heat transfer coefficient from fluidized bed to horizontally immersed tubes with and without fins.     

粗糙颗粒动理学及流化床内气固流动的数值模拟

基于气体分子动理学和颗粒动理学理论,考虑颗粒旋转流动对颗粒碰撞能量交换和耗散的影响,建立粗糙颗粒动理学。采用Chapman-Enskog颗粒速度分布函数,提出了颗粒相应力、热通量和颗粒碰撞能量耗散计算模型。采用欧拉-欧拉气固双相流模型,数值模拟鼓泡流化床内气体-颗粒两相流动特性。模拟结果得到了床内颗粒相速度和脉动速度分布,与Yuu等实验结果相吻合。分析不同的切向弹性恢复系数对颗粒相拟总温的变化规律,结果表明在低颗粒浓度时颗粒拟总温随切向弹性恢复系数而增加。