Numerical simulation and experimental verification of bubble size distribution in an air dense medium fluidized bed

Bubble size distribution is the basic apparent performance and obvious characteristics in the air dense medium fluidized bed (ADMFB). The approaches of numerical simulation and experimental verification were combined to conduct the further research on the bubble generation and movement behavior. The results show that ADMFB could display favorable expanded characteristics after steady fluidization. With different particle size distributions of magnetite powder as medium solids, we selected an appropriate prediction model for the mean bubble diameter in ADMFB. The comparison results indicate that the mean bubble diameters along the bed heights are 35 mm < D b < 66 mm and 40 mm < D b < 69 mm with the magnetite powder of 0.3 mm+0.15mm and 0.15mm+0.074mm, respectively. The prediction model provides good agreements with the experimental and simulation data. Based on the optimal operating gas velocity distribution, the mixture of magnetite powder and <1mm fine coal as medium solids were utilized to carry out the separation experiment on 6-50mm raw coal. The results show that an optimal separation density d P of 1.73g/cm 3 with a probable error E of 0.07g/cm 3 and a recovery efficiency of 99.97% is achieved, which indicates good separation performance by applying ADMFB.

A model for predicting bubble rise velocity in a pulsed gas solid fluidized bed

Bed stability, and especially the bed density distribution, is affected by the behavior of bubbles in a gas solid fluidized bed. Bubble rise velocity in a pulsed gas-solid fluidized bed was studied using photographic and computational fluid dynamics methods. The variation in bubble rise velocity was investigated as a function of the periodic pulsed air flow. A predictive model of bubble rise velocity was derived: ub=ψ(Ut+Up-Umf)+kp(gdb)(1/2). The software of Origin was used to fit the empirical coefficients to give ψ = 0.4807 and kp = 0.1305. Experimental verification of the simulations shows that the regular change in bubble rise velocity is accurately described by the model. The correlation coefficient was 0.9905 for the simulations and 0.9706 for the experiments.

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 finegrid,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(MS3 DATA)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.

Numerical simulation on the fluidized bed gasification and CaO dechlorination of refuse derived fuel

A three-dimensional numerical model verified by previous experimental data is developed to simulate the fluidized bed gasification of refuse derived fuel （RDF）. The CaO dechlorination model obtained by the thermal gravity analysis （TGA） is coupled to investigate the process of CaO dechlorination. An Eulerian-Eulerian method is adopted to simulate the gas-solid flow and self-developed chemical reaction modules are used to simulate chemical reactions. Flow patterns, gasification results and dechlorination efficiency are obtained by numerical simulation. Meanwhile, simulations are performed to evaluate the effects of Ca/Cl molar ratio and temperature on dechlorination efficiency. The simulation results show that the presence of bubbles in the gasifier lowers the CaO dechlorination efficiency. Increasing the Ca/Cl molar ratio can enhance the dechlorination efficiency. However, with the temperature increasing, the dechlorination efficiency increases initially and then decreases. The optimal Ca/Cl molar ratio is in the range of 3. 0 to 3. 5 and the optimal temperature is 923K.

Applicability of Markov chain-based stochastic model for bubbling fluidized beds

A Markov chain-based stochastic model （MCM） is developed to simulate the movement of particles in a 2D bubbling fluidized bed （BFB）. The state spaces are determined by the discretized physical cells of the bed, and the transition probability matrix is directly calculated by the results of a discrete element method （DEM） simulation. The Markov property of the BFB is discussed by the comparison results calculated from both static and dynamic transition probability matrices. The static matrix is calculated based on the Markov chain while the dynamic matrix is calculated based on the memory property of the particle movement. Results show that the difference in the trends of particle movement between the static and dynamic matrix calculation is very small. Besides, the particle mixing curves of the MCM and DEM have the same trend and similar numerical values, and the details show the time averaged characteristic of the MCM and also expose its shortcoming in describing the instantaneous particle dynamics in the BFB.

