Modeling of particle transport and combustion phenomena in a large-scale circulating fluidized bed boiler using a hybrid Euler-Lagrange approach

The constantly developing fiuidized combustion technology has become competitive with a conventional pulverized coal （PC） combustion. Circulating fluidized bed （CFB） boilers can be a good alternative to PC boilers due to their robustness and lower sensitivity to the fuel quality. However, appropriate engineering tools that can be used to model and optimize the construction and operating parameters of a CFB boiler still require development. This paper presents the application of a relatively novel hybrid Euler-Lagrange approach to model the dense gas-solid flow combined with a combustion process in a large-scale indus- trial CFB boiler. In this work, this complex flow has been resolved by applying the ANSYS FLUENT 14.0 commercial computational fluid dynamics （CFD） code. To accurately resolve the multiphase flow, the original CFD code has been extended by additional user-defined functions. These functions were used to control the boiler mass load, particle recirculation process （simplified boiler geometry）, and interphase hydrodynamic properties. This work was split into two parts. In the first part, which is referred to as pseudo combustion, the combustion process was not directly simulated. Instead, the effect of the chemi- cal reactions was simulated by modifying the density of the continuous phase so that it corresponded to the mean temperature and composition of the flue gases, In this stage, the particle transport was simu- lated using the standard Euler-Euler and novel hybrid Euler-Lagrange approaches, The obtained results were compared against measured data, and both models were compared to each other. In the second part, the numerical model was enhanced by including the chemistry and physics of combustion. To the best of the authors＇ knowledge, the use of the hybrid Euler-Lagrange approach to model combustion is a new engineering application of this model, In this work, the combustion process was modeled for air-fuel combustion. The simulation results were compared with experimental data.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics (CFD) modeling of the gas–solids two-phase flow in a circulating fluidized bed (CFB) riser is carried out. The Eularian–Eularian method with the kinetic theory of granular flow is used to solve the gas–solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory (Johnson & Jackson, 1987) with specularity coefficient and particle–wall restitution coefficient. The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operating conditions of the CFB riser.     

含液颗粒材料液固耦合分析的离散颗粒模型及特征线SPH法

对含液颗粒材料流固耦合分析建议了一个基于离散颗粒模型与特征线SPH法的显式拉格朗日-欧拉无网格方案。在已有的用以模拟固体颗粒集合体的离散颗粒模型[1]基础上,将颗粒间间隙内的流体模型化为连续介质,对其提出并推导了基于特征线的SPH法。数值例题显示了所建议方案在模拟颗粒材料与间隙流相互作用的能力和性能以及间隙流体对颗粒结构承载能力及变形的影响。     

颗粒流动力学及其离散模型评述

颗粒流是由众多颗粒组成的具有内在相互作用的非经典介质流动. 自然界常见颗粒流都是密集流, 颗粒间接触形成力链, 诸多力链相互交接构成支撑整个颗粒流重量和外载荷的网络, 其局部构型及强度在外载荷下演化, 是颗粒流摩擦特性和接触应力的来源.本文介绍球形颗粒间无粘连作用时的Hertz法向接触理论和Mindlin-Deresiewicz切向接触理论. Campbell依据是否生成较为稳定的力链把颗粒流分为弹性流和惯性流两大类, 其中弹性-准静态流和惯性-碰撞流分别对应准静态流和快速流, 作为两种极端流动情况通常处理成连续体, 分别采用摩擦塑性模型和动理论予以描述, 但是表征接触力链的颗粒弹性参数并不出现这两个模型和理论框架中, 如何进一步考虑颗粒弹性参数将非常困难. 目前离散动力学方法逐渐成为复现其复杂颗粒流动现象、提取实验不可能获得的内部流动信息进而综合起来探索颗粒流问题的一种有效工具, 其真实性强于连续介质理论的描述. 软球模型对颗粒间接触力简化处理, 忽略了切向接触力对法向接触力及其加载历史的依赖, 带来了法向和切向刚度系数如何标度等更艰难的物理问题, 但由于计算强度小而广泛应用于工程问题中. 硬球模型不考虑颗粒接触变形, 因而不能描述颗粒流内在接触应变等物理机理, 仅适用于快速颗粒流, 这不仅仅是由于两体碰撞的限制. 因此基于颗粒接触力学的离散颗粒动力学模型是崭新的模型,适用于准静态流到快速流整个颗粒流态的模拟, 可以细致考虑接触形变及接触力的细节,建立更为合理的颗粒流本构关系, 进而有力的促进颗粒流这一非经典介质流动的研究.      

