Comparison of the effect of pre-treatment and catalysts on liquid quality from fast pyrolysis of biomass

The overall objective of this work was to compare the effect of pre-treatment and catalysts on the quality of liquid products from fast pyrolysis of biomass. This study investigated the upgrading of bio-oil in terms of its quality as a bio-fuel and/or source of chemicals. Bio-oil used directly as a biofuel for heat or power needs to be improved particularly in terms of temperature sensitivity, oxygen content, chemical instability, solid content, and heating values. Chemicals produced from bio-oil need to be able to meet product specifications for market acceptability. There were two main objectives in this research. The first was to examine the influence of pre-treatment of biomass on the fast pyrolysis process and liquid quality. The relationship between the method of pre-treatment of biomass feedstock to fast pyrolysis oil quality was studied. The thermal decomposition behaviour of untreated and pretreated feedstocks was studied by using a TGA (thermogravimetric analysis) and a Py-GC/MS (pyroprobe-gas chromatography/mass spectrometry). Laboratory scale reactors (100g/h, 300g/h, 1kg/h) were used to process untreated and pretreated feedstocks by fast pyrolysis. The second objective was to study the influence of numerous catalysts on fast pyrolysis liquids from wheat straw. The first step applied analytical pyrolysis (Py-GC/MS) to determine which catalysts had an effect on fast pyrolysis liquid, in order to select catalysts for further laboratory fast pyrolysis. The effect of activation, temperature, and biomass pre-treatment on catalysts were also investigated. Laboratory experiments were also conducted using the existing 300g/h fluidised bed reactor system with a secondary catalytic fixed bed reactor. The screening of catalysts showed that CoMo was a highly active catalyst, which particularly reduced the higher molecular weight products of fast pyrolysis. From these screening tests, CoMo catalyst was selected for larger scale laboratory experiments. With reference to the effect of pre-treatment work on fast pyrolysis process, a significant effect occurred on the thermal decomposition of biomass, as well as the pyrolysis products composition, and the proportion of key components in bio-oil. Torrefaction proved to have a mild influence on pyrolysis products, when compared to aquathermolysis and steam pre-treatment.

Performance of five different electrospray ionisation sources in conjunction with rapid monolithic column liquid chromatography and fast MS/MS scanning

The performances of five different ESI sources coupled to a polystyrene-divinylbenzene monolithic column were compared in a series of LC-ESI-MS/MS analyses of Escherichia coli outer membrane proteins. The sources selected for comparison included two different modifications of the standard electrospray source, a commercial low-flow sprayer, a stainless steel nanospray needle and a coated glass Picotip. Respective performances were judged on sensitivity and the number and reproducibility of significant protein identifications obtained through the analysis of multiple identical samples. Data quality varied between that of a ground silica capillary, with 160 total protein identifications, the lowest number of high quality peptide hits obtained (3012), and generally peaks of lower intensity; and a stainless steel nanospray needle, which resulted in increased precursor ion abundance, the highest-quality peptide fragmentation spectra (5414) and greatest number of total protein identifications (259) exhibiting the highest MASCOT scores (average increase in score of 27.5% per identified protein). The data presented show that, despite increased variability in comparative ion intensity, the stainless steel nanospray needle provides the highest overall sensitivity. However, the resulting data were less reproducible in terms of proteins identified in complex mixtures -- arguably due to an increased number of high intensity precursor ion candidates.

A study of liquid-liquid extraction in a sieve plate column

