Modelling and simulation of biomass gasification in a circulating fluidized bed reactor

Computational Fluid Dynamics (CFD) has found great acceptance among the engineering community as a tool for research and design of processes that are practically difficult or expensive to study experimentally. One of these processes is the biomass gasification in a Circulating Fluidized Bed (CFB). Biomass gasification is the thermo-chemical conversion of biomass at a high temperature and a controlled oxygen amount into fuel gas, also sometime referred to as syngas. Circulating fluidized bed is a type of reactor in which it is possible to maintain a stable and continuous circulation of solids in a gas-solid system. The main objectives of this thesis are four folds: (i) Develop a three-dimensional predictive model of biomass gasification in a CFB riser using advanced Computational Fluid Dynamic (CFD) (ii) Experimentally validate the developed hydrodynamic model using conventional and advanced measuring techniques (iii) Study the complex hydrodynamics, heat transfer and reaction kinetics through modelling and simulation (iv) Study the CFB gasifier performance through parametric analysis and identify the optimum operating condition to maximize the product gas quality. Two different and complimentary experimental techniques were used to validate the hydrodynamic model, namely pressure measurement and particle tracking. The pressure measurement is a very common and widely used technique in fluidized bed studies, while, particle tracking using PEPT, which was originally developed for medical imaging, is a relatively new technique in the engineering field. It is relatively expensive and only available at few research centres around the world. This study started with a simple poly-dispersed single solid phase then moved to binary solid phases. The single solid phase was used for primary validations and eliminating unnecessary options and steps in building the hydrodynamic model. Then the outcomes from the primary validations were applied to the secondary validations of the binary mixture to avoid time consuming computations. Studies on binary solid mixture hydrodynamics is rarely reported in the literature. In this study the binary solid mixture was modelled and validated using experimental data from the both techniques mentioned above. Good agreement was achieved with the both techniques. According to the general gasification steps the developed model has been separated into three main gasification stages; drying, devolatilization and tar cracking, and partial combustion and gasification. The drying was modelled as a mass transfer from the solid phase to the gas phase. The devolatilization and tar cracking model consist of two steps; the devolatilization of the biomass which is used as a single reaction to generate the biomass gases from the volatile materials and tar cracking. The latter is also modelled as one reaction to generate gases with fixed mass fractions. The first reaction was classified as a heterogeneous reaction while the second reaction was classified as homogenous reaction. The partial combustion and gasification model consisted of carbon combustion reactions and carbon and gas phase reactions. The partial combustion considered was for C, CO, H2 and CH4. The carbon gasification reactions used in this study is the Boudouard reaction with CO2, the reaction with H2O and Methanation (Methane forming reaction) reaction to generate methane. The other gas phase reactions considered in this study are the water gas shift reaction, which is modelled as a reversible reaction and the methane steam reforming reaction. The developed gasification model was validated using different experimental data from the literature and for a wide range of operating conditions. Good agreement was observed, thus confirming the capability of the model in predicting biomass gasification in a CFB to a great accuracy. The developed model has been successfully used to carry out sensitivity and parametric analysis. The sensitivity analysis included: study of the effect of inclusion of various combustion reaction; and the effect of radiation in the gasification reaction. The developed model was also used to carry out parametric analysis by changing the following gasifier operating conditions: fuel/air ratio; biomass flow rates; sand (heat carrier) temperatures; sand flow rates; sand and biomass particle sizes; gasifying agent (pure air or pure steam); pyrolysis models used; steam/biomass ratio. Finally, based on these parametric and sensitivity analysis a final model was recommended for the simulation of biomass gasification in a CFB riser.

