Hydrogen production by steam reforming of DME in a large scale CFB reactor. Part I:computational model and predictions

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether 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). The DME-SR reactions scheme and kinetics in the presence of a bifunctional catalyst of CuO/ZnO/Al2O3+ZSM-5 were incorporated in the model using in-house developed user-defined function. The model was validated by comparing the predictions with experimental data from the literature. The results revealed for the first time detailed CFB reactor hydrodynamics, gas residence time, temperature distribution and product gas composition at a selected operating condition of 300 °C and steam to DME mass ratio of 3 (molar ratio of 7.62). The spatial variation in the gas species concentrations suggests the existence of three distinct reaction zones but limited temperature variations. The DME conversion and hydrogen yield were found to be 87% and 59% respectively, resulting in a product gas consisting of 72 mol% hydrogen. In part II of this study, the model presented here will be used to optimize the reactor design and study the effect of operating conditions on the reactor performance and products.

A computational model of hydrogen production by steam reforming of dimethyl ether in a large scale CFB reactor. Part II:parametric analysis

This study presents a computational parametric analysis of DME steam reforming in a large scale Circulating Fluidized Bed (CFB) reactor. The Computational Fluid Dynamic (CFD) model used, which is based on Eulerian-Eulerian dispersed flow, has been developed and validated in Part I of this study [1]. The effect of the reactor inlet configuration, gas residence time, inlet temperature and steam to DME ratio on the overall reactor performance and products have all been investigated. The results have shown that the use of double sided solid feeding system remarkable improvement in the flow uniformity, but with limited effect on the reactions and products. The temperature has been found to play a dominant role in increasing the DME conversion and the hydrogen yield. According to the parametric analysis, it is recommended to run the CFB reactor at around 300 °C inlet temperature, 5.5 steam to DME molar ratio, 4 s gas residence time and 37,104 ml gcat -1 h-1 space velocity. At these conditions, the DME conversion and hydrogen molar concentration in the product gas were both found to be around 80%.

A mathematical model of coal-fired fluidised bed boilers

A mathematical model of a large coal-fired fluidized bed boiler for power generation is synthesised. The effect of variations in the main parameters of the model on variables such as the background carbon concentrations in the bed, and the transient response of heat evolution are studied. The mechanisms of solids mixing within the bed, combustion and the flow of heat to the boiler tubes are shown to result in a characteristic dynamic response, knowledge of which is essential for the proper control and regulation of a practical system.

Flowback of solids through distribution plates of gas fluidized beds and associated phenomena

This work is concerned with a study of certain phenomena related to the performance and design of distributors in gas fluidized beds with particular regard to flowback of solid particles. The work to be described is divided into two parts. I. In Part one, a review of published material pertaining to distribution plates, including details from the patent specifications, has been prepared. After a chapter on the determination of the incipient fluidizing velocity, the following aspects of multi-orifice distributor plates in gas fluidized beds have been studied: (i) The effect of the distributor on bubble formation related to the way in which even distribution of bubbles on the top surface of the fluidized bed is obtained, e.g. the desirable pressure drop ratio ?PD/?PB for the even distribution of gas across the bed. Ratios of distributor pressure drop ?PD to bed pressure drop at which stable fluidization occurs show reasonable agreement with industrial practice. There is evidence that larger diameter beds tend to be less stable than smaller diameter beds when these are operated with shallow beds. Experiments show that in the presence of the bed the distributor pressure drop is reduced relative to the pressure drop without the bed, and this pressure drop in the former condition is regarded as the appropriate parameter for the design of the distributor. (ii) Experimental measurements of bubble distribution at the surface has been used to indicate maldistribution within the bed. Maldistribution is more likely at low gas flow rates and with distributors having large fractional free area characteristics (i.e. with distributors having low pressure drops). Bubble sizes obtained from this study, as well as those of others, have been successfully correlated. The correlation produced implies the existence of a bubble at the surface of an orifice and its growth by the addition of excess gas from the fluidized bed. (iii) For a given solid system, the amount of defluidized particles stagnating on the distributor plate is influenced by the orifice spacing, bed diameter and gas flow rate, but independent of the initial bed height and the way the orifices are arranged on the distributor plate. II. In Part two, solids flowback through single and multi-orifice distributors in two-dimensional and cylindrical beds of solids fluidized with air has been investigated. Distributors equipped with long cylindrical nozzles have also been included in the study. An equation for the prediction of free flowback of solids through multi-orifice distributors has been derived. Under fluidized conditions two regimes of flowback have been differentiated, namely Jumping and weeping. Data in the weeping regime have been successfully correlated. The limiting gas velocity through the distributor orifices at which flowback is completely excluded is found to be indepnndent of bed height, but a function of distributor design and physical properties of gas and solid used. A criterion for the prediction of this velocity has been established. The decisive advantage of increasing the distributor thickness or using nozzles to minimize solids flowback in fluidized beds has been observed and the opportunity taken to explore this poorly studied subject area. It has been noted, probably for the first time, that with long nozzles, there exists a critical nozzle length above which uncontrollable downflow of solids occurs. A theoretical model for predicting the critical length of a bundle of nozzles in terms of gas velocity through the nozzles has been set up. Theoretical calculations compared favourably with experiments.