Characterization of Pressure Signals in Fluidized Beds Loaded with Large Particles Using Wigner Distribution Analysis: Feasibility of Diagnosis of Agglomeration

An experimental verification is reported on the early predicting index of agglomeration in bubbling fluidized bed. Coarse quartz sand, which has the same density but larger diameter than the bed material, was used to simulate the initial agglomerated particle. Wigner distribution was used to analyze the pressure fuctuation of the tested bed, and the average amplitude of local domain frequency （LDF） and local peak weighted average frequency （LPWA） under different operating conditions were measured and compared. The results showed that the LDF is sensitive to the agglomeration phenomena and had quick response to the incipient agglomeration in fluidized beds. It can be concluded from the results that these two parameters could be taken as the characteristic indexes to the agglomeration in fuidized beds.

An analysis approach of mass and energy balance in a dual-reactor circulating fluidized bed system

An analysis approach considering gas-solids hydrodynamics,reaction kinetics and reacting species nonuniformity together in a dual-reactor system is presented for better understanding its mass and energy balance.It was achieved by a 3-dimensional comprehensive hydrodynamics and reaction model for the dual-reactor system,which was developed from the successfully verified 3-dimensional comprehensive combustion model for one circulating fluidized bed(CFB)system(Xu and Cheng,2019).The developed model and analysis approach was successfully used on a 1 MW circulating fluidized bed–bubbling fluidized bed(CFB-BFB)dual-reactor system.Results showed the sensible and chemical energy between two reactors as well as the energy distributions in each reactor were balanced and they agreed well with the experimental measurements.The analysis approach indicated energy balance had a close relationship with the mass transfer in the CFB-BFB dual-reactor system.It may be applied in a design and operation optimization for a dual-reactor system.

Investigation into the operation of an autothermal two-section subbituminous coal fluidized bed gasifier

Using a newly developed experimental setup,the features and advantages of an autothermal single-casing atmospheric sub-bituminous coal fluidized bed air-blown gasifier,combining a combustion and gasification section,and mixing the dispersed phase(inert material,char)and heat exchange between them through an annular transfer device,have been revealed.To increase the efficiency of the gasifier,an experimental-computational method was developed find the conditions for optimal operation,combining changing the annular flow's geometry and regulating the primary air for gasification.A simple and reliable multizone thermodynamic calculation model makes it possible to predict the composition of char and syngas in the gasification section with acceptable accuracy.This method confirmed that a two-section fluidized bed gasifier can provide efficient gasification of solid fuels and is suitable for use in small-scale cogeneration plants.Syngas with a heating value of 3.6-4.5 MJ/m^(3)and CGE of 38.2%-42.3%was obtained in the experimental setup without optimizing the primary air flow rate.With optimization,the indicators increased to the heating value of syngas of 5.20-5.34 MJ/m^(3)and CGE of 42.5%-50.0%.With heat regeneration of 0.8,CGE increases to 70%.

Comparative analysis on gas–solid drag models in MFIX-DEM simulations of bubbling fluidized bed