Computationsl modelling of dimethyl ether separation and steam reforming in fluidized bed reactors

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether (DME) gas adsorptive separation and steam reforming (DME-SR) in a large scale Circulating Fluidized Bed (CFB) reactor. The CFD model is based on Eulerian-Eulerian dispersed flow and solved using commercial software (ANSYS FLUENT). Hydrogen is currently receiving increasing interest as an alternative source of clean energy and has high potential applications, including the transportation sector and power generation. Computational fluid dynamic (CFD) modelling has attracted considerable recognition in the engineering sector consequently leading to using it as a tool for process design and optimisation in many industrial processes. In most cases, these processes are difficult or expensive to conduct in lab scale experiments. The CFD provides a cost effective methodology to gain detailed information up to the microscopic level. The main objectives in this project are to: (i) develop a predictive model using ANSYS FLUENT (CFD) commercial code to simulate the flow hydrodynamics, mass transfer, reactions and heat transfer in a large scale dual fluidized bed system for combined gas separation and steam reforming processes (ii) implement a suitable adsorption models in the CFD code, through a user defined function, to predict selective separation of a gas from a mixture (iii) develop a model for dimethyl ether steam reforming (DME-SR) to predict hydrogen production (iv) carry out detailed parametric analysis in order to establish ideal operating conditions for future industrial application. The project has originated from a real industrial case problem in collaboration with the industrial partner Dow Corning (UK) and jointly funded by the Engineering and Physical Research Council (UK) and Dow Corning. The research examined gas separation by adsorption in a bubbling bed, as part of a dual fluidized bed system. The adsorption process was simulated based on the kinetics derived from the experimental data produced as part of a separate PhD project completed under the same fund. The kinetic model was incorporated in FLUENT CFD tool as a pseudo-first order rate equation; some of the parameters for the pseudo-first order kinetics were obtained using MATLAB. The modelling of the DME adsorption in the designed bubbling bed was performed for the first time in this project and highlights the novelty in the investigations. The simulation results were analysed to provide understanding of the flow hydrodynamic, reactor design and optimum operating condition for efficient separation. Bubbling bed validation by estimation of bed expansion and the solid and gas distribution from simulation agreed well with trends seen in the literatures. Parametric analysis on the adsorption process demonstrated that increasing fluidizing velocity reduced adsorption of DME. This is as a result of reduction in the gas residence time which appears to have much effect compared to the solid residence time. The removal efficiency of DME from the bed was found to be more than 88%. Simulation of the DME-SR in FLUENT CFD was conducted using selected kinetics from literature and implemented in the model using an in-house developed user defined function. The validation of the kinetics was achieved by simulating a case to replicate an experimental study of a laboratory scale bubbling bed by Vicente et al [1]. Good agreement was achieved for the validation of the models, which was then applied in the DME-SR in the large scale riser section of the dual fluidized bed system. This is the first study to use the selected DME-SR kinetics in a circulating fluidized bed (CFB) system and for the geometry size proposed for the project. As a result, the simulation produced the first detailed data on the spatial variation and final gas product in such an industrial scale fluidized bed system. The simulation results provided insight in the flow hydrodynamic, reactor design and optimum operating condition. The solid and gas distribution in the CFB was observed to show good agreement with literatures. The parametric analysis showed that the increase in temperature and steam to DME molar ratio increased the production of hydrogen due to the increased DME conversions, whereas the increase in the space velocity has been found to have an adverse effect. Increasing temperature between 200 oC to 350 oC increased DME conversion from 47% to 99% while hydrogen yield increased substantially from 11% to 100%. The CO2 selectivity decreased from 100% to 91% due to the water gas shift reaction favouring CO at higher temperatures. The higher conversions observed as the temperature increased was reflected on the quantity of unreacted DME and methanol concentrations in the product gas, where both decreased to very low values of 0.27 mol% and 0.46 mol% respectively at 350 °C. Increasing the steam to DME molar ratio from 4 to 7.68 increased the DME conversion from 69% to 87%, while the hydrogen yield increased from 40% to 59%. The CO2 selectivity decreased from 100% to 97%. The decrease in the space velocity from 37104 ml/g/h to 15394 ml/g/h increased the DME conversion from 87% to 100% while increasing the hydrogen yield from 59% to 87%. The parametric analysis suggests an operating condition for maximum hydrogen yield is in the region of 300 oC temperatures and Steam/DME molar ratio of 5. The analysis of the industrial sponsor’s case for the given flow and composition of the gas to be treated suggests that 88% of DME can be adsorbed from the bubbling and consequently producing 224.4t/y of hydrogen in the riser section of the dual fluidized bed system. The process also produces 1458.4t/y of CO2 and 127.9t/y of CO as part of the product gas. The developed models and parametric analysis carried out in this study provided essential guideline for future design of DME-SR at industrial level and in particular this work has been of tremendous importance for the industrial collaborator in order to draw conclusions and plan for future potential implementation of the process at an industrial scale.