Influence of drag laws on pressure and bed material recirculation rate in a cold flow model of an 8 MW dual fluidized bed system by means of CPFD

A cold flow model of an 8 MW dual fluidized bed (DFB) system is simulated using the commercial computational particle fluid dynamics (CPFD) software package Barracuda. The DFB system comprises a bubbling bed connected to a fast fluidized bed with the bed material circulating between them. As the hydrodynamics in hot DFB plants are complex because of high temperatures and many chemical reaction processes, cold flow models are used. Performing numerical simulations of cold flows enables a focus on the hydrodynamics as the chemistry and heat and mass transfer processes can be put aside. The drag law has a major influence on the hydrodynamics, and therefore its influence on pressure, particle distribution, and bed material recirculation rate is calculated using Barracuda and its results are compared with experimental results. The drag laws used were energy-minimization multiscale (EMMS), Ganser, Turton–Levenspiel, and a combination of Wen–Yu/Ergun. Eleven operating points were chosen for that study and each was calculated with the aforementioned drag laws. The EMMS drag law best predicted the pressure and distribution of the bed material in the different parts of the DFB system. For predicting the bed material recirculation rate, the Ganser drag law showed the best results. However, the drag laws often were not able to predict the experimentally found trends of the bed material recirculation rate. Indeed, the drag law significantly influences the hydrodynamic outcomes in a DFB system and must be chosen carefully to obtain meaningful simulation results. More research may enable recommendations as to which drag law is useful in simulations of a DFB system with CPFD.     

Numerical modeling of fluid-particle interaction in granular media

Fluid-particle interaction underpins important behavior of granular media. Particle-scale simulation may help to provide key microscopic information governing the interaction and offer better understanding of granular media as a whole. This paper presents a coupled computational fluid dynamics and discrete element method (CFD-DEM) approach for this purpose. The granular particle system is modeled by DEM, while the fluid flow is simulated by solving the locally averaged Navier–Stokes equation with CFD. The coupling is considered by exchanging such interaction forces as drag force and buoyancy force between the DEM and CFD. The approach is benchmarked by two classic geomechanics problems for which analytical solutions are available, and is further applied to the prediction of sand heap formation in water through hopper flow. It is demonstrated that the key characteristic of granular materials interacting with pore water can be successfully captured by the proposed method.     

A three-dimensional yield-criterion-based flow model for avalanches

In this paper, a flow model for avalanches based on the three-dimensional yield criterion is presented in an attempt to allow the relaxation of the assumption of lateral confinement pressure that is adopted in the traditional three-dimensional Savage–Hutter model (S–H model). One of the advantages of this model is that a simplified constitutive relationship for granular flow, which could reveal the internal mechanism of avalanches, is adopted. Additionally, another advantage is that the strength parameters used in the proposed model are readily available for natural materials. The flow properties of avalanches are influenced by the generalized friction coefficient, which is a parameter that can be assessed by introducing the three-dimensional yield criterion. By comparing the results obtained by numerical simulations using the model proposed in this paper and laboratory experiments, a reasonably good agreement can be reached with regard to the prediction of the moving process of avalanches.     

Numerical simulation of the interaction between flow and flexible nets

A numerical approach is proposed to simulate the interaction between flow and flexible nets in steady current. The numerical approach is based on the joint use of the porous-media model and the lumped-mass model. The configuration of flexible nets can be simulated using the lumped-mass model and the flow field around fishing nets can be simulated using the porous-media model. Using an appropriate iterative scheme, the fluid–structure interaction problem can be solved and the steady flow field around flexible nets can be obtained. In order to validate the numerical models, the numerical results were compared with the data obtained from corresponding physical model tests. The comparisons show that the numerical results are in good agreement with the experimental data. Using the proposed numerical approach, this paper presents the flow field around a single flexible net and two flexible nets with a spacing distance. Both the configuration of the flexible nets and the flow velocity results are in accordance with those of the corresponding physical model tests.     

