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.

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.

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.

Pyrolysis of Napier grass in a fixed bed reactor:effect of operating conditions on product yields and characteristic

This study presents a report on pyrolysis of Napier grass stem in a fixed bed reactor. The effects of nitrogen flow (20 to 60 mL/min), and reaction temperature (450 to 650 degrees C) were investigated. Increasing the nitrogen flow from 20 to 30 mL/min increased the bio-oil yield and decreased both bio-char and non-condensable gas. 30 mL/min nitrogen flow resulted in optimum bio-oil yield and was used in the subsequent experiments. Reaction temperatures between 450 and 600 degrees C increased the bio-oil yield, with maximum yield of 32.26 wt% at 600 degrees C and a decrease in the corresponding bio-char and non-condensable gas. At 650 degrees C, reductions in the bio-oil and bio-char yields were recorded while the non-condensable gas increased. Water content of the bio-oil decreased with increasing reaction temperature, while density and viscosity increased. The observed pH and higher heating values were between 2.43 to 2.97, and 25.25 to 28.88 MJ/kg, respectively. GC-MS analysis revealed that the oil was made up of highly oxygenated compounds and requires upgrading. The bio-char and non-condensable gas were characterized, and the effect of reaction temperature on the properties was evaluated. Napier grass represents a good source of renewable energy when all pyrolysis products are efficiently utilized.

Coupling of Heck and hydrogenation reactions in a continuous compact reactor

A continuous multi-step synthesis of 1,2-diphenylethane was performed sequentially in a structured compact reactor. This process involved a Heck C-C coupling reaction followed by the addition of hydrogen to perform reduction of the intermediate obtained in the first step. Both of the reactions were catalysed by microspherical carbon-supported Pd catalysts. Due to the integration of the micro-heat exchanger, the static mixer and the mesoscale packed-bed reaction channel, the compact reactor was proven to be an intensified tool for promoting the reactions. In comparison with the batch reactor, this flow process in the compact reactor was more efficient as: (i) the reaction time was significantly reduced (ca. 7 min versus several hours), (ii) no additional ligands were used and (iii) the reaction was run at lower operational pressure and temperature. Pd leached in the Heck reaction step was shown to be effectively recovered in the following hydrogenation reaction section and the catalytic activity of the system can be mostly retained by reverse flow operation. © 2009 Elsevier Inc. All rights reserved.

Gas distribution in shallow large diameter packed beds

It is important to maintain a uniform distribution of gas and liquid in large diameter packed columns to maintain mass transfer efficiency on scaling up. This work presents measurements and methods of evaluating maldistributed gas flow in packed columns. Little or no previous work has been done in this field. A gas maldistribution number, F, was defined, based on point to point velocity variations in the gas emerging from the top of packed beds. f has a minimum value for a uniformly distributed flow and much larger values for maldistributed flows. A method of testing the quality of vapour distributors is proposed, based on "the variation of f with packed height. A good gas distributor requires a short packed depth to give a good gas distribution. Measurements of gas maldistribution have shown that the principle of dynamic similarity is satisfied if two geometrically similar beds are operated at the same Reynold's number. The validity of f as a good measure of gas maldistribution, and the principle of dynamic similarity are tested statistically by Multi-Factor Analysis of the variance, and visually by the response "surfaces technique. Pressure distribution has been measured in a model of a large diameter packed bed, and shown to be associated with the velocity of the gas in a tangential feed pipe. Two simplified theoretical models are proposed to describe the flow of gases through packed beds and to support the principle of dynamic similarity. These models explain why the packed bed itself causes the flow of gas to become more uniformly distributed. A 1.2m. diameter scaled-down model was constructed geometrically similar to a 7.3m. diameter vacuum crude distillation column. The previously known internal cylinder gas distributor was tested. Three new distributors suitable for use in a large diameter column were developed and tested, these are: Internal Cylinder with Slots and Cross Baffles, Internal Cylinder with Guides in the Annulus, Internal Cylinder with Internal Cross Baffles - It has been shown that this is an excellent distributor.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors. Part B: Heat, momentum and mass transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the heat, momentum and mass transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user-defined function (UDF). The study completes the fast pyrolysis modelling in bubbling fluidised bed reactors.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors, Part A: Eulerian computation of momentum transport in bubbling fluidised beds

The fluid–particle interaction inside a 150 g/h fluidised bed reactor is modelled. The biomass particle is injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase, according to the literature. FLUENT 6.2 has been used as the modelling framework of the simulations with a completely revised drag model, in the form of user defined function (UDF), to calculate the forces exerted on the particle as well as its velocity components. 2-D and 3-D simulations are tested and compared. The study is the first part of a complete pyrolysis model in fluidised bed reactors.

