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.

Nuclear magnetic resonance spectroscopy and ultrasound

The work described in this thesis is directed to the examination of the hypothesis that ultrasound may be used to perturb molecular motion in the liquid phase. These changes can then be detected by nuclear magnetic resonance (NMR) in spin-lattice and spin-spin relaxation times. The objective being to develop a method capable of reducing the pulsed NMR acquisition times of slowly relaxing nuclei. The thesis describes the theoretical principles underlying both NMR spectroscopy and ultrasonics with particular attention being paid to factors that impinge on testing the above hypothesis. Apparatus has been constructed to enable ultrasound at frequencies between 1 and 10 mega-hertz with a variable power up to 100W/cm-2 to be introduced in the NMR sample. A broadband high frequency generator is used to drive PZT piezo-electric transducer via various transducer to liquid coupling arrangements. A commercial instrument of 20 kilo-hertz has also been employed to test the above hypothesis and also to demonstrate the usefulness of ultrasound in sonochemistry. The latter objective being, detection of radical formation in monomer and polymer ultrasonic degradation. The principle features of the results obtained are: Ultrasonic perturbation of T1 is far smaller for pure liquids than is for mixtures. The effects appear to be greater on protons (1H) than on carbon-13 nuclei (13C) relaxation times. The observed effect of ultrasonics is not due to temperature changes in the sample. As the power applied to the transducer is progressively increased T1 decreases to a minimum and then increases. The T1's of the same nuclei in different functional groups are influenced to different extents by ultrasound. Studies of the 14N resonances from an equimolar mixture of N, N-dimethylformamide and deuterated chloroform with ultrasonic frequencies at 1.115, 6, 6.42 and 10 MHz show that as the frequency is increased the NMR signal to noise ratio decreases to zero at the Larmor frequency of 6.42 MHz and then again rises. This reveals the surprising indication that an effect corresponding to nuclear acoustic saturation in the liquid may be observable. Ultrasonic irradiation of acidified ammonium chloride solution at and around 6.42 MHz appears to cause distinctive changes in the proton-nitrogen J coupling resonance at 89.56 MHz. Ultrasonic irradiation of N, N-dimethylacetamide at 2 KHz using the lowest stable power revealed the onset of coalescence in the proton spectrum. The corresponding effect achieved by direct heating required a temperature rise of approximately 30oC. The effects of low frequency (20 KHz) on relaxation times appear to be nil. Detection of radical formation proved difficult but is still regarded as the principle route for monomer and polymer degradation. The initial hypothesis is considered proven with the results showing significant changes in the mega-hertz region and none at 20 KHz.

Solids mixing and segregation in a gas spout-fluid bed-effect of the distributor design

This work is concerned with the assessment of a newer version of the spout-fluid bed where the gas is supplied from a common plenum and the distributor controls the operational phenomenon. Thus the main body of the work deals with the effect of the distributor design on the mixing and segregation of solids in a spout-filled bed. The effect of distributor design in the conventional fluidised bed and of variation of the gas inlet diameter in a spouted bed were also briefly investigated for purpose of comparison. Large particles were selected for study because they are becoming increasingly important in industrial fluidised beds but have not been thoroughly investigated. The mean particle diameters of the fraction ranged from 550 to 2400 mm, and their specific gravity from 0.97 to 2.45. Only work carried out with binary systems is reported here. The effect of air velocity, particle properties, bed height, the relative amount of jetsam and flotsam and initial conditions on the steady-state concentration profiles were assessed with selected distributors. The work is divided into three sections. Sections I and II deal with the fluidised bed and spouted bed systems. Section III covers the development of the spout-filled bed and its behaviour with reference to distributor design and it is shown how benefits of both spouting and fluidising phenomena can be exploited. In the fluidisation zone, better mixing is achieved by distributors which produce a large initial bubble diameter. Some common features exist between the behaviour of unidensity jetsam-rich systems and different density flotsam-rich systems. The shape factor does not seem to have an affect as long as it is only restricted to the minor component. However, in the case of the major component, particle shape significantly affects the final results. Studies of aspect ratio showed that there is a maximum (1.5) above which slugging occurs and the effect of the distributor design is nullified. A mixing number was developed for unidensity spherical rich systems, which proved to be extremely useful in quantifying the variation in mixing and segregation with changes in distributor design.

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.

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.

Protein deposition at polymer surfaces:the bulk and surface structure-property effects of hydrogels