Fast pyrolysis of biomass in a circulating fluidised bed reactor

The objective of this work was to design, construct, test and operate a novel circulating fluid bed fast pyrolysis reactor system for production of liquids from biomass. The novelty lies in incorporating an integral char combustor to provide autothermal operation. A reactor design methodology was devised which correlated input parameters to process variables, namely temperature, heat transfer and gas/vapour residence time, for both the char combustor and biomass pyrolyser. From this methodology a CFB reactor was designed with integral char combustion for 10 kg/h biomass throughput. A full-scale cold model of the CFB unit was constructed and tested to derive suitable hydrodynamic relationships and performance constraints. Early difficulties encountered with poor solids circulation and inefficient product recovery were overcome by a series of modifications. A total of 11 runs in a pyrolysis mode were carried out with a maximum total liquids yield of 61.50% wt on a maf biomass basis, obtained at 500°C and with 0.46 s gas/vapour residence time. This could be improved by improved vapour recovery by direct quenching up to an anticipated 75 % wt on a moisture-and-ash-free biomass basis. The reactor provides a very high specific throughput of 1.12 - 1.48 kg/hm2 and the lowest gas-to-feed ratio of 1.3 - 1.9 kg gas/kg feed compared to other fast pyrolysis processes based on pneumatic reactors and has a good scale-up potential. These features should provide significant capital cost reduction. Results to date suggest that the process is limited by the extent of char combustion. Future work will address resizing of the char combustor to increase overall system capacity, improvement in solid separation and substantially better liquid recovery. Extended testing will provide better evaluation of steady state operation and provide data for process simulation and reactor modeling.

Fluidized bed gasification of biomass

A fluidized bed process development unit of 0.8 m internal diameter was designed on basis of results obtained from a bench scale laboratory unit. For the scaling up empirical models from the literature were used. The process development unit and peripheral equipment were constructed, assembled and commissioned, and instruments were provided for data acquisition. The fluidization characteristics of the reactor were determined and were compared to the design data. An experimental programme was then carried out and mass and energy balances were made for all the runs. The results showed that the most important independent experimental parameter was the air factor, with an optimum at 0.3. The optimum higher heating value of the gas produced was 6.5 MJ/Nm3, while the thermal efficiency was 70%. Reasonably good agreement was found between the experimental results, theoretical results from a thermodynamic model and data from the literature. It was found that the attainment of steady state was very sensitive to a continuous and constant feedstock flowrate, since the slightest variation in feed flow resulted in fluctuations of the gas quality. On the basis of the results a set of empirical relationships was developed, which constitutes an empirical model for the prediction of the performance of fluidized bed gasifiers. This empirical model was supplemented by a design procedure by which fluidized bed gasifiers can be designed and constructed. The design procedure was then extended to cover feedstock feeding and gas cleaning in a conceptual design of a fluidized bed gasification facility. The conceptual design was finally used to perform an economic evaluation of a proposed gasification facility. The economics of this plant (retrofit application) were favourable.

Mass transfer coefficient for drying of moist particulate in a bubbling fluidized bed

