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).

Bioreaction and separation in batch chromatographic columns

The objective of this work has been to study the behaviour and performance of a batch chromatographic column under simultaneous bioreaction and separation conditions for several carbohydrate feedstocks. Four bioreactions were chosen, namely the hydrolysis of sucrose to glucose and fructose using the enzyme invertase, the hydrolysis of inulin to fructose and glucose using inulinase, the hydrolysis of lactose to glucose and galactose using lactase and the isomerization of glucose to fructose using glucose isomerase. The chromatographic columns employed were jacketed glass columns ranging from 1 m to 2 m long and the internal diameter ranging from 0.97 cm to 1.97 cm. The stationary phase used was a cation exchange resin (PUROLITE PCR-833) in the Ca2+ form for the hydrolysis and the Mg2+ form for the isomerization reactions. The mobile phase used was a diluted enzyme solution which was continuously pumped through the chromatographic bed. The substrate was injected at the top of the bed as a pulse. The effect of the parameters pulse size, the amount of substrate solution introduced into the system corresponding to a percentage of the total empty column volume (% TECV), pulse concentration, eluent flowrate and the enzyme activity of the eluent were investigated. For the system sucrose-invertase complete conversions of substrate were achieved for pulse sizes and pulse concentrations of up to 20% TECV and 60% w/v, respectively. Products with purity above 90% were obtained. The enzyme consumption was 45% of the amount theoretically required to produce the same amount of product as in a conventional batch reactor. A value of 27 kg sucrose/m3 resin/h for the throughput of the system was achieved. The systematic investigation of the factors affecting the performance of the batch chromatographic bioreactor-separator was carried out by employing a factorial experimental procedure. The main factors affecting the performance of the system were the flowrate and enzyme activity. For the system inulin-inulinase total conversions were also obtained for pulses sizes of up to 20 % TECV and a pulse concentration of 10 % w/v. Fructose rich fractions with 100 % purity and representing up to 99.4 % of the total fructose generated were obtained with an enzyme consumption of 32 % of the amount theoretically required to produce the same amount of product in a conventional batch reactor. The hydrolysis of lactose by lactase was studied in the glass columns and also in an SCCR-S unit adapted for batch operation, in co-operation with Dr. Shieh, a fellow researcher in the Chemical Engineering and Applied Chemistry Department at Aston University. By operating at up to 30 % w/v lactose feed concentrations complete conversions were obtained and the purities of the products generated were above 90%. An enzyme consumption of 48 % of the amount theoretically required to produce the same amount of product in a conventional batch reactor was achieved. On working with the system glucose-glucose isomerase, which is a reversible reaction, the separation obtained with the stationary phase conditioned in the magnesium form was very poor although the conversion obtained was compatible with those for conventional batch reactors. By working with a mixed pulse of enzyme and substrate, up to 82.5 % of the fructose generated with a purity of 100 % was obtained. The mathematical modelling and computer simulation of the batch chromatographic bioreaction-separation has been performed on a personal computer. A finite difference method was used to solve the partial differential equations and the simulation results showed good agreement with the experimental results.

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.

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).

Scale-up of human mesenchymal stem cell culture:current technologies and future challenges

Regenerative medicine technologies have the potential to revolutionise human healthcare. However, whilst science has revealed the potential, and early products have shown the power of such therapies, there is now a need for the long-term supply of human stem cells in sufficient numbers to create reproducible and cost effective therapeutic products. The industrial platforms to be developed for human cell culture are in some ways analogous to those already developed for biopharmaceutical production using mammalian cells at large scales. However, there are a number of unique challenges that need to be addressed, largely because the quality of the cell is paramount, rather than the proteins that they express. © 2013 Elsevier Ltd.

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.

Analysis of powder caking in multicomponent powders using atomic force microscopy to examine particle properties

Atomic force microscopy has been used to study the surface properties of model spray dried powders. Phase imaging, nanoindentation and force modulation microscopy have differentiated between the different surface material properties of the particles, revealing a regular dispersion of soft, oil rich areas distributed across the particles' surface. Humidity and temperature cycling effects on the caking behavior of the particles have also been investigated, with significant morphology changes and onset of caking found to occur within relatively short periods of time.

Caking behaviour of spray-dried milk powders:using scanning probe microscopy to study nanoscale surface properties and material composition

Spray drying is widely used to manufacture many powdered products, with the drying process parameters having significant influence over the final powder's surface properties and propensity for unwanted caking. In most cases caking experiments are performed on bulk powders, but especially in multi-component powders, it is often difficult to interpret these results, where interaction effects between particles can be complex. Here the technique of scanning probe microscopy is used to characterize the nanoscale properties of spray dried model milk powders in order to investigate the surface properties of the powders.

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).

The use of atomic force microscopy in investigating particle caking mechanisms

Spray-dried materials are being used increasingly in industries such as food, detergent and pharmaceutical manufacture. Spray-dried sodium carbonate is an important product that has a great propensity to cake; its moisture-sorption properties are very different to the crystalline and amorphous species, with a great affinity for atmospheric moisture. This work demonstrates how the noncontact surface analysis of individual particles using atomic force microscopy can highlight the possible mechanisms of unwanted agglomeration. The nondestructive nature of this method allows cycling of localised humidity in situ and repeated scanning of the same particle area. The resulting topography and phase scans showed that humidity cycling caused changes in the distribution of material phases that were not solely dependent on topographical changes. © 2013 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.

Gas-liquid simulation of an airlift bubble column reactor

The simulation of two-phase flow for an experimental airlift reactor (32-l volume) using commercially available software from Fluent Incorporated is presented here (http://www.fluent.co.uk). Data from the simulation is compared with the experimental data obtained by the tracking of a magnetic particle and analysis of the pressure drop to determine the gas hold-up. Comparisons between vertical velocity and gas hold-up were made for a series of experiments where the superficial gas velocity in the riser was adjusted between 0.01 and 0.075 m s-1. © 2003 Elsevier B.V. All rights reserved.

Gas-liquid-solid flow modelling in a bubble column

An alternative approach to the modelling of solid-liquid and gas-liquid-solid flows for a 5:1 height to width aspect ratio bubble column is presented here. A modified transport equation for the volume fraction of a dispersed phase has been developed for the investigation of turbulent buoyancy driven flows (Chem. Eng. Proc., in press). In this study, a modified transport equation has been employed for discrete phase motion considering both solid-liquid and gas-liquid-solid flows. The modelling of the three-phase flow in a bubble column was achieved in the following case: injecting a slug of solid particles into the column for 10 s at a velocity of 0.1 m s-1 and then the gas phase flow was initiated with a superficial gas velocity of 0.02 cm s-1. © 2003 Elsevier B.V. All rights reserved.

The modelling of buoyancy driven flow in bubble columns

Using the analogy between lateral convection of heat and the two-phase flow in bubble columns, alternative turbulence modelling methods were analysed. The k-ε turbulence and Reynolds stress models were used to predict the buoyant motion of fluids where a density difference arises due to the introduction of heat or a discrete phase. A large height to width aspect ratio cavity was employed in the transport of heat and it was shown that the Reynolds stress model with the use of velocity profiles including the laminar flow solution resulted in turbulent vortices developing. The turbulence models were then applied to the simulation of gas-liquid flow for a 5:1 height to width aspect ratio bubble column. In the case of a gas superficial velocity of 0.02 m s-1 it was determined that employing the Reynolds stress model yielded the most realistic simulation results. © 2003 Elsevier B.V. All rights reserved.