CFD modelling of the fast pyrolysis of biomass in fluidised bed reactors: modelling the impact of biomass shrinkage

The fluid–particle interaction and the impact of shrinkage on pyrolysis of biomass inside a 150 g/h fluidised bed reactor is modelled. Two 500 View the MathML sourcem in diameter biomass particles are injected into the fluidised bed with different shrinkage conditions. The two different conditions consist of (1) shrinkage equal to the volume left by the solid devolatilization, and (2) shrinkage parameters equal to approximately half of particle volume. The effect of shrinkage is analysed in terms of heat and momentum transfer as well as product yields, pyrolysis time and particle size considering spherical geometries. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Heat transfer from the bubbling bed to the discrete biomass particle, as well as biomass reaction kinetics are modelled according to the literature. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of user defined function (UDF).

A CFD approach on the effect of particle size on char entrainment in bubbling fluidised bed reactors

The fluid – particle interaction inside a 41.7 mg s-1 fluidised bed reactor is modelled. Three char particles of sizes 500 µm, 250 µm, and 100 µm are injected into the fluidised bed and the momentum transport from the fluidising gas and fluidised sand is modelled. Due to the fluidising conditions and reactor design the char particles will either be entrained from the reactor or remain inside the bubbling bed. The particle size is the factor that differentiates the particle motion inside the reactor and their efficient entrainment out of it. A 3-Dimensional simulation has been performed with a completele revised momentum transport model for bubble three-phase flow according to the literature as an extension to the commercial finite volume code FLUENT 6.2.

Computational modelling of the impact of particle size to the heat transfer coefficient between biomass particles and a fluidised bed

The fluid–particle interaction and the impact of different heat transfer conditions on pyrolysis of biomass inside a 150 g/h fluidised bed reactor are modelled. Two different size biomass particles (350 µm and 550 µm in diameter) are injected into the fluidised bed. The different biomass particle sizes result in different heat transfer conditions. This is due to the fact that the 350 µm diameter particle is smaller than the sand particles of the reactor (440 µm), while the 550 µm one is larger. The bed-to-particle heat transfer for both cases is calculated according to the literature. Conductive heat transfer is assumed for the larger biomass particle (550 µm) inside the bed, while biomass–sand contacts for the smaller biomass particle (350 µm) were considered unimportant. The Eulerian approach is used to model the bubbling behaviour of the sand, which is treated as a continuum. Biomass reaction kinetics is modelled according to the literature using a two-stage, semi-global model which takes into account secondary reactions. The particle motion inside the reactor is computed using drag laws, dependent on the local volume fraction of each phase. FLUENT 6.2 has been used as the modelling framework of the simulations with the whole pyrolysis model incorporated in the form of User Defined Function (UDF).

Gas distribution in shallow packed beds

Packed beds have many industrial applications and are increasingly used in the process industries due to their low pressure drop. With the introduction of more efficient packings, novel packing materials (i.e. adsorbents) and new applications (i.e. flue gas desulphurisation); the aspect ratio (height to diameter) of such beds is decreasing. Obtaining uniform gas distribution in such beds is of crucial importance in minimising operating costs and optimising plant performance. Since to some extent a packed bed acts as its own distributor the importance of obtaining uniform gas distribution has increased as aspect ratios (bed height to diameter) decrease. There is no rigorous design method for distributors due to a limited understanding of the fluid flow phenomena and in particular of the effect of the bed base / free fluid interface. This study is based on a combined theoretical and modelling approach. The starting point is the Ergun Equation which is used to determine the pressure drop over a bed where the flow is uni-directional. This equation has been applied in a vectorial form so it can be applied to maldistributed and multi-directional flows and has been realised in the Computational Fluid Dynamics code PHOENICS. The use of this equation and its application has been verified by modelling experimental measurements of maldistributed gas flows, where there is no free fluid / bed base interface. A novel, two-dimensional experiment has been designed to investigate the fluid mechanics of maldistributed gas flows in shallow packed beds. The flow through the outlet of the duct below the bed can be controlled, permitting a rigorous investigation. The results from this apparatus provide useful insights into the fluid mechanics of flow in and around a shallow packed bed and show the critical effect of the bed base. The PHOENICS/vectorial Ergun Equation model has been adapted to model this situation. The model has been improved by the inclusion of spatial voidage variations in the bed and the prescription of a novel bed base boundary condition. This boundary condition is based on the logarithmic law for velocities near walls without restricting the velocity at the bed base to zero and is applied within a turbulence model. The flow in a curved bed section, which is three-dimensional in nature, is examined experimentally. The effect of the walls and the changes in gas direction on the gas flow are shown to be particularly significant. As before, the relative amounts of gas flowing through the bed and duct outlet can be controlled. The model and improved understanding of the underlying physical phenomena form the basis for the development of new distributors and rigorous design methods for them.