Experiments on drying of moist particles by ambient air were carried out to measure the mass transfer coefficient in a bubbling fluidized bed. Fine glass beads of mean diameter 125?µm were used as the bed material. Throughout the drying process, the dynamic material distribution was recorded by electrical capacitance tomography (ECT) and the exit air condition was recorded by a temperature/humidity probe. The ECT data were used to obtain qualitative and quantitative information on the bubble characteristics. The exit air moisture content was used to determine the water content in the bed. The measured overall mass transfer coefficient was in the range of 0.0145–0.021?m/s. A simple model based on the available correlations for bubble-cloud and cloud-dense interchange (two-region model) was used to predict the overall mass transfer coefficient. Comparison between the measured and predicted mass transfer coefficient have shown reasonable agreement. The results were also used to determine the relative importance of the two transfer regions.

Extended surface fluidized bed heat transfer

Experirnental data and theoretical calculation on the heat transfer performance of extended surface submerged: in shallow air fluidized beds ~ less than 150 mm, are presented. Energy t;ransferrence from the bed material was effected by water cooled tubes passing through the fins. The extended surface tested was either manufactured from square or radial copper fins silver soldered to a circular basic tube or commercially supplied, being of the crimped or extruded helical fin type. Performances are compared, for a wide range of geometric variables, bed configurations and fluidized materials, with plain and oval tubes operating under similar experimental conditions. A statistical analysis of all results, using a regression technique, has shown the relative importance of each significant variable. The bed to surface heat transfer coefficients are higher than those reported in earlier published work using finned tubes in much deeper beds and the heat transfer to the whole of the extended surface is at least as good as that previously reported for un-finned tubes. The improved performance is attributed partly to the absence of large bubbles in shallow beds and it is suggested that the improved circulation of the solids when constrained in the narrow passages between adjacent fins may be a contributory factor. Flow visualisation studies between a perspex extended surface and a fluidized bed using air at ambient temperatures, have demonstrated the effect of too small a fin spacing. Fin material and the bonding to the basic tube are more important in the optimisation of performance than in conventional convective applications because of the very much larger heat fluxes involved. A theoretical model of heat flow for a radial fin surface, provides data concerning the maximum heat transfer and minimum metal required to fulfil a given heat exchange duty. Results plotted in a series of charts aim at assisting the designer of shalJow fluidized beds.

Some studies relating to the design of shallow, fluidized bed boilers

The concept of shallow fluidized bed boilers is defined and a preliminary working design for a gas-fired package boiler has been produced. Those areas of the design requiring further study have been specified. Experimental investigations concerning these areas have been carried out. A two-dimensional, conducting paper analog has been developed for the specific purpose of evaluating sheet fins. The analog has been generalised and is presented as a simple means of simulating the general, two-dimensional Helmholtz equation. By recording the transient response of spherical, calorimetric probes when plunged into heated air-fluidized beds, heat transfer coefficients have been measured at bed temperatures up to 1 100°C. A correlation fitting all the data to within ±10% has been obtained. A model of heat transfer to surfaces immersed in high temperature beds has been proposed. The model solutions are, however, only in qualitative agreement with the experimental data. A simple experimental investigation has revealed that the effective, radial, thermal conductivities of shallow fluidized beds are an order of magnitude lower than the axial conductivities. These must, consequently, be taken into account when considering heat transfer to surfaces immersed within fluidized beds. Preliminary work on pre-mixed gas combustion and some further qualitative experiments have been used as the basis for discussing the feasibility of combusting heavy fuel oils within shallow beds. The use of binary beds, within which the fuel could be both gasified and subsequently burnt, is proposed. Finally, the consequences of the experimental studies on the initial design are considered, and suggestions for further work are made.

Fin spacing interaction in fluidized bed heat exchangers