The primary objective of this research has been to investigate the interfacial phenomenon of protein adsorption in relation to the bulk and surface structure-property effect s of hydrogel polymers. In order to achieve this it was first necessary to characterise the bulk and surface properties of the hydrogels, with regard to the structural chemistry of their component monomers. The bulk properties of the hydrogels were established using equilibrium water content measurements, together with water-binding studies by differential scanning calorimetry (D.S.C.). Hamilton and captive air bubble-contact angle techniques were employed to characterise the hydrogel-water interface and from which by a mathematical derivation, the interfacial free energy (ðsw) and the surface free energy components (ð psv, ðdsv, ðsv) were obtained. From the adsorption studies using the radio labelled iodinated (125I) proteins of human serum albumin (H.S.A.) and human fibrinogen (H.Fb.), it was Found that multi-layered adsorption was occurring and that the rate and type of this adsorption was dependent on the physico-chemical behaviour of the adsorbing protein (and its bulk concentration in solution), together with the surface energetics of the adsorbent polymer. A potential method for the invitro evaluation of a material's 'biocompatibility' was also investigated, based on an empirically observed relationship between the adsorption of albumin and fibrinogen and the 'biocompatibility' of polymeric materials. Furthermore, some consideration was also given to the biocompatibility problem of proteinaceous deposit formation on hydrophilic soft' contact lenses and in addition a number of potential continual wear contact lens formulations now undergoing clinical trials,were characterised by the above techniques.

Protein-protein interactions in ovalbumin solutions studied by small angle scattering: effect of ionic strength, and the chemical nature of cations

The influence of ionic strength and of the chemical nature of cations on the protein-protein interactions in ovalbumin solution was studied using small-angle X-ray and neutron scattering (SAXS/SANS). The globular protein ovalbumin is found in dimeric form in solutions as suggested by SANS/SAXS experiments. Due to the negative charge of the proteins at neutral pH, the protein-protein interactions without any salt addition are dominated by electrostatic repulsion. A structure factor related to screened Coulombic interactions together with an ellipsoid form factor was used to fit the scattering intensity. A monovalent salt (NaCl) and a trivalent salt (YCl3) were used to study the effect of the chemical nature of cations on the interaction in protein solutions. Upon addition of NaCl, with ionic strength below that of physiological conditions (150 mM), the effective interactions are still dominated by the surface charge of the proteins and the scattering data can be understood using the same model. When yttrium chloride was used, a reentrant condensation behavior, i.e., aggregation and subsequent redissolution of proteins with increasing salt concentration, was observed. SAXS measurements reveal a transition from effective repulsion to attraction with increasing salt concentration. The solutions in the reentrant regime become unstable after long times (several days). The results are discussed and compared with those from bovine serum albumin (BSA) in solutions.

Industrial policy, FDI and employment:still a 'missing a strategy'

This paper contrasts the effects of trade, inward FDI and technological development upon the demand for skilled and unskilled workers in the UK. By focussing on industry level data panel data on smaller firms, the paper also contrasts these effects with those generated by large scale domestic investment. The analysis is placed within the broader context of shifts in British industrial policy, which has seen significant shifts from sectoral to horizontal measures and towards stressing the importance of SMEs, clusters and new technology, all delivered at the regional scale. This, however, is contrasted with continued elements of British and EU regional policy which have emphasised the attraction of inward investment in order to alleviate regional unemployment. The results suggest that such policies are not naturally compatible; that while both trade and FDI benefit skilled workers, they have adverse effects on the demand for unskilled labour in the UK. At the very least this suggests the need for a range of policies to tackle various targets (including in this case unemployment and social inclusion) and the need to integrate these into a coherent industrial strategy at various levels of governance, whether regional and/or national. This has important implications for the form of any ‘new’ industrial policy.

Liposome formulation of poorly water soluble drugs:optimisation of drug loading and ESEM analysis of stability

Liposomes due to their biphasic characteristic and diversity in design, composition and construction, offer a dynamic and adaptable technology for enhancing drug solubility. Starting with equimolar egg-phosphatidylcholine (PC)/cholesterol liposomes, the influence of the liposomal composition and surface charge on the incorporation and retention of a model poorly water soluble drug, ibuprofen was investigated. Both the incorporation and the release of ibuprofen were influenced by the lipid composition of the multi-lamellar vesicles (MLV) with inclusion of the long alkyl chain lipid (dilignoceroyl phosphatidylcholine (C 24PC)) resulting in enhanced ibuprofen incorporation efficiency and retention. The cholesterol content of the liposome bilayer was also shown to influence ibuprofen incorporation with maximum ibuprofen incorporation efficiency achieved when 4 μmol of cholesterol was present in the MLV formulation. Addition of anionic lipid dicetylphosphate (DCP) reduced ibuprofen drug loading presumably due to electrostatic repulsive forces between the carboxyl group of ibuprofen and the anionic head-group of DCP. In contrast, the addition of 2 μmol of the cationic lipid stearylamine (SA) to the liposome formulation (PC:Chol - 16 μmol:4 μmol) increased ibuprofen incorporation efficiency by approximately 8%. However further increases of the SA content to 4 μmol and above reduced incorporation by almost 50% compared to liposome formulations excluding the cationic lipid. Environmental scanning electron microscopy (ESEM) was used to dynamically follow the changes in liposome morphology during dehydration to provide an alternative assay of liposome stability. ESEM analysis clearly demonstrated that ibuprofen incorporation improved the stability of PC:Chol liposomes as evidenced by an increased resistance to coalescence during dehydration. These finding suggest a positive interaction between amphiphilic ibuprofen molecules and the bilayer structure of the liposome. © 2004 Elsevier B.V. All rights reserved.