Element characteristic tolerance for semi-batch fixed bed biomass pyrolysis

Biomass pyrolysis to bio-oil is one of the promising sustainable fuels. In this work, relation between biomass feedstock element characteristic and crude bio-oil production yield and lower heating value was explored. The element characteristics considered in this study include moisture, ash, fix carbon, volatile matter, C, H, N, O, S, cellulose, hemicellulose, and lignin content. A semi-batch fixed bed reactor was used for biomass pyrolysis with heating rate of 30 °C/min from room temperature to 600 °C and the reactor was held at 600 °C for 1 h before cooling down. Constant nitrogen flow (1bar) was provided for anaerobic condition. Sago and Napier glass were used in the study to create different element characteristic of feedstock by altering mixing ratio. Comparison between each element characteristic to crude bio-oil yield and low heating value was conducted. The result suggested potential key element characteristic for pyrolysis and provide a platform to access the feedstock element acceptance range.

Application of staged O2 feeding in the oxidative dehydrogenation of ethylbenzene to styrene over Al2O3 and P2O5/SiO2 catalysts

Drastic improvements in styrene yield and selectivity were achieved in the oxidative dehydrogenation of ethylbenzene by staged feeding of O2. Six isothermal packed bed reactors were used in series with intermediate feeding of O2, while all EB was fed to the first reactor, diluted with helium or CO2 (1:5 molar ratio), resulting in total O2:EB molar feed ratios of 0.2-0.6. The two catalyst samples, γ-Al 2O3 and 5P/SiO2, that were applied both benefitted from this operation mode. The ethylbenzene conversion per stage and the selectivity to styrene were significantly improved. The production of COX was effectively reduced, while the selectivity to other side products remained unchanged. Compared with co-feeding at a total O 2:EB molar feed ratio of 0.6, by staged feeding the EB conversion (+15% points for both catalysts), ST selectivity (+4% points for both samples) and O2 (ST) selectivity (+9% points for γ-Al2O 3 and +17% points for 5P/SiO2) all improved. The ethylbenzene conversion over 5P/SiO2 can be increased from 18% to 70% by increasing the number of reactors from 1 to 6 with each reactor a total amount of O2 of 0.1 without the loss of ST selectivity (93%). For 5P/SiO2 a higher temperature (500 C vs. 450 C for Al 2O3) is required. Essentially more catalyst (5P/SiO 2) was required to achieve full O2 conversion in each reactor. Staged feeding of O2 does not eliminate the existing issues of the catalyst stability both in time-on stream and as a function of the number of catalyst regenerations (5P/SiO2), or the relatively moderate performance (relatively low styrene selectivity for γ-Al2O 3). © 2014 Elsevier B.V.

Coke formation in the oxidative dehydrogenation of ethylbenzene to styrene by TEOM

A packed bed microbalance reactor setup (TEOM-GC) is used to investigate the formation of coke as a function of time-on-stream on γ-Al2O3 and 3P/SiO2 catalyst samples under different conditions for the ODH reaction of ethylbenzene to styrene. All samples show a linear correlation of the styrene selectivity and yield with the initial coverage of coke. The COX production increases with the coverage of coke. On the 3 wt% P/SiO2 sample, the initial coke build-up is slow and the coke deposition rate increases with time. On alumina-based catalyst samples, a fast initial coke build-up takes place, decreasing with time-on-stream, but the amount of coke does not stabilize. A higher O2 : EB feed ratio results in more coke, and a higher temperature results in less coke. This coking behaviour of Al2O3 can be described by existing "monolayer-multilayer" models. Further, the coverage of coke on the catalyst varies with the position in the bed. For maximal styrene selectivity, the optimal coverage of coke should be sufficient to convert all O2, but as low as possible to prevent selectivity loss by COX production. This is in favour of high temperature and low O2 : EB feed ratios. The optimal coke coverage depends in a complex way on all the parameters: temperature, the O2 : EB feed ratio, reactant concentrations, and the type of starting material. This journal is