In this study,the open-source software MFIX-DEM simulations of a bubbling fluidized bed(BFB)are applied to assess nine drag models according to experimental and direct numerical simulation(DNS)results.The influence of superficial gas velocity on gas–solid flow is also examined.The results show that according to the distribution of time-averaged particle axial velocity in y direction,except for Wen–Yu and Tenneti–Garg–Subramaniam(TGS),other drag models are consistent with the experimental and DNS results.For the TGS drag model,the layer-by-layer movement of particles is observed,which indicates the particle velocity is not correctly predicted.The time domain and frequency domain analysis results of pressure drop of each drag model are similar.It is recommended to use the drag model derived from DNS or fine grid computational fluid dynamics–discrete element method(CFD-DEM)data first for CFD-DEM simulations.For the investigated BFB,the superficial gas velocity less than 0.9 m·s^(-1) should be adopted to obtain normal hydrodynamics.

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 hydrodynamic disparity between homogeneous and heterogeneous fluidization is accounted for by the heterogeneity index,H_(d),which can be affected by choice of bubble diameter correlation.How this choice of bubble diameter correlation influences the model prediction calls for further fundamental research.This article incorporated seven different bubble diameter correlations into EMMS/bubbling drag model and studied their effects on H_(d).The performance of these correlations has been compared with the correlation used previously by EMMS/bubbling drag model.We found that some of the correlations predicted lower Hd by order of a magnitude than the correlation used by the original EMMS/bubbling drag.Based on such analysis,we proposed a modification in the EMMS drag model for bubbling and turbulent fluidized beds.A computational fluid dynamics(CFD)simulation using two-fluid model with the modified EMMS/bubbling drag model was performed for two bubbling and one turbulent fluidized beds.Voidage distribution,time averaged solid concentration and axial solid concentration profiles were studied and compared with the previous version of the EMMS/bubbling drag model and experimental data.We found that the right choice of bubble diameter correlations can significantly improve the results for CFD simulations.

Hydrodynamics and OCM reaction performance in a bubbling fluidized bed reactor and a riser reactor

The reasonable reactor design is of great importance for increasing the C_(2) yield(C2H4 and C2H6)of the oxidative coupling of methane(OCM),and the OCM reactor should remove the heat released in reactions quickly and efficiently and minimize the consecutive reaction of ethylene to carbon oxides.The fluidized bed reactor is characterized by excellent heat transfer,superior mass transport,and large handling capacity,while fewer studies focused on large-scale fluidized bed reactors for the OCM reaction.Therefore,large cold-model experiments and computational fluid dynamics simulations were conducted to investigate hydrodynamics and the OCM reaction performance in a large-scale bubbling fluidized bed(BFB)and a large-scale riser.In the BFB reactor,consecutive reactions of ethylene are acute because of the strong gas back-mixing,high solids holdup,and non-uniform solids distribution.While the consecutive reactions of ethylene are negligible due to the plug flow structure and low solids holdup in the riser reactor.Further,both reactors can achieve isothermal operation for the OCM process.The C_(2) selectivity of 45.4% and C_(2) yield of 21.1% are obtained in the riser reactor,increasing by 20.3% and 5.8% individually than that in the BFB reactor.This study provides useful information and reference to the OCM reactor designandcommercialization.

Scale-up of bubbling fluidized beds with continuous particle flow based on particle-residence-time distribution

Few studies have investigated scale-up of the residence-time distribution （RTD） of particles in bubbling fluidized beds （BFBs） with continuous particle flow. Two approaches were investigated in this study： first, using well-known scaling laws that require changes in particle properties and gas velocity; second, using a simple approach keeping the same particles and gas velocity for different beds. Our theoretical analysis indicates it is possible to obtain similar RTDs in different BFBs with scaling laws if the plug-flow residence time （tpiug） is changed as m^0.5, where m is the scaling ratio of the bed; however, neither approach can ensure similar RTDs if tplug is kept invariant. To investigate RTD variations using two approaches without changing tplug, we performed experiments in three BFBs. The derivatives dE（θ）/dθ （where E（θ） is the dimensionless RTD density function and θ is the dimensionless time） in the early stage of the RTDs always varied with m 1, which was attributed to the fact that the particle movement in the early stage were mainly subject to dispersion. Using the simple approach, we obtained similar RTDs by separately treating the RTDs in the early and post-stages. This approach guarantees RTD similarity and provides basic rules for designing BFBs.

PRESSURE FLUCTUATIONS IN GAS-SOLIDS FLUIDIZED BEDS