Total solution of Gibilaro and Rowe model for a segregating fluidised bed

This study re-examines the one-dimensional equilibrium model of Gibilaro and Rowe (1974) for a segregating gas fluidized bed. The model was based on volumetric jetsam concentration and divided the bed contents into bulk and wake phases, taking account of bulk and wake flux, segregation, exchange between the bulk and wake phases, and axial mixing. Due to the complex nature of the model and its unstable solution, the lack of computing power at the time prevented the authors from doing little more than the analytical solutions to specific cases of this model. This paper provides a numerical total solution and allows the effect of the respective parameters to be compared for the first time. There is also a comparison with experimental results, which showed a reasonable agreement.

Modelling of transient heat transfer in annealing furnaces

This is a study of heat transfer in a lift-off furnace which is employed in the batch annealing of a stack of coils of steel strip. The objective of the project is to investigate the various factors which govern the furnace design and the heat transfer resistances, so as to reduce the time of the annealing cycle, and hence minimize the operating costs. The work involved mathematical modelling of patterns of gas flow and modes of heat transfer. These models are: Heat conduction and its conjectures in the steel coils;Convective heat transfer in the plates separating the coils in the stack and in other parts of the furnace; and Radiative and convective heat transfer in the furnace by using the long furnace model. An important part of the project is the development of numerical methods and computations to solve the transient models. A limited number of temperature measurements was available from experiments on a test coil in an industrial furnace. The mathematical model agreed well with these data. The model has been used to show the following characteristics of annealing furnaces, and to suggest further developments which would lead to significant savings: - The location of the limiting temperature in a coil is nearer to the hollow core than to the outer periphery. - Thermal expansion of the steel tends to open the coils, reduces their thermal conductivity in the radial direction, and hence prolongs the annealing cycle. Increasing the tension in the coils and/or heating from the core would overcome this heat transfer resistance. - The shape and dimensions of the convective channels in the plates have significant effect on heat convection in the stack. An optimal design of a channel is shown to be of a width-to-height ratio equal to 9. - Increasing the cooling rate, by using a fluidized bed instead of the normal shell and tube exchanger, would shorten the cooling time by about 15%, but increase the temperature differential in the stack. - For a specific charge weight, a stack of different-sized coils will have a shorter annealing cycle than one of equally-sized coils, provided that production constraints allow the stacking order to be optimal. - Recycle of hot flue gases to the firing zone of the furnace would produce a. decrease in the thermal efficiency up to 30% but decreases the heating time by about 26%.

Sequential pyrolysis of willow SRC at low and high heating rates:implications for selective pyrolysis

The main aim of the work is to investigate sequential pyrolysis of willow SRC using two different heating rates (25 and 1500 °C/min) between 320 and 520 °C. Thermogravimetric analysis (TGA) and pyrolysis - gas chromatography - mass spectroscopy (Py-GC-MS) have been used for this analysis. In addition, laboratory scale processing has been undertaken to compare product distribution from fast and slow pyrolysis at 500 °C. Fast pyrolysis was carried out using a 1 kg/h continuous bubbling fluidized bed reactor, and slow pyrolysis using a 100 g batch reactor. Findings from this study show that heating rate and pyrolysis temperatures have a significant influence on the chemical content of decomposition products. From the analytical sequential pyrolysis, an inverse relationship was seen between the total yield of furfural (at high heating rates) and 2-furanmethanol (at low heating rates). The total yield of 1,2-dihydroxybenzene (catechol) was found to be significant higher at low heating rates. The intermediates of catechol, 2-methoxy-4-(2-propenyl)phenol (eugenol); 2-methoxyphenol (guaiacol); 4-Hydroxy-3,5-dimethoxybenzaldehyde (syringaldehyde) and 4-hydroxy-3-methoxybenzaldehyde (vanillin), were found to be highest at high heating rates. It was also found that laboratory scale processing alters the pyrolysis bio-oil chemical composition, and the proportions of pyrolysis product yields. The GC-MS/FID analysis of fast and slow pyrolysis bio-oils reveals significant differences. © 2011 Elsevier Ltd. All rights reserved.

