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.

Mechanisms of convective drying in thick beds of solids

The literature on heat and mass transfer mechanisms in the convective drying of thick beds of solids has been critically reviewed. Related mathematical models of heat transfer are also considered. Experimental and theoretical studies were made of the temperature distribution within beds, and of drying rates, with various materials undergoing convective drying. The experimental work covered thick beds of hygroscopic and non-hygroscopic materials (glass beads of different diameters, polystyrene pellets, activated alumina and wood powder) at air temperatures of 54°C to 84°C. Tests were carried out in a laboratory drying apparatus comprising a wind tunnel through which the air, of controlled temperature and humidity, was passed over a sample suspended from a balance. Thermocouples were inserted at different depths within the sample bed. The temperature distribution profiles for both hygroscopic and non-hygroscopic beds exhibited a clear difference between the temperatures at the surface and bottom during the constant rate period. An effective method was introduced for predicting the critical moisture content. During the falling rate the profiles showed the existence of a receding evaporation plane; this divided the system into a hotter dry zone in the upper section and a wet zone near the bottom. A graphical procedure was established to predict accurately the position of the receding evaporation front at any time. A new mathematical model, based on the receding evaporation front phenomenon, was proposed to predict temperature distributions throughout a bed during drying. Good agreement was obtained when the model was validated by comparing its predictions with experimental data. The model was also able to predict the duration of each drying stage. In experiments using sample trays of different diameters, the drying rate was found to increase with a decrease in the effective length of the bed surface. During the constant rate period with trays of a small effective length, i.e. less than 0.08 m, an 'inversion' in temperature distribution occurred in the bed; the bottom temperature increased and became greater than that of the surface. Experimental measurements were verified in several ways to ensure this phenomenon was real. Theoretical explanations are given for both the effective length and temperature inversion phenomena.

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.

Coalescence of secondary dispersions in packed beds

A study has been made of the coalescence of secondary dispersions in beds of monosized glass ballotini. The variables investigated were superficial velocity, bed depth, ballotini size and dispersed phase concentration. Equipment was designed to generate a toluene ln water dispersion with phase ratios from 0.1 - 1.0 v/v % and whose mean drop size was determined using a Coulter Counter. The coalesced drops were sized by photography and the mean diameter of the effluent drops was determined using a Malvern Particle Size Analyser. Previous models describing single phase flow in porous media are reviewed and it was found that the experimental data obtained in this study is best represented by the Carman-Kozeny equations. Relative permeability correlations were used to predict the saturation profiles across the bed from measured two phase pressure drop data. Theoretical comparison of drop capture mechanisms indicated that direct and indirect interception are predominant. The total capture efficiency for the bed can also be evaluated using Spielman and Fitzpatrick's correlation.The resulting equation is used to predict the initial, local drop capture rate in a coalescer. A mathematical description of the saturation profiles is formulated and verified by the saturation profiles obtained by relative permeability. Based on the Carman-Kozeny equation, an expression is derived analytically to .predict the two phase pressure drop using the parameters which characterise the saturation profiles. By specifying the local saturation at the inlet face for a given velocity and phase ratio, good agreement between experimental pressure drop data and the model predictions was obtained. An attempt to predict the exit drop size has been made using an analogy for flow through non cylindrical channels.

Mechanisms of secondary dispersions separation in particulate beds

The mechanisms by which drops of secondary liquid dispersion ie. <100μ m, are collected, coalesced and transferred have been studied in particulate beds of different sizes and heights of glass ballotini. The apparatus facilitated different coalescer cell arrangements. The liquid-liquid system was toluene/de-ionised water. The inlet drop size distribution was measured by microscopy and using the Malvern Particle Size analyser; the outlet dispersion was sized by photography. The effect of packed height and packing size upon critical velocity, pressure drop and coalescence efficiency have been investigated. Single and two phase flow pressure drops across the packing were correlated by modified Blake-Kozeny equations. Two phase pressure drop was correlated by two equations, one for large ballotini sizes (267μm - 367μm), the other for small ballotini sizes (93μm- 147.5μm). The packings were efficient coalescers up to critical velocities of 3 x 10-2 m/s to 5 x 10-2 m/s. The saturation was measured across the bed using relative permeability and a mathematical model developed which related this profile to measured pressure drops. Filter coefficients for the range of packing studied were found to be accurately predicted from a modified queueing drop model. 

The separation of secondary oil-water dispersions in particulate beds