Pressure fluctuation data measured in a series of fluidized beds with diameters of 0.05, 0.1, 0.29, 0.60 and 1.56 m showed that the maximum amplitude or standard deviation increased with increasing the superficial gas velocity and static bed height for relatively shallow beds and became insensitive to the increase in static bed height in relatively deep beds. The amplitude appeared to be less dependent on the measurement location in the dense bed. Predictions based on bubble passage, bubble eruption at the upper bed surface and bed oscillation all failed to explain all observed trends and underestimated the amplitude of pressure fluctuations, suggesting that the global pressure fluctuations in gas-solids bubbling fluidized beds are the superposition of local pressure variations, bed oscillations and pressure waves generated from the bubble formation in the distributor region, bubble coalescence during their rise and bubble eruption at the upper bed surface.

High-density circulating fluidized bed gasifier for advanced IGCC/IGFC-Advantages and challenges

Coal-fired Integrated Gasification Combined Cycle （IGCC） and Integrated coal Gasification Fuel-cell Com- bined cycle （IGFC） are being developed as high-efficiency electric power generation technology. However, the highest theoretical gross thermal efficiency of the conventional IGCC]IGFC is still below 52~. In order to obtain higher power generation efficiency, an advanced IGCC （A-IGCC） or advanced IGFC （A-IGFC） sys- tem making use of the exergy recuperation concept by recycling waste heat from gas turbine or fuel cells for steam gasification of coal and biomass was proposed in our laboratory, Corresponding to this system, a novel high-density triple-bed combined circulating fluidized bed （TBCFB） gasifier, composed of a downer pyrolyzer, a bubbling fluidized bed char gasifier, and a riser combustor, was proposed to replace traditional gasifiers such as the entrained flow bed gasifier. The new system is expected to more effectively utilize the waste heat from gas turbines or fuel cells and the heat produced by the combustion of the unreacted char in the riser combustor for pyrolysis and gasification of coal and biomass. In this short review, the advantages and future challenges in the development of high-density TBCFB gasifier are presented and discussed.

Effect of bed size on the gas-solid flow characterized by pressure fluctuations in bubbling fluidized beds

Pressure fluctuations in four bubbling fluidized beds having different bed sizes (three square cross-sections of 5, 10, and 15 cm in side length, and one rectangular cross-section of 2 × 10 cm2) were measured at four axial positions (P1, P2, P3, and P4). Several characteristic indicators of the flow specifically of the pressure were calculated. In terms of these characteristic indicators, the effect of bed size on flow behavior was investigated. The results show that in the fully fluidized state, the pressure drop is slightly higher in smaller beds, but the pressure drops in the 10- and 15-cm beds are close. The 15-cm bed has the lowest pressure-fluctuation amplitude. The amplitudes at P1 and P2 in the lower part of the bed are very close for bed sizes below 10 cm, but the amplitude at P3 near the bed surface increases with decreasing bed size. No general trend was observed regarding the effect of bed size on skewness and kurtosis of the pressure for all four axial heights. For the average, standard deviation, skewness, and kurtosis of the pressure at P4, the values are close for the two small beds (2 × 10 and 5 × 5 cm2) and the two large beds (10 × 10 and 15 × 15 cm2), and hence the effect of bed size separates the beds into two groups. In the fully fluidized state, for P1, P2, and P3, the Kolmogorov entropy and the dominant frequency both increase with increasing bed size, but in the pseudo-2D bed both are between the values for the 5- and 10-cm beds.

CFD simulation for reduction of pyrolusite in fluidized beds

In this study,a CFD model coupled with heterogeneous flow structure,mass transfer equations,and chemical reaction kinetics is established to forecast the pyrolusite reduction reaction behavior.Compared with the previous studies which ignore the volume change of solids phase,the influence of volume shrinkage on reaction and flow behavior is explored in this research.Volume shrinkage of pyrolusite is proved to be non-negligible in predicting the conversion rate.The negligence of volume shrinkage leads to the overestimation of conversion rate for its inaccurate estimation of surface area for reaction.Besides,the influence of volume shrinkage on the reaction is found smaller in the scaled-up reactor.