Lignin pyrolysis: results from an international collaboration

An international study of fast pyrolysis of lignin was undertaken. Fourteen laboratories in eight different countries contributed. Two lignin samples were distributed to the laboratories for analysis and bench-scale process testing in fast pyrolysis. Analyses included proximate and ultimate analysis, thermogravimetric analysis, and analytical pyrolysis. The bench-scale test included bubbling fluidized-bed reactors and entrained-flow systems. Based on the results of the various analyses and tests it was concluded that a concentrated lignin (estimated at about 50% lignin and 50% cellulose) behaved like a typical biomass, producing a slightly reduced amount of a fairly typical bio-oil, while a purified lignin material was difficult to process in the fast pyrolysis reactors and produced a much lower amount of a different kind of bio-oil. It was concluded that for highly concentrated lignin feedstocks new reactor designs will be required other than the typical fluidized-bed fast pyrolysis systems.

Study on the pyrolytic behaviour of xylan-based hemicellulose using TG-FTIR and Py-GC-FTIR

Two sets of experiments, categorized as TG–FTIR and Py–GC–FTIR, are employed to investigate the mechanism of the hemicellulose pyrolysis and the formation of main gaseous and bio-oil products. The “sharp mass loss stage” and the corresponding evolution of the volatile products are examined by the TG–FTIR graphs at the heating rate of 3–80 K/min. A pyrolysis unit, composed of fluidized bed reactor, carbon filter, vapour condensing system and gas storage, is employed to investigate the products of the hemicellulose pyrolysis under different temperatures (400–690 °C) at the feeding flow rate of 600 l/h. The effects of temperature on the condensable products are examined thoroughly. The possible routes for the formation of the products are systematically proposed from the primary decomposition of the three types of unit (xylan, O-acetylxylan and 4-O-methylglucuronic acid) and the secondary reactions of the fragments. It is found that the formation of CO is enhanced with elevated temperature, while slight change is observed for the yield of CO2 which is the predominant products in the gaseous mixture.

Numerical Comparison of the Drag Models of Granular Flows Applied to the Fast Pyrolysis of Biomass

The paper presents a comparison between the different drag models for granular flows developed in the literature and the effect of each one of them on the fast pyrolysis of wood. The process takes place on an 100 g/h lab scale bubbling fluidized bed reactor located at Aston University. FLUENT 6.3 is used as the modeling framework of the fluidized bed hydrodynamics, while the fast pyrolysis of the discrete wood particles is incorporated as an external user defined function (UDF) hooked to FLUENT’s main code structure. Three different drag models for granular flows are compared, namely the Gidaspow, Syamlal O’Brien, and Wen-Yu, already incorporated in FLUENT’s main code, and their impact on particle trajectory, heat transfer, degradation rate, product yields, and char residence time is quantified. The Eulerian approach is used to model the bubbling behavior of the sand, which is treated as a continuum. Biomass reaction kinetics is modeled according to the literature using a two-stage, semiglobal model that takes into account secondary reactions.

Time scale analysis for fluidizedbedmeltgranulation-II: binder spreading rate

The spreading time of liquid binder droplet on the surface a primary particle is analyzed for Fluidized Bed Melt Granulation (FBMG). As discussed in the first paper of this series (Chua et al., in press) the droplet spreading rate has been identified as one of the important parameters affecting the probability of particles aggregation in FBMG. In this paper, the binder droplet spreading time has been estimated using Computational Fluid Dynamic modeling (CFD) based on Volume of Fluid approach (VOF). A simplified analytical solution has been developed and tested to explore its validity for predicting the spreading time. For the purpose of models validation, the droplet spreading evolution was recorded using a high speed video camera. Based on the validated model, a generalized correlative equation for binder spreading time is proposed. For the operating conditions considered here, the spreading time for Polyethylene Glycol (PEG1500) binder was found to fall within the range of 10-2 to 10-5 s. The study also included a number of other common binders used in FBMG. The results obtained here will be further used in paper III, where the binder solidification rate is discussed.

