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.     

Validation of a numerical model for the simulation of an electrostatic powder coating process

Numerical modeling of a complete powder coating process is carried out to understand the gas-particle two-phase flow field inside a powder coating booth and results of the numerical simulations are compared with experimental data to validate the numerical results. The flow inside the coating booth is modeled as a three-dimensional turbulent continuous gas flow with solid powder particles as a discrete phase. The continuous gas flow is predicted by solving Navier–Stokes equations using a standard k–ε turbulence model with non-equilibrium wall functions. The discrete phase is modeled based on a Lagrangian approach. In the calculation of particle propagation, a particle size distribution obtained through experiments is applied. The electrostatic field, including the effect of space charge due to free ions, is calculated with the use of the user defined scalar transport equations and user defined scalar functions in the software package, FLUENT, for the electrostatic potential and charge density.     

Extending EMMS-based models to CFB boiler applications

Recently, EMMS-based models are being widely applied in simulations of high-throughput circulating fluidized beds (CFBs) with fine particles. Its use for low flux systems, such as CFB boiler (CFBB), still remains unexplored. In this work, it has been found that the original definition of cluster diameter in EMMS model is unsuitable for simulations of the CFB boiler with low solids flux. To remedy this, we propose a new model of cluster diameter. The EMMS-based drag model (EMMS/matrix model) with this revised cluster definition is validated through the computational fluid dynamics (CFD) simulation of a CFB boiler.     

Virtual experimentation through 3D full-loop simulation of a circulating fluidized bed

Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-male (EMMS) model was used to simulate a semi-industry scale circulating fiuidized bed (CFB).Three-dimensional(3D), time-dependent simulation of a full-loop CFB revealed that the axial profiles of cross-sectionally averaged solid volume fraction,and the radial profiles of solid axial velocity and solid volume fraction were in reasonable agreement with experimental data.Based on this agreement,database derived from experiments not yet accomplished was replenished with such simulations, and fluid regime diagrams and pressure balance around the CFB loop were derived accordingly. This work presents an integrated viewpoint on CFB and unfolds a fresh paradigm fur CFB modeling, which can be expected to help resolve certain issues long in dispute but hard for experiments.     

Virtual experimentation through 3D full-loop simulation of a circulating fluidized bed

Eulerian granular multiphase model with a drag coefficient correction based on the energy-minimization multi-scale （EMMS） model was used to simulate a semi-industry scale circulating fluidized bed （CFB）. Three-dimensional （3D）, time-dependent simulation of a full-loop CFB revealed that the axial profiles of cross-sectionally averaged solid volume fraction, and the radial profiles of solid axial velocity and solid volume fraction were in reasonable agreement with experimental data. Based on this agreement, database derived from experiments not yet accomplished was replenished with such simulations, and fluid regime diagrams and pressure balance around the CFB loop were derived accordingly. This work presents an integrated viewpoint on CFB and unfolds a fresh paradigm for CFB modeling, which can be expected to help resolve certain issues long in dispute but bard for experiments.     

Simulation of mass balance behavior in a large-scale circulating fluidized bed reactor

We determine using a compound model the influence of the mass of granular matter on the behavior of a supercritical circulating fluidized bed (CFB) reactor. Population balance enables a stationary-regime modeling of the mass flow of granular matter inside a CFB unit in a large-scale. The simulation includes some important dynamic processes of gas-particle flows in fluidized bed such as attrition, fragmentation, elutriation, and fuel combustion. Numerical calculations with full boiler loading were performed of operational parameters such as furnace temperature, furnace pressure, feeding materials mass flows, and excess air ratio. Furthermore, three bed inventory masses were adopted as experimental variables in the simulation model of mass balance. This approach enables a sensitivity study of mass flows of granular matter inside a CFB facility. Some computational results from this population balance model obtained for a supercritical CFB reactor are presented that show consistency with the operational data for large-scale CFB units.     

Experimental and numerical study of pseudo-2D circulating fluidized beds