3D full-loop simulation and experimental verification of gas-solid flow hydrodynamics in a dense circulating fluidized bed

Because of their advantages of high efficiency and low cost, numerical research methods for large-scale circulating fluidized bed （CFB） apparatus are gaining ever more importance. This article presents a numer- ical study of gas-solid flow dynamics using the Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-scale （EMMS） model. A three-dimensional, full-loop, time-dependent simulation of the hydrodynamics of a dense CFB apparatus is performed. The process parameters （e.g., operating and initial conditions） are provided in accordance with the real experiment to enhance the accuracy of the simulation. The axial profiles of the averaged solid volume fractions and the solids flux at the outlet of the cyclone are in reasonable agreement with experimental data, thereby verifying the applicability of the mathematical and physical models. As a result, the streamline in the riser and standpipe as well as the solids distribution contours at the cross sections is analyzed. Computational fluid dynamics （CFD） serves as a basis for CFB modeling to help resolve certain issues long in dispute but difficult to address experimentally. The results of this study provide the basis of a general approach to describing dynamic simulations of gas-solid flows.     

Modeling Variable Density Flow and Solute Transport in Porous Medium: 2. Re‐Evaluation of the Salt Dome Flow Problem

Case 5, Level 1 of the international HYDROCOIN groundwater flow modeling project is an example of idealized flow over a salt dome. The groundwater flow is strongly coupled to solute transport since density variations in this example are large (20%).Several independent teams simulated this problem using different models. Results obtained by different codes can be contradictory. We develop a new numerical model based on the mixed hybrid finite elements approximation for flow, which provides a good approximation of the velocity, and the discontinuous finite elements approximation to solve the advection equation, which gives a good approximation of concentration even when the dispersion tensor is very small. We use the new numerical model to simulate the salt dome flow problem.In this paper we study the effect of molecular diffusion and we compare linear and nonlinear dispersion equations. We show the importance of the discretization of the boundary condition on the extent of recirculation and the final salt distribution. We study also the salt dome flow problem with a more realistic dispersion (very small dispersion tensor). Our results are different to prior works with regard to the magnitude of recirculation and the final concentration distribution. In all cases, we obtain recirculation in the lower part of the domain, even for only dispersive fluxes at the boundary. When the dispersion tensor becomes very small, the magnitude of recirculation is small. Swept forward displacement could be reproduced by using finite difference method to compute the dispersive fluxes instead of mixed hybrid finite elements.     

Fractional order spring/spring-pot/actuator element in a multibody system: Application of an expansion formula

Fractional order models of a spring/spring-pot and spring/spring-pot/actuator element connected into a multibody system are proposed in order to represent smart materials and components in adaptronic systems by introducing new tuning parameter. The models are introduced into dynamic equations via generalized forces and using the Lagrange's equations of the second kind in covariant form. Generalized forces are derived by taking into account fractional order derivatives in force–displacement relations and by using the principle of virtual work. The numerical scheme for solving fractional order differential equations proposed in Atanacković and Stanković (2008) is used in order to approximate fractional order derivative of a composite function appearing in the presented fractional order model. Numerical example for the multibody system with three degrees of freedom is presented. The results obtained for generalized forces are compared for different values of parameters in the fractional order derivative model.     

Numerical modeling of downslope flows of different rheology

Naturally occurring flow along a long homogeneous slope is numerically simulated. It is taken into account that the flow is able to capture the slope material and to entrain it into motion. The flow depth and velocity increase with time at the expense of the capture. The medium in motion is simulated using different rheological models including those of Herschel & Bulkley and Cross, as well as the power-law fluid model. For all the models the time dependences of the total depth and the mean flow velocity are obtained. The slope inclination effect on the dynamic flow parameters is studied. For the Herschel–Bulkley model the yield strength effect is also investigated. On the basis of the numerical calculations some assumptions are made and then used to derive asymptotic formulas for the bottom material entrainment rate at large times from the entrainment onset for all the above-listed rheological models.     

Mechanism of the Impact of Particle Size Distribution to Bed-Inventory Overturn for Pant-Leg Circulating Fluidized Bed