Time scale analysis for fluidizedbedmeltgranulation III: binder solidification rate

In series I and II of this study ([Chua et al., 2010a] and [Chua et al., 2010b]), we discussed the time scale of granule–granule collision, droplet–granule collision and droplet spreading in Fluidized Bed Melt Granulation (FBMG). In this third one, we consider the rate at which binder solidifies. Simple analytical solution, based on classical formulation for conduction across a semi-infinite slab, was used to obtain a generalized equation for binder solidification time. A multi-physics simulation package (Comsol) was used to predict the binder solidification time for various operating conditions usually considered in FBMG. The simulation results were validated with experimental temperature data obtained with a high speed infrared camera during solidification of ‘macroscopic’ (mm scale) droplets. For the range of microscopic droplet size and operating conditions considered for a FBMG process, the binder solidification time was found to fall approximately between 10-3 and 10-1 s. This is the slowest compared to the other three major FBMG microscopic events discussed in this series (granule–granule collision, granule–droplet collision and droplet spreading).

Fast pyrolysis and nitrogenolysis of biomass and biogenic residues:production of a sustainable slow release fertiliser

The production of agricultural and horticultural products requires the use of nitrogenous fertiliser that can cause pollution of surface and ground water and has a large carbon footprint as it is mainly produced from fossil fuels. The overall objective of this research project was to investigate fast pyrolysis and in-situ nitrogenolysis of biomass and biogenic residues as an alternative route to produce a sustainable solid slow release fertiliser mitigating the above stated problems. A variety of biomasses and biogenic residues were characterized by proximate analysis, ultimate analysis, thermogravimetric analysis (TGA) and Pyrolysis – Gas chromatography – Mass Spectroscopy (Py–GC–MS) for their potential use as feedstocks using beech wood as a reference material. Beech wood was virtually nitrogen free and therefore suitable as a reference material as added nitrogen can be identified as such while Dried Distillers Grains with Solubles (DDGS) and rape meal had a nitrogen content between 5.5wt.% and 6.1wt.% qualifying them as high nitrogen feedstocks. Fast pyrolysis and in-situ nitrogenolysis experiments were carried out in a continuously fed 1kg/h bubbling fluidized bed reactor at around 500°C quenching the pyrolysis vapours with isoparaffin. In-situ nitrogenolysis experiments were performed by adding ammonia gas to the fast pyrolysis reactor at nominal nitrogen addition rates between 5wt.%C and 20wt.%C based on the dry feedstock’s carbon content basis. Mass balances were established for the processing experiments. The fast pyrolysis and in-situ nitrogenolysis products were characterized by proximate analysis, ultimate analysis and GC– MS. High liquid yields and good mass balance closures of over 92% were obtained. The most suitable nitrogen addition rate for the in-situ nitrogenolysis experiments was determined to be 12wt.%C on dry feedstock carbon content basis. However, only a few nitrogen compounds that were formed during in-situ nitrogenolysis could be identified by GC–MS. A batch reactor process was developed to thermally solidify the fast pyrolysis and in-situ nitrogenolysis liquids of beech wood and Barley DDGS producing a brittle solid product. This was obtained at 150°C with an addition of 2.5wt% char (as catalyst) after a processing time of 1h. The batch reactor was also used for modifying and solidifying fast pyrolysis liquids derived from beech wood by adding urea or ammonium phosphate as post processing nitrogenolysis. The results showed that this type of combined approach was not suitable to produce a slow release fertiliser, because the solid product contained up to 65wt.% of highly water soluble nitrogen compounds that would be released instantly by rain. To complement the processing experiments a comparative study via Py–GC–MS with inert and reactive gas was performed with cellulose, hemicellulose, lignin and beech wood. This revealed that the presence of ammonia gas during analytical pyrolysis did not appear to have any direct impact on the decomposition products of the tested materials. The chromatograms obtained showed almost no differences between inert and ammonia gas experiments indicating that the reaction between ammonia and pyrolysis vapours does not occur instantly. A comparative study via Fourier Transformed Infrared Spectroscopy of solidified fast pyrolysis and in-situ nitrogenolysis products showed that there were some alterations in the spectra obtained. A shift in frequencies indicating C=O stretches typically related to the presence of carboxylic acids to C=O stretches related to amides was observed and no double or triple bonded nitrogen was detected. This indicates that organic acids reacted with ammonia and that no potentially harmful or non-biodegradable triple bonded nitrogen compounds were formed. The impact of solid slow release fertiliser (SRF) derived from pyrolysis and in-situ nitrogenolysis products from beech wood and Barley DDGS on microbial life in soils and plant growth was tested in cooperation with Rothamsted Research. The microbial incubation tests indicated that microbes can thrive on the SRFs produced, although some microbial species seem to have a reduced activity at very high concentrations of beech wood and Barley DDGS derived SRF. The plant tests (pot trials) showed that the application of SRF derived from beech wood and barley DDGS had no negative impact on germination or plant growth of rye grass. The fertilizing effect was proven by the dry matter yields in three harvests after 47 days, 89 days and 131 days. The findings of this research indicate that in general a slow release fertiliser can be produced from biomass and biogenic residues by in-situ nitrogenolysis. Nevertheless the findings also show that additional research is necessary to identify which compounds are formed during this process.