A new drag model for TFM simulation of gas-solid bubbling fluidized beds with Geldart-B particles

In this work, a new drag model for TFM simulation in gas-solid bubbling fluidized beds was proposed, and a set of equations was derived to determine the meso-scale structural parameters to calculate the drag characteristics of Geldart-B particles under low gas velocities. In the new model, the meso-scale structure was characterized while accounting for the bubble and meso-scale structure effects on the drag coefficient. The Fluent software, incorporating the new drag model, was used to simulate the fluidization behavior. Experiments were performed in a Plexiglas cylindrical fluidized bed consisting of quartz sand as the solid phase and ambient air as the gas phase. Comparisons based on the solids hold-up inside the fluidized bed at different superficial gas velocities, were made between the 2D Cartesian simulations, and the experimental data, showing that the results of the new drag model reached much better agreement with exoerimental data than those of the Gidasoow dra~ model did.

An exploratory study of three-dimensional MP-PIC-based simulation of bubbling fluidized beds with and without baffles

In this study, the flow characteristics of Geldart A particles in a bobbling fluidized bed with and without perforated plates were simulated by the multiphase particle-in-cell （MP-PlC）-based Eolerian-Lagrangian method. A modified structure-based drag model was developed based on our previous work. Other drag models including the Parker and Wen-Yo-Ergon drag models were also employed to investigate the effects of drag models on the simulation results. Although the modified structure-based drag model better predicts the gas-solid flow dynamics of a baffle-free bubbling fluidized bed in comparison with the experimental data, none of these drag models predict the gas-solid flow in a baffled bobbling floidized bed sufficiently well because of the treatment of baffles in the Barracuda software. To improve the simulation accuracy, future versions of Barracuda should address the challenges of incorporating the bed height and the baffles.

Experimental study on reactor performance of gas-solids low-velocity fluidized beds

Reactor performance of bubbling fluidized bed(BFB)and turbulent fluidized bed(TFB)was carefully examined and systematically compared using catalytic ozone decomposition as a model reaction,based on a complete mapping of local flow structures and spatial distributions of ozone conversion and solids holdup.TFB clearly has a higher conversion and shows better reactor performance than BFB as a result of the vigorously turbulent flow and the relatively homogeneous gas–solids mixing in TEB.Besides,the intensive interaction between gas and solids in TFB leads to greater gas–solids contact efficiency of TFB over that of BFB.Due to gas bypassing and backmixing caused by bubbling behaviours and two-phase structure,BFB deviates significantly from a plug flow reactor and sometimes from a continuously stirred tank reactor.The flow structures essentially dictate the reactor performance in the low-velocity fluidized beds.

Effect of agitation on the characteristics of air dense medium fluidization

In order to study the effect of agitation on the characteristics of air dense medium fluidization, we designed and constructed an agitation device. Analyses were then conducted on the fluidization characteristics curves, the bed density stability and the average bubble rise velocity Uaunder different agitation conditions. The results indicated that a lower bed pressure drop(without considering lower gas velocity in a fixed bed stage) and higher minimum fluidized velocity are achieved with increasing agitation speed.The height d(distance between the lower blades and air distribution plate) at which the agitation paddle was located had a considerable effect on the stability of the bed density at 9.36 cm/s < U < 10.70 cm/s. The higher the value of d, the better the stability, and the standard deviation of the bed density fluctuation r dropped to 0.0364 g/cm^3 at the ideal condition of d = 40 mm. The agitation speed also had a significant influence on the fluidization performance, and r was only 0.0286 g/cm^3 at an agitation speed of N = 75 r/min. The average bubble rise velocity decreased significantly with increasing agitation speed under the operating condition of 1.50 cm/s < U–U_(mf)< 3.50 cm/s. This shows that appropriate agitation contributes to a significant improvement in the fluidization quality in a fluidized bed, and enhances the separation performance of a fluidized bed.