The N-terminus of CD14 acts to bind apoptotic cells and confers rapid-tethering capabilities on non-myeloid cells:CD14 and rapid tethering of apoptotic cells

Cell death and removal of cell corpses in a timely manner is a key event in both physiological and pathological situations including tissue homeostasis and the resolution of inflammation. Phagocytic clearance of cells dying by apoptosis is a complex sequential process comprising attraction, recognition, tethering, signalling and ultimately phagocytosis and degradation of cell corpses. A wide range of molecules acting as apoptotic cell-associated ligands, phagocyte-associated receptors or soluble bridging molecules have been implicated within this process. The role of myeloid cell CD14 in mediating apoptotic cell interactions with macrophages has long been known though key molecules and residues involved have not been defined. Here we sought to further dissect the function of CD14 in apoptotic cell clearance. A novel panel of THP-1 cell-derived phagocytes was employed to demonstrate that CD14 mediates effective apoptotic cell interactions with macrophages in the absence of detectable TLR4 whilst binding and responsiveness to LPS requires TLR4. Using a targeted series of CD14 point mutants expressed in non-myeloid cells we reveal CD14 residue 11 as key in the binding of apoptotic cells whilst other residues are reported as key for LPS binding. Importantly we note that expression of CD14 in non-myeloid cells confers the ability to bind rapidly to apoptotic cells. Analysis of a panel of epithelial cells reveals that a number naturally express CD14 and that this is competent to mediate apoptotic cell clearance. Taken together these data suggest that CD14 relies on residue 11 for apoptotic cell tethering and it may be an important tethering molecule on so called 'non-professional' phagocytes thus contributing to apoptotic cell clearance in a non-myeloid setting. Furthermore these data establish CD14 as a rapid-acting tethering molecule, expressed in monocytes, which may thus confer responsiveness of circulating monocytes to apoptotic cell derived material. © 2013 Thomas et al.

The deoxidation of steel with simple and complex deoxidisers

The deoxidation of steel with complex deoxidisers was studied at 1550°C and compared with silicon, aluminium and silicon/aluminium alloys as standards. The deoxidation alloy systems, Ca/Si/Al, Mg/Si/Al and Mn/Si/Al, were chosen for the low liquidus temperatures of many of their oxide mixtures and the potential deoxidising power of their constituent elements. Product separation rates and compositional relationships following deoxidation were examined. Silicon/aluminium alloy deoxidation resulted in the product compositions and residual oxygen contents expected from equilibrium and stoichiometric considerations, but with the Ca/Si/Al and Mg/Si/Al alloys the volatility of calcium and magnesium prevented them participating in the final solute equilibrium, despite their reported solubility in liquid iron. Electron-probe microanalysis of the products showed various concentrations of lime and magnesia, possibly resulting from reaction between the metal vapours and dissolved oxygen.The consequent reduction of silica activity in the products due to the presence of CaO and hgO produced an indirect effect of calcium and magnesium on the residual oxygen content. Product separation rates, indicated by vacuum fusion analyses, were not significantly influenced by calcium and magnesium but the rapid separation of products having a high Al2O3Si02 ratio was confirmed. Manganese participated in deoxidation, when present either as an alloying element in the steel or as a deoxidation alloy constituent. The compositions of initial oxide products were related to deoxidation alloy compositions. Separated products which were not alumina saturated, dissolved crucible material to achieve saturation. The melt equilibrated with this slag and crucible by diffusion to determine the residual oxygen content. MnO and SiO2 activities were calculated, and the approximate values of MnO deduced for the compositions obtained. Separation rates were greater for products of high interfacial tension. The rates calculated from a model based on Stoke's Law, showed qualitative agreement with experimental data when corrected for coalescence effects.

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.

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.