A comparative study of straw, perennial grasses and hardwoods in terms of fast pyrolysis products

The aim of this study is to characterise and compare fast pyrolysis product yields from straw, high yielding perennial grasses and hardwoods. Feedstocks selected for this study include: wheat straw (Triticum aestivum), switch grass (Panicum virgatum), miscanthus (Miscanthus x giganteus), willow short rotation coppice (Salix viminalis) and beech wood (Fagus sylvatica). The experimental work is divided into two sections: analytical (TGA and Py-GC-MS) and laboratory scale processing using a continuously fed bubbling fluidized bed reactor with a capacity of up to 1 kg/h. Pyrolysis-gas chromatography-mass spectrometry (Py-GC-MS) has been used to quantify pyrolysis products and simulate fast pyrolysis heating rates, in order to study potential key light and medium volatile decomposition products found in these feedstocks. Py-GC-MS quantification results show that the highest yields of furfural (0.57 wt.%), 2-furanmethanol (0.18 wt.%), levoglucosan (0.73 wt.%), 1,2-benzenediol (0.27 wt.%) and 2-methoxy-4-vinylphenol (0.38 wt.%) were found in switch grass, and that willow SRC produced the highest yield of phenol (0.33 wt.%). The bio-oil higher heating value was highest for switch grass (22.3 MJ/kg). Water content within the bio-oil is highest in the straw and perennial grasses and lowest in the hardwood willow SRC. The high bio-oil and char heating value and low water content found in willow SRC, makes this crop an attractive energy feedstock for fast pyrolysis processing, if the associated production costs and harvest yields can be maintained at current reported values. The bio-oil from switch grass has the highest potential for the production of high value chemicals. © 2013 Elsevier Ltd. All rights reserved.

Novel biomass conversion routes:ammonia from biomass, and marine macroalgae for energy

This thesis presents a techno-economic investigation of the generation of electricity from marine macroalgae (seaweed) in the UK (Part 1), and the production of anhydrous ammonia from synthesis gas (syngas) generated from biomass gasification (Part 2). In Part 1, the study covers the costs from macroalgae production to the generation of electricity via a CHP system. Seven scenarios, which varied the scale and production technique, were investigated to determine the most suitable scale of operation for the UK. Anaerobic digestion was established as the most suitable technology for macroalgae conversion to CHP, based on a number of criteria. All performance and cost data have been taken from published literature. None of the scenarios assessed would be economically viable under present conditions, although the use of large-scale electricity generation has more potential than small-scale localised production. Part 2 covers the costs from the delivery of the wood chip feedstock to the production of ammonia. Four cases, which varied the gasification process used and the scale of production, were investigated to determine the most suitable scale of operation for the UK. Two gasification processes were considered, these were O2-enriched air entrained flow gasification and Fast Internal Circulating Fluidised Bed. All performance and cost data have been taken from published literature, unless otherwise stated. Large-scale (1,200 tpd) ammonia production using O2-enriched air entrained flow gasification was determined as the most suitable system, producing the lowest ammonia-selling price, which was competitive to fossil fuels. Large-scale (1,200 tpd) combined natural gas/biomass syngas ammonia production also generated ammonia at a price competitive to fossil fuels.