We present experimental investigations and numerical simulations of a pseudo-2D riser. Experiments were performed for various airflow rates, particle types/diameters, and particle size distributions. Pressure distributions along the wall of the riser were measured. Additional measurements from a smaller pseudo-2D riser (Kallio et al., 2009; Shah et al., 2012) were used to analyze horizontal solids volume fraction profiles. The experimental data were compared with simulation results carried out using an Euler–Euler approach. A mesh sensitivity study was conducted for numerical simulations and effects associated with simplifying real 3D geometry to a 2D model were examined. In addition, the effect of using an algebraic equation to represent the granular temperature versus a full partial differential equation also was examined for numerical simulations. Results showed small but significant near-wall sensitivity of the flow variables to mesh size. Substantial differences in mean pressure, solids distribution, and solid velocities were obtained, when 2D and 3D simulation results were compared. Finally, applying the simplified granular temperature equation for turbulent fluidization and for dilute-phase transport can lead to incorrect predictions in models.     

Modeling and simulation of circulating fluidized bed reactors applied to a carbonation/calcination loop

A fluid dynamic model for a gas-solid circulating fluidized bed （CFB） designed using two coupled riser reactors is developed and implemented numerically with code programmed in Matlab. The fluid dynamic model contains heat and species mass balances to calculate temperatures and compositions for a carbonation/calcination loop process. Because of the high computational costs required to resolve the three-dimensional phenomena, a model representing a trade-offbetween computational time requirements and accuracy is developed. For dynamic processes with a solid flux between the two reactor units that depends on the fluid dynamics of both risers, a dynamic one-dimensional two-fluid model is sufficient. A two-fluid model using the constant particle viscosity closure for the stress term is used for the solid phase, and an algebraic turbulence model is applied to the gas phase. The numerical model implementa- tion is based on the finite volume method with a staggered grid scheme. The exchange of solids between the reactor units constituting the circulating fluidized bed （solid flux） is implemented through additional mass source/sink terms in the continuity equations of the two phases, For model validation, a relevant experimental analysis provided in the literature is reproduced by the numerical simulations, The numerical analysis indicates that sufficient heat integration between the two reactor units is important for the performance of the circulating fluidized bed system, The two-fluid model performs fairly well for this chemical process operated in a CFB designed as two coupled riser reactors. Further analysis and optimization of the solution algorithms and the reactor coupling strategy is warranted.     

Modeling and simulation of circulating fluidized bed reactors applied to a carbonation/calcination loop

A fluid dynamic model for a gas-solid circulating fluidized bed (CFB) designed using two coupled riser reactors is developed and implemented numerically with code programmed in Matlab. The fluid dynamic model contains heat and species mass balances to calculate temperatures and compositions for a carbonation/calcination loop process.Because of the high computational costs required to resolve the three-dimensional phenomena, a model representing a trade-off between computational time requirements and accuracy is developed. For dynamic processes with a solid flux between the two reactor units that depends on the fluid dynamics of both risers, a dynamic one-dimensional two-fluid model is sufficient.A two-fluid model using the constant particle viscosity closure for the stress term is used for the solid phase, and an algebraic turbulence model is applied to the gas phase. The numerical model implementation is based on the finite volume method with a staggered grid scheme. The exchange of solids between the reactor units constituting the circulating fluidized bed (solid flux) is implemented through additional mass source/sink terms in the continuity equations of the two phases.For model validation, a relevant experimental analysis provided in the literature is reproduced by the numerical simulations. The numerical analysis indicates that sufficient heat integration between the two reactor units is important for the performance of the circulating fluidized bed system.The two-fluid model performs fairly well for this chemical process operated in a CFB designed as two coupled riser reactors. Further analysis and optimization of the solution algorithms and the reactor coupling strategy is warranted.     

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.     

稠密气固两相湍流流动的实验和数值模拟

基于气固两相流动模型计算循环流化床内稠密气固两相流湍流动,颗粒动理学方法模拟颗粒相湍动能,ＳＧＳ模型模拟气相湍流,采用γ－射线密度计和非等速取样管测量局部颗粒浓度和流率,利用ＦＦＴ方法计算颗粒浓度功率谱密度。模拟计算得到上升管内气相和固相速度和浓度分布等。同时数值模拟与Ｔｓｕｊｉ等和Ｋｎｏｗｌｔｏｎ等试验结果进行了比较,结果表明数值模拟计算与实验结果相吻合。     

A state-of-the-art review on transition between the fast fluidization and pneumatic transport of Geldart group B particles