An improved two particle sizes numerical model based on a uniform size model was established to investigate the influence of the average particle size on bed-inventory overturn inside a pant-leg circulating fluidized bed (CFB). The new model successfully simulated the dynamic performance of a pant-leg CFB that as average particle size shrank, the pant-leg CFB tended to overturn, while the uniform size model showed a contradicted trend. The success was attributed to the difference of the flow pattern with different particle size. The smaller particles tend to stay the upper furnace after fluidized by the primary air flow while the larger particles tend to fall back to the bottom soon after being carried to upper furnace by the primary air flow and smaller particles. As pointed out in our previous work, the lateral mass transfer resulted in a lateral pressure difference at the upper finance and inhibited further lateral mass transfer, which was regarded as a self-balancing process. The quick fall down of the large particles somehow weaken the lateral pressure different built-up at the upper furnace. Therefore, as the average particle size shrink, the weaken effect of the large particles on self-balancing ability of the pant-leg CFB increase, resulting a more tendency for bed-inventory overturn. It was such a characteristic behavior of the large particles, which was neglected in the uniform particle size model, caused the difference between the results of the two models described above.     

Evaluation of the effect of wall boundary conditions on numerical simulations of circulating fluidized beds

A computational fluid dynamics （CFD） modeling of the gas-solids two-phase flow in a circulating fluidized bed （CFB） riser is carried out. The Eularian-Eularian method with the kinetic theory of granular flow is used to solve the gas-solids two-phase flow in the CFB riser. The wall boundary condition of the riser is defined based on the Johnson and Jackson wall boundary theory （Johnson ＆ Jackson, 1987） with specularity coefficient and particle-wall restitution coefficient.The numerical results show that these two coefficients in the wall boundary condition play a major role in the predicted solids lateral velocity, which affects the solid particle distribution in the CFB riser. And the effect of each of the two coefficients on the solids distribution also depends on the other one. The generality of the CFD model is further validated under different operatin~ conditions of the CFB riser.     

On the rheology of dilative granular media: Bridging solid- and fluid-like behavior

A new rate-dependent plasticity model for dilative granular media is presented, aiming to bridge the seemingly disparate solid- and fluid-like behavioral regimes. Up to date, solid-like behavior is typically tackled with rate-independent plasticity models emanating from Mohr–Coulomb and Critical State plasticity theory. On the other hand, the fluid-like behavior of granular media is typically treated using constitutive theories amenable to viscous flow, e.g., Bingham fluid. In our proposed model, the material strength is composed of a dilation part and a rate-dependent residual strength. The dilatancy strength plays a key role during solid-like behavior but vanishes in the fluid-like regime. The residual strength, which in a classical plasticity model is considered constant and rate-independent, is postulated to evolve with strain rate. The main appeal of the model is its simplicity and its ability to reconcile the classic plasticity and rheology camps. The applicability and capability of the model are demonstrated by numerical simulation of granular flow problems, as well as a classical shear banding problem, where the performance of the continuum model is compared to discrete particle simulations and physical experiment. These results shed much-needed light onto the mechanics and physics of granular media at various shear rates.     

SPH-based simulation of granular collapse on an inclined bed

The paper presents a numerical model for simulating a granular flow and its deposition on an inclined bed. A granular material is described as an elastic–plastic continuum and its constitutive law, namely Hooke's law, is discretized on the basis of the Smoothed Particle Hydrodynamics (SPH) method. In the equation of motion, however, the artificial viscosity, which is widely used in SPH, is not applied. The diffusive term derived from Hooke's law is introduced with a diffusion coefficient that varies depending on the stress and strain rate based on the Drucker–Prager yield function. The model is verified and validated through two numerical tests. It is shown that the basic elastic–perfectly plastic characteristics are reproduced with a simple shearing test. The effects of the diffusion coefficient and spatial resolution are investigated to show the validity of the model. In the simulation of the gravitational collapse of a granular column on an inclined bed, the performance of the model from the final deposition profile, the time history of the front position of the granular flow, the maximum runout distance, and the velocity profile are investigated for several cases of basal inclinations. The calculated results show good agreement with the experimental results.     

Numerical investigation of a perturbed swirling annular two-phase jet