The work presented in this thesis is concerned with the heat transfer performance of a single horizontal bare tube and a variety of finned tubes immersed in a shallow air fluidized bed. Results of experimental investigations with the bare tube indicate that the tube position in the bed influences its performance narticularly where fine bed materials are used. In some cases the maximum heat transfer is obtained with the tube in the particle cloud just above the dense phase fluidized bed - a phenomenon that has not been previously observed. This was attributed to the unusual particle circulation in shallow beds. The data is also presented in dimensionless correlations which may be useful for design purposes. A close approximation to the bare tube data can be obtained by using thetransient heating of a spherical robe and this provides a valuable way of accumulating a lot of data very rapidly. The experimental data on finned tubes shows that a fin spacing less than twenty times the average particle diameter can cause a significant reduction in heat transfer due to the interaction which takes place between the particles and the surface of the fins. Furthermore, evidence is provided to show that particle shape plays an important part in the interaction with spherical particles being superior to angular particles at low fin spacing/particle diameter ratio. The finned tube data is less sensitive to tube position in the bed than bare tubes and the best performance is when the tube is positioned at the distributor.A reduction in bed depth decreases the thermal performance of the finned tube but in many practical installations the reduction in pressure drop might more than comnensate for the reduced heat flux. Information is also provided on the theoretical uerformance of fins and the effect of the root contact area between the fins and the tube was investigated.

Time scale analysis for fluidized bed melt granulation I: granule-granule and granule-droplet collision rates

Fluidized bed spray granulators (FBMG) are widely used in the process industry for particle size growth; a desirable feature in many products, such as granulated food and medical tablets. In this paper, the first in a series of four discussing the rate of various microscopic events occurring in FBMG, theoretical analysis coupled with CFD simulations have been used to predict granule–granule and droplet–granule collision time scales. The granule–granule collision time scale was derived from principles of kinetic theory of granular flow (KTGF). For the droplet–granule collisions, two limiting models were derived; one is for the case of fast droplet velocity, where the granule velocity is considerable lower than that of the droplet (ballistic model) and another for the case where the droplet is traveling with a velocity similar to the velocity of the granules. The hydrodynamic parameters used in the solution of the above models were obtained from the CFD predictions for a typical spray fluidized bed system. The granule–granule collision rate within an identified spray zone was found to fall approximately within the range of 10-2–10-3 s, while the droplet–granule collision was found to be much faster, however, slowing rapidly (exponentially) when moving away from the spray nozzle tip. Such information, together with the time scale analysis of droplet solidification and spreading, discussed in part II and III of this study, are useful for probability analysis of the various event occurring during a granulation process, which then lead to be better qualitative and, in part IV, quantitative prediction of the aggregation rate.

Mass transfer in fluidized bed drying of moist particulate

Bubbling fluidized bed technology is one of the most effective mean for interaction between solid and gas flow, mainly due to its good mixing and high heat and mass transfer rate. It has been widely used at a commercial scale for drying of grains such as in pharmaceutical, fertilizers and food industries. When applied to drying of non-pours moist solid particles, the water is drawn-off driven by the difference in water concentration between the solid phase and the fluidizing gas. In most cases, the fluidizing gas or drying agent is air. Despite of the simplicity of its operation, the design of a bubbling fluidized bed dryer requires an understanding of the combined complexity in hydrodynamics and the mass transfer mechanism. On the other hand, reliable mass transfer coefficient equations are also required to satisfy the growing interest in mathematical modelling and simulation, for accurate prediction of the process kinetics. This chapter presents an overview of the various mechanisms contributing to particulate drying in a bubbling fluidized bed and the mass transfer coefficient corresponding to each mechanism. In addition, a case study on measuring the overall mass transfer coefficient is discussed. These measurements are then used for the validation of mass transfer coefficient correlations and for assessing the various assumptions used in developing these correlations.

A priori prediction of aggregation efficiency and rate constant for fluidized bed melt granulation