The fluidization state in the circulating fluidized bed (CFB) boiler is crucial to its stable and safe operation. However, up to now, the research field has not reached unanimity on whether the fluidization regime that the upper furnace of the boiler operates in is the fast fluidization or pneumatic transport. To this end, this paper reviewed relevant research on the transition between the fast fluidization and pneumatic transport of Geldart group B particles, including the flow characteristics of the fast fluidization, the transition condition between the fast fluidization and pneumatic transport, the determination methods of the transport velocity u_tr and saturation carrying capacity G_s1 and the influencing factors on these two parameters. Previous research findings can provide certain guidelines for the design and optimization of the CFB boiler, and result in plenty of prediction correlations for u_tr and G_s1. Nonetheless, owing to insufficient data available on Geldart group B particles, especially the ones obtained under high temperature or pressure conditions and in large-scale CFB apparatuses, the existing correlations are not well suited for the prediction of u_tr and G_s1 of Geldart group B particles. Thus, further efforts are urgently demanded on the fast fluidization transition of Geldart group B particles.     

Anti-wear beam effects on gas–solid hydrodynamics in a circulating fluidized bed

Anti-wear beams installed on water walls of circulating fluidized bed (CFB) boilers are one of the most effective ways to protect against water-wall erosion. Beam effects from, for example, beam size and superficial gas velocity were investigated on gas–solid hydrodynamics in a CFB test rig using CFD simulations and experimental methods. The downward flow of the wall layer solids is observed to be disrupted by the beam but is then restored some distance further downstream. When falling solids from the wall layer hit the anti-wear beam, the velocity of the falling solids decreases rapidly. A fraction of the solids accumulates on the beam. Below the beams, the falling solids have reduced velocities but upward-moving solids were observed on the wall. The effect of the beam increases with width and superficial gas velocity. Wear occurs mainly above the beam and its variation with width is different above to below the beam. There is an optimum width that, when combined with beam height, results in less erosion.     

Modelling of solids distribution in stirred tanks: analysis of simulation strategies and comparison with experimental data

The predictive capabilities of CFD techniques as applied to solid–liquid stirred vessels are investigated. The distribution of solid particles was simulated in three baffled stirred tanks agitated with single and multiple impellers. Suspensions of glass beads of different diameters and average concentration up to 6?vol. % in water were studied. The simulations of solid–liquid suspensions in the stirred vessels were performed by using fully predictive approaches. Eulerian multiphase models were adopted for modelling the solid–liquid flow, coupled with three different extensions of the standard k-? model to the case of multiphase flows. The simulated particle axial concentration profiles are compared with experimental data and critically discussed. The most successful simulation strategy and one possible implementation are described.     

Particle size distribution in CPFD modeling of gas–solid flows in a CFB riser

A computational particle fluid dynamics (CPFD) numerical method to model gas–solid flows in a circulating fluidized bed (CFB) riser was used to assess the effects of particle size distribution (PSD) on solids distribution and flow. We investigated a binary PSD and a polydisperse PSD case. Our simulations were compared with measured solids concentrations and velocity profiles from experiments, as well as with a published Eulerian-Eulerian simulation. Overall flow patterns were similar for both simulation cases, as confirmed by experimental measurements. However, our fine-mesh CPFD simulations failed to predict a dense bottom region in the riser, as seen in other numerical studies. Above this bottom region, distributions of particle volume fraction and particle vertical velocity were consistent with our experiments, and the simulated average particle diameter decreased as a power function with riser height. Interactions between particles and walls also were successfully modeled, with accurate predictions for the lateral profiles of particle vertical velocity. It was easy to implement PSD into the CPFD numerical model, and it required fewer computational resources compared with other models, especially when particles with a polydisperse PSD were present in the heterogeneous flow.     

Experimental investigation on heat transfer and frictional characteristics of vertical upward rifled tube in supercritical CFB boiler

Water wall design is a key issue for supercritical Circulating Fluidized Bed (CFB) boiler. On account of the good heat transfer performance, rifled tube is applied in the water wall design of a 600 MW supercritical CFB boiler in China. In order to investigate the heat transfer and frictional characteristics of the rifled tube with vertical upward flow, an in-depth experiment was conducted in the range of pressure from 12 to 30 MPa, mass flux from 230 to 1200 kg/(m² s), and inner wall heat flux from 130 to 720 kW/m². The wall temperature distribution and pressure drop in the rifled tube were obtained in the experiment. The normal, enhanced and deteriorated heat transfer characteristics were also captured. In this paper, the effects of pressure, inner wall heat flux and mass flux on heat transfer characteristics are analyzed, the heat transfer mechanism and the frictional resistance performance are discussed, and the corresponding empirical correlations are presented. The experimental results show that the rifled tube can effectively prevent the occurrence of Departure from Nucleate Boiling (DNB) and keep the tube wall temperature in a permissible range under the operating condition of supercritical CFB boiler.     