The literature relating to haze formation, methods of separation, coalescence mechanisms, and models by which droplets <100 μm are collected, coalesced and transferred, have been reviewed with particular reference to particulate bed coalescers. The separation of secondary oil-water dispersions was studied experimentally using packed beds of monosized glass ballotini particles. The variables investigated were superficial velocity, bed depth, particle size, and the phase ratio and drop size distribution of inlet secondary dispersion. A modified pump loop was used to generate secondary dispersions of toluene or Clairsol 350 in water with phase ratios between 0.5-6.0 v/v%.Inlet drop size distributions were determined using a Malvern Particle Size Analyser;effluent, coalesced droplets were sized by photography. Single phase flow pressure drop data were correlated by means of a Carman-Kozeny type equation. Correlations were obtained relating single and two phase pressure drops, as (ΔP2/μc)/ΔP1/μd) = kp Ua Lb dcc dpd Cine A flow equation was derived to correlate the two phase pressure drop data as, ΔP2/(ρcU2) = 8.64*107 [dc/D]-0.27 [L/D]0.71 [dp/D]-0.17 [NRe]1.5 [e1]-0.14 [Cin]0.26  In a comparison between functions to characterise the inlet drop size distributions a modification of the Weibull function provided the best fit of experimental data. The general mean drop diameter was correlated by: q_p q_p p_q /β      Γ ((q-3/β) +1) d qp = d fr  .α        Γ ((P-3/β +1 The measured and predicted mean inlet drop diameters agreed within ±15%. Secondary dispersion separation depends largely upon drop capture within a bed. A theoretical analysis of drop capture mechanisms in this work indicated that indirect interception and London-van der Waal's mechanisms predominate. Mathematical models of dispersed phase concentration m the bed were developed by considering drop motion to be analogous to molecular diffusion.The number of possible channels in a bed was predicted from a model in which the pores comprised randomly-interconnected passage-ways between adjacent packing elements and axial flow occured in cylinders on an equilateral triangular pitch. An expression was derived for length of service channels in a queuing system leading to the prediction of filter coefficients. The insight provided into the mechanisms of drop collection and travel, and the correlations of operating parameters, should assist design of industrial particulate bed coalescers.

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.

The microbial ecology of cow slurry in vermiculture beds

The suitability of cow slurry as a substrate for vermicomposting by Eisenia fetida was investigated. Particular attention was given to the effects of the earthworm on the decomposition and stabilisation of the slurry; and to the interactions between E. fetida and the microflora of the substrate. Assessment of the chemical and microbiological changes in cow slurry stored under forced aeration, and subsequently in shallow trays, showed that neither method was suitable for the treatment of slurry. A comparison of two methods of vermicomposting showed that top-feeding of slurry was more efficient in promoting earthworm growth and cocoon production than the mixing of slurry with solid materials. Management practices were found to have an important influence on the efficiency of the process. An investigation o:f the effect of E. fetida. on the decomposition of slurry indicated that the presence of this earthworm enhanced the stabilisation of the substrate and increased the plant-available nitrogen content. Specific nutritional interactions were observed between E. fetida and micro-organisms in sand/cellulose microcosms. The earthworms were found to be feeding directly upon the cells of certain micro-organisms. Other species were found to be toxic to E. fetida.. A technique was developed :for the production of axenic E. fetida., and the use of such earthworms in :feeding experiments confirmed the importance of some micro-organisms in earthworm nutrition. The seeding of vermiculture beds with one such micro-organism stimulated earthworm growth and consumption of the substrate. Vermicomposted mixtures of cow slurry and spent mushroom compost were shown to have potential application as casing materials in mushroom cultivation. The findings of this study indicate the suitability of vermicomposting as a method for the stabilisation of intensively-produced cow slurry, and give some indication of the importance of micro-organisms in the nutrition of E. fetida.

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.

The flow behaviour of shallow fluidised beds

The flow behaviour of shallow gas-fluidised beds was studied. experimentally using a rotational viscometer, and an inclined open channel. Initially, tests were carried out with the viscometer in order to establish qualitative trends in the flow properties of a variety of materials over a wide range of fluidising conditions. Also, a technique was developed which enabled quantitative viscosity data to be extracted from the experimental results. The flow properties were found to be sensitive to the size, size-range and density of the fluidised material, the type of distributor used, and the moisture content of the fluidising gas. Tests in beds up to 120 mm deep showed that the fluidity of the bed improves with reduction in depth; and indicated a range of flow behaviour from shear-thinning to Newtonian, depending chiefly on fluidising velocity .. Later, an apparatus was built which provided for a steady, continuous flow of fluidised material down an inclined open channel of 3m length x 0.15m square, up to a mass flowrate of 10 kg/s (35 ton/hr). This facility has enabled data to be obtained that is of practical value in industrial applications; which is otherwise difficult in view of the present limited understanding of the true mechanism of fluidised flow. A correlation has been devised, based on analogy with laminar liquid flow, which describes the channel flow behaviour with reasonable accuracy over the whole range of shear-rates used. 1he channeI results indicated that at low fluidiising velocities the flow was adversely affected by settlement of a stagnant layer of particles on to the distributor, which gave rise to increased flow resistance. Conversely, at higher fluidising velocities the resistance at the distributor appeared to be less than at the walls. In view of this, and also because of the disparity in shear-rates between the two types of apparatus, it is not possible as yet to predict exactly the flow behaviour in an open channel from small-scale viscometer tests.