A swirling annular gas–liquid two-phase jet flow system has been investigated by solving the compressible, time-dependent, non-dimensional Navier–Stokes equations using highly accurate numerical methods. The mathematical formulation for the flow system is based on an Eulerian approach with mixed-fluid treatment while an adjusted volume of fluid method is utilised to account for the gas compressibility. Surface tension effects are captured by a continuum surface force model. Swirling motion is applied at the inlet while a small helical perturbation is also applied to initiate the instability. Three-dimensional spatial direct numerical simulation has been performed with parallelisation of the code based on domain decomposition. The results show that the flow is characterised by a geometrical recirculation zone adjacent to the nozzle exit and by a central recirculation zone further downstream. Swirl enhances the flow instability and vorticity and promotes liquid dispersion in the cross-streamwise directions. A dynamic precessing vortex core is developed demonstrating that the growth of such a vortex in annular configurations can be initiated even at low swirl numbers, in agreement with experimental findings. Analysis of the averaged results revealed the existence of a geometrical recirculation zone and a swirl induced central recirculation zone in the flow field.     

颗粒流体系统的格子Boltzmann数值方法研究进展

介绍Euler-Lagrange框架下基于格子Boltzmann方法LBM (Lattice Boltzmann Method)发展的两种不同层次(即不同时-空尺度和精度)的颗粒流体系统离散模拟方法,即格子Boltzmann颗粒解析直接数值模拟(LB-based PR-DNS)方法和格子Boltzmann离散颗粒模拟(LB-based DPS)方法,总结了Euler-Euler框架下基于格子Boltzmann双流体模型(LB-based TFM)方面的探索研究。LB-based PR-DNS方法中颗粒尺寸远大于格子步长,能够直接解析出流体在颗粒表面的流动以及颗粒所受完整的动力学信息;LB-based DPS方法中格子步长远大于颗粒直径,其在计算精度、时间耗费和计算效率之间能达到很好的平衡,可获得流体的宏观平均流动及颗粒的运动轨迹信息。LB-based DNS和DPS是探索颗粒流体系统的有力手段,但LB-based TFM应用于模拟颗粒流体系统仍需进一步探索。     

Investigation on effect of drag models on flow behavior of power-law fluid–solid two-phase flow in fluidized bed

In this study, a Eulerian-Eulerian two-fluid model combined with the kinetic theory of granular flow is adopted to simulate power-law fluid–solid two-phase flow in the fluidized bed. Two new power-law liquid–solid drag models are proposed based on the rheological equation of power-law fluid and pressure drop. One called model A is a modified drag model considering tortuosity of flow channel and ratio of the throat to pore, and the other called model B is a blending drag model combining drag coefficients of high and low particle concentrations. Predictions are compared with experimental data measured by Lali et al., where the computed porosities from model B are closer to the measured data than other models. Furthermore, the predicted pressure drop rises as liquid velocity increases, while it decreases with the increase of particle size. Simulation results indicate that the increases of consistency coefficient and flow behavior index lead to the decrease of drag coefficient, and particle concentration, granular temperature, granular pressure, and granular viscosity go down accordingly.     

Numerical evaluation of turbulence models for dense to dilute gas–solid flows in vertical conveyor

A two-fluid model (TFM) of multiphase flows based on the kinetic theory and small frictional limit boundary condition of granular flow was used to study the behavior of dense to dilute gas–solid flows in vertical pneumatic conveyor. An axisymmetric 2-dimensional, vertical pipe with 5.6 m length and 0.01 m internal diameter was chosen as the computation domain, same to that used for experimentation in the literature. The chosen particles are spherical, of diameter 1.91 mm and density 2500 kg/m³. Turbulence interaction between the gas and particle phases was investigated by Simonin's and Ahmadi's models and their numerical results were validated for dilute to dense conveying of particles. Flow regimes transition and pressure drop were predicted. Voidage and velocity profiles of each phase were calculated in radial direction at different lengths of the conveying pipe. It was found that the voidage has a minimum, and gas and solid velocities have maximum values along the center line of the conveying pipe and pressure drop has a minimum value in transition from dense slugging to dilute stable flow regime. Slug length and pressure fluctuation reduction were predicted with increasing gas velocity, too. It is shown that solid phase turbulence plays a significant role in numerical prediction of hydrodynamics of conveyor and the capability of particles turbulence models depends on tuning parameters of slip-wall boundary condition.