CFD-DEM simulation of the hole cleaning process in a deviated well drilling: The effects of particle shape

We investigate the effect of particle shape on the transportation mechanism in well-drilling using a three-dimensional model that couples computational fluid dynamics (CFD) with the discrete element method (DEM). This numerical method allows us to incorporate the fluid–particle interactions (drag force, contact force, Saffman lift force, Magnus lift force, buoyancy force) using momentum exchange and the non-Newtonian behavior of the fluid. The interactions of particle−particle, particle−wall, and particle−drill pipe are taken into account with the Hertz–Mindlin model. We compare the transport of spheres with non-spherical particles (non-smooth sphere, disc, and cubic) constructed via the multi-sphere method for a range of fluid inlet velocities and drill pipe inclination angles. The simulations are carried out for laboratory-scale drilling configurations. Our results demonstrate good agreement with published experimental data. We evaluate the fluid–particle flow patterns, the particle velocities, and the particle concentration profiles. The results reveal that particle sphericity plays a major role in the fluid–solid interaction. The traditional assumption of an ideal spherical particle may cause inaccurate results.     

CFD-simulation of a circulating fluidized bed riser

In the current work, a model of the fluid mechanics in the riser of a circulating fluidized bed （CFB） has been implemented using computational fluid dynamics （CFD）. The model developed shall be used in future as the basis of 3D-reactor model for the simulation of large scale CFB combustors. The two-fluid model （TFM） approach is used to represent the fluid mechanics involved in the flow. The computational implementation is accomplished by the commercial software FLUENT. Different closure formulations are tested on a simplified geometry. Two different turbulence formulations, namely the swirl modified RNG k-e model and the Realizable k-e model, are tested in combination with two different approaches to solid phase turbulence, namely the dispersion and per phase approach. One focus of the current work is put on the study of different drag correlations. Besides the drag correlations by Syamlal et al. [Syamlal, M., Rogers, W., ＆ O＇Brien, T. J. （1993）. MFIX documentation theory guide. Technical Report DOE/METC-9411004, U.S. Department of Energy （DOE）. Morgantown Energy Technology Center： Morgantown, WV] and Gidaspow [Gidaspow, D. （1994）. Multiphaseflow andfluidization. New York： Academic Press] the EMMS model has been used to determine the momentum exchange between the two phases. The resulting formulation is then used to simulate a 1-m × 0.3-m cold CFB setup and is validated by experimental results [Schlichtharle, P. （2000）. Fluid dynamics and mixing of solids and gas in the bottom zone of circulating fluidized beds. Unoublished doctoral dissertation, Technische Universitaet Hamburg-Harburg, Shaker Verlag： Aachen].     

Two-dimension estimate of effective properties of solid with random voids

A discrete model for the estimation of the effect of random voids on the structural properties of a two-dimensional solid is presented. Triangular void cells are used to simulate randomly located micro-cracks.The proposed model is solved using Cell Method, a recent numerical method that allows a direct discrete formulation of balance equations. Both heterogeneities of the structure and stress concentrations due to voids distribution are taken into account by the model. Following an introduction, some aspects of Cell Method for plane elasticity that are relevant for this paper will be briefly recalled and results from simulations in the elastic field will be discussed. Next, the proposed model will be extended to plastic field and more simulations will be presented. Results show that this model can be successfully employed to assess the structural response of a 2D solid with randomly distributed voids.     

Experimental study and numerical simulation of the characteristics of a percussive gas–solid separator

The presence of solid particles in the flow of hypersonic wind tunnels damages the appearance of the experiment models in the wind tunnel and influences the accuracy of experimental results. The design of a highly efficient gas–solid separator was therefore undertaken. Particle trajectory imaging methods were used to measure trajectories under different conditions. The flow field and particle movement characteristics for different head angles (HAs) and separation tooth angles (STAs), inlet velocities, and the exhaust gas outlet pressures in the separator, were calculated using simulations based on the discrete phase model. The particle separation efficiency, pressure loss, and flow loss resulting from different structural parameters were also studied. In line with experimental observations, the characteristic angle of particle movements in the separator and the separation efficiency of the separator were found to increase with decreasing HA and with increasing STA. Separation efficiency improves with increasing inlet velocity and with increasing negative pressure of the exhaust gas outlet; however, the corresponding pressure loss and the flow rate of the waste gas also increased.