The literature relating to sieve plate liquid extraction columns and relevant hydrodynamic phenomena have been surveyed. Mass transfer characteristics during drop formation, rise and coalescence, and related models were also reviewed. Important design parameters i.e. flooding, dispersed phase hold-up, drop size distribution, mean drop size, coalescence/flocculation zone height beneath a plate and jetting phenomena were investigated under non-mass transfer and mass transfer conditions in a 0.45m diameter, 2.3m high sieve plate column. This column had provision for four different plate designs, and variable plate spacing and downcomer heights, and the system used was Clairsol `350' (dispersed) - acetone - deionised water (continuous) with either direction of mass transfer. Drop size distributions were best described by the functions proposed by Gal-or, and then Mugele-Evans. Using data from this study and the literature, correlations were developed for dispersed phase hold-up, mean drop size in the preferred jetting regime and in the non-jetting regime, and coalescence zone height. A method to calculate the theoretical overall mass transfer coefficient allowing for the range of drop sizes encountered in the column gave the best fit to experimental data. This applied the drop size distribution diagram to estimate the volume percentage of stagnant, circulating and oscillating drops in the drop population. The overall coefficient Kcal was then calculated as the fractional sum of the predicted individual single drop coefficients and their proportion in the drop population. In a comparison between the experimental and calculated overall mass transfer coefficients for cases in which all the drops were in the oscillating regime (i.e. 6.35mm hole size plate), and for transfer from the dispersed(d) to continuous(c) phase, the film coefficient kd predicted from the Rose-Kintner correlation together with kc from that of Garner-Tayeban gave the best representation. Droplets from the 3.175mm hole size plate, were of a size to be mainly circulating and oscillating; a combination of kd from the Kronig-Brink (circulating) and Rose-Kintner (oscillating) correlations with the respective kc gave the best agreement. The optimum operating conditions for the SPC were identified and a procedure proposed for design from basic single drop data.

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.

Three-dimensional gas-liquid simulation of an airlift bubble column reactor

Basic hydrodynamic parameters of an airlift reactor with internal loop were estimated experimentally and simulated using commercially available CFD software from Fluent. Circulation velocity in a 32-dm(3)-airlift reactor was measured using the magnetic tracer method, meanwhile the gas hold-up was obtained by analysis of the pressure drop using the method of inverted U-tube manometers. Comparison of simulated (in two and three dimensions) and experimental data was performed at different superficial gas velocities in the riser.

Liquid-liquid extraction in a pilot scale rotating disc contactor

A study of the hydrodynamics and mass transfer characteristics of a liquid-liquid extraction process in a 450 mm diameter, 4.30 m high Rotating Disc Contactor (R.D.C.) has been undertaken. The literature relating to this type of extractor and the relevant phenomena, such as droplet break-up and coalescence, drop mass transfer and axial mixing has been revjewed. Experiments were performed using the system C1airsol-350-acetone-water and the effects of drop size, drop size-distribution and dispersed phase hold-up on the performance of the R.D.C. established. The results obtained for the two-phase system C1airso1-water have been compared with published correlations: since most of these correlations are based on data obtained from laboratory scale R.D.C.'s, a wide divergence was found. The hydrodynamics data from this study have therefore been correlated to predict the drop size and the dispersed phase hold-up and agreement has been obtained with the experimental data to within +8% for the drop size and +9% for the dispersed phase hold-up. The correlations obtained were modified to include terms involving column dimensions and the data have been correlated with the results obtained from this study together with published data; agreement was generally within +17% for drop size and within +14% for the dispersed phase hold-up. The experimental drop size distributions obtained were in excellent agreement with the upper limit log-normal distributions which should therefore be used in preference to other distribution functions. In the calculation of the overall experimental mass transfer coefficient the mean driving force was determined from the concentration profile along the column using Simpson's Rule and a novel method was developed to calculate the overall theoretical mass transfer coefficient Kca1, involving the drop size distribution diagram to determine the volume percentage of stagnant, circulating and oscillating drops in the sample population. Individual mass transfer coefficients were determined for the corresponding droplet state using different single drop mass transfer models. Kca1 was then calculated as the fractional sum of these individual coefficients and their proportions in the drop sample population. Very good agreement was found between the experimental and theoretical overall mass transfer coefficients. Drop sizes under mass transfer conditions were strongly dependant upon the direction of mass transfer. Drop Sizes in the absence of mass transfer were generally larger than those with solute transfer from the continuous to the dispersed phase, but smaller than those with solute transfer in the opposite direction at corresponding phase flowrates and rotor speed. Under similar operating conditions hold-up was also affected by mass transfer; it was higher when solute transfered from the continuous to the dispersed phase and lower when direction was reversed compared with non-mass transfer operation.

The pre-treatment and pyrolysis of biomass for the production of liquids for fuels and speciality chemicals

Fast pyrolysis of biomass is a significant technology for producing pyrolysis liquids [also known as bio-oil], which contain a number of chemicals. The pyrolysis liquid can be used as a fuel, can be produced solely as a source of chemicals or can have some of the chemicals extracted and the residue used as a fuel. There were two primary objectives of this work. The first was to determine the fast pyrolysis conditions required to maximise the pyrolysis liquid yield from a number of biomass feedstocks. The second objective was to selectively increase the yield of certain chemicals in the pyrolysis liquid by pre-treatment of the feedstock prior to pyrolysis. For a particular biomass feedstock the pyrolysis liquid yield is affected by the reactor process parameters. It has been found that, providing the other process parameters are restricted to the values shown below, reactor temperature is the controlling parameter. The maximum pyrolysis liquid yield and the temperature at which it occurs has been found by a series of pyrolysis experiments over the temperature range 400-600°C. high heating rates > 1000°C/s; pyrolysis vapour residence times <2 seconds; pyrolysis vapour temperatures >400 but <500°C; rapid quenching of the product vapours. Pre-treatment techniques have been devised to modify the chemical composition and/or structure of the biomass in such a way as to influence the chemical composition of the pyrolysis liquid product. The pre-treatments were divided into two groups, those that remove material from the biomass and those which add material to the biomass. Component removal techniques have selectively increased the yield of levoglucosan from 2.45 to 18.58 mf wt.% [dry feedstock basis]. Additive techniques have selectively increased the yield of hydroxyacetaldehyde from 7.26 to 11.63 mf w.% [dry feedstock basis]. Techno-economic assessment has been carried out on an integrated levoglucosan production process [incorporating pre-treatment, pyrolysis and chemical extraction stages] to assess which method of chemical production is the more cost effective. It has been found that it is better to pre-treat the biomass in order to increase the yield of specific chemicals in the pyrolysis liquid and hence improve subsequent chemicals extraction.

Computer simulation of liquid flow patterns on distillation trays

This thesis describes work carried out to improve the fundamental modelling of liquid flows on distillation trays. A mathematical model is presented based on the principles of computerised fluid dynamics. It models the liquid flow in the horizontal directions allowing for the effects of the vapour through the use of an increased liquid turbulence, modelled by an eddy viscosity, and a resistance to liquid flow caused by the vapour being accelerated horizontally by the liquid. The resultant equations are similar to the Navier-Stokes equations with the addition of a resistance term.A mass-transfer model is used to calculate liquid concentration profiles and tray efficiencies. A heat and mass transfer analogy is used to compare theoretical concentration profiles to experimental water-cooling data obtained from a 2.44 metre diameter air-water distillation simulation rig. The ratios of air to water flow rates are varied in order to simulate three pressures: vacuum, atmospheric pressure and moderate pressure.For simulated atmospheric and moderate pressure distillation, the fluid mechanical model constantly over-predicts tray efficiencies with an accuracy of between +1.7% and +11.3%. This compares to -1.8% to -10.9% for the stagnant regions model (Porter et al. 1972) and +12.8% to +34.7% for the plug flow plus back-mixing model (Gerster et al. 1958). The model fails to predict the flow patterns and tray efficiencies for vacuum simulation due to the change in the mechanism of liquid transport, from a liquid continuous layer to a spray as the liquid flow-rate is reduced. This spray is not taken into account in the development of the fluid mechanical model. A sensitivity analysis carried out has shown that the fluid mechanical model is relatively insensitive to the prediction of the average height of clear liquid, and a reduction in the resistance term results in a slight loss of tray efficiency. But these effects are not great. The model is quite sensitive to the prediction of the eddy viscosity term. Variations can produce up to a 15% decrease in tray efficiency. The fluid mechanical model has been incorporated into a column model so that statistical optimisation techniques can be employed to fit a theoretical column concentration profile to experimental data. Through the use of this work mass-transfer data can be obtained.

Use of computational fluid dynamics in the design of bubble column reactors

Investigations into the modelling techniques that depict the transport of discrete phases (gas bubbles or solid particles) and model biochemical reactions in a bubble column reactor are discussed here. The mixture model was used to calculate gas-liquid, solid-liquid and gasliquid-solid interactions. Multiphase flow is a difficult phenomenon to capture, particularly in bubble columns where the major driving force is caused by the injection of gas bubbles. The gas bubbles cause a large density difference to occur that results in transient multi-dimensional fluid motion. Standard design procedures do not account for the transient motion, due to the simplifying assumptions of steady plug flow. Computational fluid dynamics (CFD) can assist in expanding the understanding of complex flows in bubble columns by characterising the flow phenomena for many geometrical configurations. Therefore, CFD has a role in the education of chemical and biochemical engineers, providing the examples of flow phenomena that many engineers may not experience, even through experimentation. The performance of the mixture model was investigated for three domains (plane, rectangular and cylindrical) and three flow models (laminar, k-e turbulence and the Reynolds stresses). mThis investigation raised many questions about how gas-liquid interactions are captured numerically. To answer some of these questions the analogy between thermal convection in a cavity and gas-liquid flow in bubble columns was invoked. This involved modelling the buoyant motion of air in a narrow cavity for a number of turbulence schemes. The difference in density was caused by a temperature gradient that acted across the width of the cavity. Multiple vortices were obtained when the Reynolds stresses were utilised with the addition of a basic flow profile after each time step. To implement the three-phase models an alternative mixture model was developed and compared against a commercially available mixture model for three turbulence schemes. The scheme where just the Reynolds stresses model was employed, predicted the transient motion of the fluids quite well for both mixture models. Solid-liquid and then alternative formulations of gas-liquid-solid model were compared against one another. The alternative form of the mixture model was found to perform particularly well for both gas and solid phase transport when calculating two and three-phase flow. The improvement in the solutions obtained was a result of the inclusion of the Reynolds stresses model and differences in the mixture models employed. The differences between the alternative mixture models were found in the volume fraction equation (flux and deviatoric stress tensor terms) and the viscosity formulation for the mixture phase.

Liquid-phase adsorption of n-paraffins on type 5A molecular sieves

As a basis for the commercial separation of normal paraffins a detailed study has been made of factors affecting the adsorption of binary liquid mixtures of high molecular weight normal paraffins (C12, C16, and C20) from isooctane on type 5A molecular sieves. The literature relating to molecular sieve properties and applications, and to liquid-phase adsorption of high molecular weight normal paraffin compounds by zeolites, was reviewed. Equilibrium isotherms were determined experimentally for the normal paraffins under investigation at temperatures of 303oK, 323oK and 343oK and showed a non-linear, favourable- type of isotherm. A higher equilibrium amount was adsorbed with lower molecular weight normal paraffins. An increase in adsorption temperature resulted in a decrease in the adsorption value. Kinetics of adsorption were investigated for the three normal paraffins at different temperatures. The effective diffusivity and the rate of adsorption of each normal paraffin increased with an increase in temperature in the range 303 to 343oK. The value of activation energy was between 2 and 4 kcal/mole. The dynamic properties of the three systems were investigated over a range of operating conditions (i.e. temperature, flow rate, feed concentration, and molecular sieve size in the range 0.032 x 10-3 to 2 x 10-3m) with a packed column. The heights of adsorption zones calculated by two independent equations (one based on a constant width, constant velocity and adsorption zone and the second on a solute material balance within the adsorption zone) agreed within 3% which confirmed the validity of using the mass transfer zone concept to provide a simple design procedure for the systems under study. The dynamic capacity of type 5A sieves for n-eicosane was lower than for n-hexadecane and n-dodecane corresponding to a lower equilibrium loading capacity and lower overall mass transfer coefficient. The values of individual external, internal, theoretical and experimental overall mass transfer coefficient were determined. The internal resistance was in all cases rate-controlling. A mathematical model for the prediction of dynamic breakthrough curves was developed analytically and solved from the equilibrium isotherm and the mass transfer rate equation. The experimental breakthrough curves were tested against both the proposed model and a graphical method developed by Treybal. The model produced the best fit with mean relative percent deviations of 26, 22, and 13% for the n-dodecane, n-hexadecane, and n-eicosane systems respectively.

The effect of flow regimes in the distillation of systems containing two liquid phases

Single phase solutions containing three components have been observed to exhibit foaminess near a single to two liquid phase boundary. It was seen, in a sintered plate column under mass transfer conditions, that distillation systems where the liquid appeared as one phase in one part of a column and two phases in another part, exhibited foaminess when the liquid concentration was near the one phase to two phase boundary. Various ternary systems have been studied in a 50 plate. 30mm i.d. Oldershaw column and it was observed that severe foaming occurred in the middle section of the column near the one liquid phase to two liquid phase boundary and no foaming occurred at the end of the column where liquid was either one phase or two phase. This is known as Ross type foam. Mass transfer experiments with Ross type ternary systems have been carried out in a perspex simulator with small and large hole diameter trays. It was observed that by removal of the more volatile component, Ross type foam did not build up on the tray. Severe entrainment of liquid was observed in all cases leading to a 'dry' tray, even with a low free area small diameter hole tray which was expected to produce a bubbly mixture. Entrainment was more severe for high gas superficial velocities and large hole diameters. This behaviour is quite different from the build up of foam observed when one liquid phase/two liquid phase Ross systems were contacted with air above a small sintered disc or with vapour in an Oldershaw distillation column. This observation explains why distillation columns processing mixtures which change from one liquid phase to two liquid phases (or vice versa) must be severely derated to avoid flooding. Single liquid phase holdups at the spray to bubbly transition were measured using a perspex simulator similar to that of Porter & Wong (17). i.e. with no liquid cross flow. A light transmission technique was used to measure the transition from spray regime to bubbly regime. The effect of tray thickness and the ratio of hole diameter to tray thickness on the transition was evaluated using trays of the same hole diameter and free area but having thickness of 2.38 mm, 4 mm, and 6.35 mm. The liquid holdup at the transition was less with the thin metal trays. This result may be interpreted by the theory of Lockett (101), which predicts the transition liquid holdup in terms of the angle of the gas iet leaving the holes in the sieve plate. All the existing correlations have been compared and none were found to be satisfactory and these correlations have been modified in view of the experimental results obtained. A new correlation has been proposed which takes into account the effect of the hole diameter to tray thickness ratio on the transition and good agreement was obtained between the experimental results and the correlated values of the liquid holdup at the transition. Results have been obtained for two immiscible liquids [kerosene and water] on trays to determine whether foaming can be eliminated by operating in the spray regime. Kerosene was added to a fixed volume of water or water was added to a fixed volume of kerosene. In both cases, there was a transition from spray to bubbly. In the water fixed system. the liquid holdup at the transition was slightly less than the pure kerosene system. Whilst for the kerosene fixed system, the transition occurred at much lower liquid holdups. Trends In the results were similar to those for single liquid phase. New correlations have been proposed for the two cases. It has been found that Ross type foams, observed in a sintered plate column and in the Oldershaw column can be eliminated by either carrying out the separation in a packed column or by the addition of defoaming additives.

The hydrodynamic behaviour and mass transfer characteristics of single droplets in a pulsed sieve plate column

The literature relating to the performance of pulsed sieve plate liquid-liquid extraction columns and the relevant hydrodynamic phenomenon have been surveyed. Hydrodynamic behaviour and mass transfer characteristics of droplets in turbulent and non-turbulent conditions have also been reviewed. Hydrodynamic behaviour, i.e. terminal and characteristic velocity of droplets, droplet size and droplet breakup processes, and mass transfer characteristics of single droplets (d≤0.6 cm) were investigated under pulsed (mixer-settler & transitional regimes) and non-pulsed conditions in a 5.0 cm diameter, 100 cm high, pulsed sieve plate column with three different sieve plate types and variable plate spacing. The system used was toluene (displaced) - acetone - distilled water. Existing photographic techniques for following and recording the droplet behaviour, and for observing the parameters of the pulse and the pulse shape were further developed and improved. A unique illumination technique was developed by which a moving droplet could be photographed using cine or video photography with good contrast without using any dye. Droplet size from a given nozzle and droplet velocity for a given droplet diameter are reduced under pulsing condition, and it was noted that this effect is enhanced in the presence of sieve plate. The droplet breakup processes are well explained by reference to an impact-breakup mechanism. New correlations to predict droplet diameter based on this mechanism are given below.vskip 1.0cm or in dimensionless groups as follows:- (We)crit= 3.12 - 1.79 (Eo)crit A correlation based on the isotropic turbulence theory was developed to calculate droplet diameter in the emulsion regime.vskip 1.0cm Experimental results show that in the mixer-settler and transitional regimes, pulsing parameters had little effect on the overall dispersed phase mass transfer coefficient during the droplet formation and unhindered travel periods.