This paper presents a predictive aggregation rate model for spray fluidized bed melt granulation. The aggregation rate constant was derived from probability analysis of particle–droplet contact combined with time scale analysis of droplet solidification and granule–granule collision rates. The latter was obtained using the principles of kinetic theory of granular flow (KTGF). The predicted aggregation rate constants were validated by comparison with reported experimental data for a range of binder spray rate, binder droplet size and operating granulator temperature. The developed model is particularly useful for predicting particle size distributions and growth using population balance equations (PBEs).

Influence of reaction atmosphere (H2O, N2, H2, CO2, CO) on fluidized-bed fast pyrolysis of biomass using detailed tar vapor chemistry in computational fluid dynamics

Secondary pyrolysis in fluidized bed fast pyrolysis of biomass is the focus of this work. A novel computational fluid dynamics (CFD) model coupled with a comprehensive chemistry scheme (134 species and 4169 reactions, in CHEMKIN format) has been developed to investigate this complex phenomenon. Previous results from a transient three-dimensional model of primary pyrolysis were used for the source terms of primary products in this model. A parametric study of reaction atmospheres (H2O, N2, H2, CO2, CO) has been performed. For the N2 and H2O atmosphere, results of the model compared favorably to experimentally obtained yields after the temperature was adjusted to a value higher than that used in experiments. One notable deviation versus experiments is pyrolytic water yield and yield of higher hydrocarbons. The model suggests a not overly strong impact of the reaction atmosphere. However, both chemical and physical effects were observed. Most notably, effects could be seen on the yield of various compounds, temperature profile throughout the reactor system, residence time, radical concentration, and turbulent intensity. At the investigated temperature (873 K), turbulent intensity appeared to have the strongest influence on liquid yield. With the aid of acceleration techniques, most importantly dimension reduction, chemistry agglomeration, and in-situ tabulation, a converged solution could be obtained within a reasonable time (∼30 h). As such, a new potentially useful method has been suggested for numerical analysis of fast pyrolysis.

Comprehensive secondary pyrolysis in fluidized-bed fast pyrolysis of biomass, a fluid dynamics based modelling effort

Homogenous secondary pyrolysis is category of reactions following the primary pyrolysis and presumed important for fast pyrolysis. For the comprehensive chemistry and fluid dynamics, a probability density functional (PDF) approach is used; with a kinetic scheme comprising 134 species and 4169 reactions being implemented. With aid of acceleration techniques, most importantly Dimension Reduction, Chemistry Agglomeration and In-situ Tabulation (ISAT), a solution within reasonable time was obtained. More work is required; however, a solution for levoglucosan (C6H10O5) being fed through the inlet with fluidizing gas at 500 °C, has been obtained. 88.6% of the levoglucosan remained non-decomposed, and 19 different decomposition product species were found above 0.01% by weight. A homogenous secondary pyrolysis scheme proposed can thus be implemented in a CFD environment and acceleration techniques can speed-up the calculation for application in engineering settings.

Geometrical optimization of a fast pyrolysis bubbling fluidized bed reactor using comutational fluid dynamics

This paper analyzes the physical phenomena that take place inside an 1 kg/h bubbling fluidized bed reactor located at Aston University and presents a geometrically modified version of it, in order to improve certain hydrodynamic and gas flow characteristics. The bed uses, in its current operation, 40 L/min of N2 at 520 °C fed through a distributor plate and 15 L/min purge gas stream, i.e., N2 at 20 °C, via the feeding tube. The Eulerian model of FLUENT 6.3 is used for the simulation of the bed hydrodynamics, while the k - ε model accounts for the effect of the turbulence field of one phase on the other. The three-dimensional simulation of the current operation of the reactor showed that a stationary bubble was formed next to the feeding tube. The size of the permanent bubble reaches up to the splash zone of the reactor, without any fluidizaton taking place underneath the feeder. The gas flow dynamics in the freeboard of the reactor is also analyzed. A modified version of the reactor is presented, simulated, and analyzed, together with a discussion on the impact of the flow dynamics on the fast pyrolysis of biomass. © 2010 American Chemical Society.