Further studies on flow distribution on distillation sieve trays

The thesis describes experimental work on the possibility of using deflection baffles in conventional distillation trays as flow straightening devices, with the view of enhancing tray efficiency. The mode of operation is based on deflecting part of the liquid momentum from the centre of the tray to the segment regions in order to drive stagnating liquid at the edges forward. The first part of the work was a detailed investigation into the two-phase flow patterns produced on a conventional sieve tray having 1 mm hole size perforations. The data provide a check on some earlier work and extend the range of the existing databank, particularly to conditions more typical of industrial operation. A critical survey of data collected on trays with different hole sizes (Hine, 1990; Chambers, 1993; Fenwick, 1996; this work) showed that the hole diameter has a significant influence on the flow regime, the size of the stagnant regions and the hydraulic and mass transfer performance. Five modified tray topologies were created with different configurations of baffles and tested extensively in the 2.44 m diameter air-water pilot distillation simulator for their efficacy in achieving uniform flow across the tray and for their impact on tray loading capacity and mass transfer efficiency. Special attention was given to the calibration of the over 100 temperature probes used in measuring the water temperature across the tray on which the heat and mass transfer analogy is based. In addition to normal tray capacity experiments, higher weir load experiments were conducted using a 'half-tray' mode in order to extend the range of data to conditions more typical of industrial operation. The modified trays show superior flow characteristics compared to the conventional tray in terms of the ability to replenish the zones of exceptionally low temperatures and high residence times at the edges of the tray, to lower the bulk liquid gradient and to achieve a more uniform flow across the tray. These superior flow abilities, however, tend to diminish with increasing weir load because of the increasing tendency for the liquid to jump over the barriers instead of flowing over them. The modified tray topologies showed no tendency to cause undue limitation to tray loading capacity. Although the improvement in the efficiency of a single tray over that of the conventional tray was moderate and in some cases marginal, the multiplier effect in a multiple tray column situation would be significant (Porter et al., 1972). These results are in good agreement with an associated CFD studies (Fischer, 1999) carried out by partners in the Advanced Studies in Distillation consortium. It is concluded that deflection baffles can be used in a conventional distillation sieve tray to achieve better liquid flow distribution and obtain enhanced mass transfer efficiency, without undermining the tray loading capacity. Unlike any other controlled-flow tray whose mechanical complexity impose stringent manufacturing and installation tolerances, the baffled-tray models are simple to design, manufacture and install and thus provide an economic method of retrofitting badly performing sieve trays both in terms of downtime and fabrication. NOTE APPENDICES 2-5 ARE ON A SEPARATE FLOPPY DISK ONLY AVAILABLE FOR CONSULTATION AT ASTON UNIVERSITY LIBRARY WITH PRIOR ARRANGEMENT

Mathematical modelling and control of agitated extraction columns

Liquid-liquid extraction has long been known as a unit operation that plays an important role in industry. This process is well known for its complexity and sensitivity to operation conditions. This thesis presents an attempt to explore the dynamics and control of this process using a systematic approach and state of the art control system design techniques. The process was studied first experimentally under carefully selected. operation conditions, which resembles the ranges employed practically under stable and efficient conditions. Data were collected at steady state conditions using adequate sampling techniques for the dispersed and continuous phases as well as during the transients of the column with the aid of a computer-based online data logging system and online concentration analysis. A stagewise single stage backflow model was improved to mimic the dynamic operation of the column. The developed model accounts for the variation in hydrodynamics, mass transfer, and physical properties throughout the length of the column. End effects were treated by addition of stages at the column entrances. Two parameters were incorporated in the model namely; mass transfer weight factor to correct for the assumption of no mass transfer in the. settling zones at each stage and the backmixing coefficients to handle the axial dispersion phenomena encountered in the course of column operation. The parameters were estimated by minimizing the differences between the experimental and the model predicted concentration profiles at steady state conditions using non-linear optimisation technique. The estimated values were then correlated as functions of operating parameters and were incorporated in·the model equations. The model equations comprise a stiff differential~algebraic system. This system was solved using the GEAR ODE solver. The calculated concentration profiles were compared to those experimentally measured. A very good agreement of the two profiles was achieved within a percent relative error of ±2.S%. The developed rigorous dynamic model of the extraction column was used to derive linear time-invariant reduced-order models that relate the input variables (agitator speed, solvent feed flowrate and concentration, feed concentration and flowrate) to the output variables (raffinate concentration and extract concentration) using the asymptotic method of system identification. The reduced-order models were shown to be accurate in capturing the dynamic behaviour of the process with a maximum modelling prediction error of I %. The simplicity and accuracy of the derived reduced-order models allow for control system design and analysis of such complicated processes. The extraction column is a typical multivariable process with agitator speed and solvent feed flowrate considered as manipulative variables; raffinate concentration and extract concentration as controlled variables and the feeds concentration and feed flowrate as disturbance variables. The control system design of the extraction process was tackled as multi-loop decentralised SISO (Single Input Single Output) as well as centralised MIMO (Multi-Input Multi-Output) system using both conventional and model-based control techniques such as IMC (Internal Model Control) and MPC (Model Predictive Control). Control performance of each control scheme was. studied in terms of stability, speed of response, sensitivity to modelling errors (robustness), setpoint tracking capabilities and load rejection. For decentralised control, multiple loops were assigned to pair.each manipulated variable with each controlled variable according to the interaction analysis and other pairing criteria such as relative gain array (RGA), singular value analysis (SVD). Loops namely Rotor speed-Raffinate concentration and Solvent flowrate Extract concentration showed weak interaction. Multivariable MPC has shown more effective performance compared to other conventional techniques since it accounts for loops interaction, time delays, and input-output variables constraints.

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.

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.

Mechanisms of drying of skin forming materials

The literature relating to evaporation from single droplets of pure liquids, and to the drying of droplets containing solids and of droplet sprays has been reviewed. The heat and mass transfer rates for a single droplet suspended from a nozzle were studied within a 42mm I.D. horizontal wind tunnel designed to supply hot dry air, to simulate conditions encountered in a practical spray dryer. A novel rotating glass nozzle was developed to facilitate direct measurements of droplet weight and core temperature. This design minimised heat conduction through the nozzle. Revised correlations were obtained for heat and mass transfer coefficients, for evaporation from pure water droplets suspended from a rotating nozzle. Nu = 2.0 + 0.27 (l/B)°-18Re°-5Pr°-83 Sh = 2.0 + 0.575 ((T0-T.)/Tomfc) -o.o4Reo.5 ^0.33 Experimental drying studies were carried out on single droplets of different types of skin-forming materials, namely, custard, gelatin, skim milk and fructose at air temperatures ranging from 19°C to 198°C. Dried crusts were recovered and examined by Scanning Electron Microscopy. Skin-forming materials were classified into three types according to the mechanisms of skin formation. In the first type (typified by droplets of custard and starch) skin formed due to gelatinisation at high temperatures. Increasing the drying temperature resulted in increased crust resistance to mass transfer due to increased granule swelling and the crust resistance was completely transferred to a skin resistance at drying temperatures > 150°C. In the second type e.g. gelatin droplets the skin formed immediately drying had taken place at any drying temperature. At drying temperature > 60° C a more resistant skin was formed. In the third type (typified by droplets of skim milk and fructose) the skin appeared on the droplet surface at a certain stage of the drying process under any drying conditions. As the drying temperature was increased the resistance of the skin to mass transfer increased. The drying rate history of any material depended upon the nature of the skin formed which, in turn, depended upon the drying conditions. A mathematical model was proposed for the drying of the first type of skin-forming material. This was based on the assumption that, once all the granules gelatinised at the gelatinisation temperature, a skin appeared instantaneously on the droplet surface. The experimentally-observed times at which the skin appeared on the droplets surfaces were in excellent agreement with those predicted from the model. The work should assist in understanding the fundamentals of paniculate drying processes, particularly when skin-formation occurs and may be a crucial factor in volatiles retention.

Flow patterns on sieve trays

Studies into gas-liquid flow patterns were carried out on commercial scale sieve trays where the ratio of froth depth to flow path length is typical of that found in practice. Experiments were conducted on a 2.44 m diameter air-water distillation simulator, in which flow patterns were investigated by direct observation, using directional flow pointers; by water cooling, to simulate mass transfer; and by height of clear liquid measurements across the tray. The flow rates used are typical of those found in practice. The approach adopted was to investigate the effect of the gas flow on the liquid flow by comparing water only flow patterns across an unperforated tray with air-water flow patterns on perforated trays. Initial gas-liquid contacting experiments on the 6.35 mm hole tray showed that, under certain conditions, the gas flow pattern beneath the test tray can have a significant effect on the tray liquid flow pattern such that gas-driven liquid circulation was produced. This was found to be a function of this particular air-water simulator design, and as far as is known this is the first time that this phenomenon has been observed. Consequently non-uniform gas flow effects were removed by modification of the gas distribution system. By eliminating gas circulation effects, the effect of the gas flow on the separation of liquid flow was similar to that obtained on the 1.0 mm hole tray (Hine, 1990). That is, flow separation occurred at the ends of the inlet downcomer which produced large circulating zones along the tray segments both on the non-perforated and perforated trays. The air when forced through the liquid, inhibited circulating flow such that it only occurred at high water inlet velocities. With the 6.35 mm hole tray, the growth and velocity of circulating flow was reduced at high superficial air velocities, and in the experiments to simulate distillation, liquid was in forward flow over most of the tray.

Effect of liquid flow patterns on distillation trays

The thesis describes experimental work on sieve trays in an air-water simulator, 2.44 m in diameter. The liquid flow pattern, for flowrates similar to those used in commercial scale distillation, was observed experimentally by water cooling experiments, in which the temperature of the water is measured at over 100 positions over the tray area. The water is cooled by the rising air which is forced through the tray. A heat and mass transfer analogy is drawn whereby the water temperature is mapped to liquid concentration in mass transfer, and the water temperature profiles reveal how liquid channelling may reduce the tray efficiency. The first experiment was to observe the flow of water only over an unperforated tray. With the exception of very low weir loads, the flow separated at the ends of the inlet downcomer. This caused liquid to flow straight across the tray between the downcomers and large circulating regions to be formed in the side regions of the tray. The effect of the air crossflow on the flow pattern was then observed on a sieve tray of 10% free area with 1 mm diameter holes (such as is used in cryogenic distillation). The flow patterns developed on the tray were similar to those produced with water only on the unperforated tray, but at low weir loads the air crossflow prevented separation of the water flow and the associated circulating regions. At higher weir loads, liquid channelling down the centre of the tray and circulation in the side regions occurred. The percentage of the tray occupied by circulating liquid depended upon the velocity of the liquid entering the tray, which was set by the weir load and size of the gap under the inlet downcomer. The water cooling experiments showed that the temperature of the water in a circulating region is much lower than in other parts of the tray, indicating that the driving force for heat transfer is reduced. In a column section where trays (and circulating areas) are mounted on top of each other, the circulating regions will cause air (or vapour) passing through them to have a reduced change in temperature or concentration leading a loss in tray efficiency.

A study of the drying of single droplets in free-flight

The literature relating to evaporation from single droplets of pure liquids and the drying of solution and slurry droplets, and of droplet sprays has been reviewed. The heat and mass transfer rates for individual droplets suspended in free-flight, were investigated using a specially-designed vertical wind tunnel, to simulate conditions in a spray drier. The technique represented a unique alternative method for investigating evaporation from unrestrained single droplets with variable residence times. The experiments covered droplets of pure liquid allowbreak (water, isopropanol) allowbreak and of significantly different solutions (sucrose, potassium sulphate) over a range of temperatures of 37oC to 97oC, initial concentrations of 5 to 40wt/wt% , and initial drop sizes of 2.8 to 4.6mm. Drop behaviour was recorded photographically and dried particles were examined by Scanning Electron Microscopy. Correlations were developed for mass transfer coefficients for pure water droplets in free-flight; (i) experiencing oscillations, rotation and deformation, Sh = -105 + 3.9 [Ta - Td/Tamb]0.18Re0.5Sc033 for Re approx. > 1380 (ii) when these movements had ceased or diminished, Sh = 2.0 + 0.71 [Ta - Td/Tamb]0.18Re0.5Sc033 for Re approx. < 1060. Data for isopropanol drops were correlated resonably well by these equations. The heat transfer data showed a similar transition range. The drying rate curves for drops of sucrose and potassium sulphate solution exhibited three distinct stages; an initial increase in the drying rate as drop temperature reduced to the wet-bulb temperature, a short constant-rate period and a falling-rate period characterised by formation of a crust which controlled the mass transfer rate. Due to drop perturbation the rates in the high Re number region were up to 5 times greater than predicted from theory for spherical droplets. In the case of sucrose solution a `skin' formed over the drop surface prior to crust formation. This provided an additional resistance to mass transfer and resulted in extended drying times and a smooth crust of low porosity. The relevance of the results to practical spray drying operations is discussed.

Behaviour of concrete frame structures under localised fire scenarios

This thesis encompasses an investigation of the behaviour of concrete frame structure under localised fire scenarios by implementing a constitutive model using finite-element computer program. The investigation phase included properties of material at elevated temperature, description of computer program, thermal and structural analyses. Transient thermal properties of material have been employed in this study to achieve reasonable results. The finite-element computer package of ANSYS is utilized in the present analyses to examine the effect of fire on the concrete frame under five various fire scenarios. In addition, a report of full-scale BRE Cardington concrete building designed to Eurocode2 and BS8110 subjected to realistic compartment fire is also presented. The transient analyses of present model included additional specific heat to the base value of dry concrete at temperature 100°C and 200°C. The combined convective-radiation heat transfer coefficient and transient thermal expansion have also been considered in the analyses. For the analyses with the transient strains included, the constitutive model based on empirical formula in a full thermal strain-stress model proposed by Li and Purkiss (2005) is employed. Comparisons between the models with and without transient strains included are also discussed. Results of present study indicate that the behaviour of complete structure is significantly different from the behaviour of individual isolated members based on current design methods. Although the current tabulated design procedures are conservative when the entire building performance is considered, it should be noted that the beneficial and detrimental effects of thermal expansion in complete structures should be taken into account. Therefore, developing new fire engineering methods from the study of complete structures rather than from individual isolated member behaviour is essential.

Technical and economic study of plantain food products

Plantain (Banana-Musa AAB) is a widely growing but commercially underexploited tropical fruit. This study demonstrates the processing of plantain to flour and extends its use and convenience as a constituent of bread, cake and biscuit. Plantain was peeled, dried and milled to produce flour. Proximate analysis was carried out on the flour to determine the food composition. Drying at temperatures below 70ºC produced light coloured plantain flour. Experiments were carried out to determine the mechanism of drying, the heat and mass transfer coefficients, effect of air velocity, temperature and cube size on the rate of drying of plantain cubes. The drying was diffusion controlled. Pilot scale drying of plantain cubes in a cabinet dryer showed no significant increase of drying rate above 70ºC. In the temperature range found most suitable for plantain drying (ie 60 to 70ºC) the total drying time was adequately predicted using a modified equation based on Fick's Law provided the cube temperature was taken to be about 5ºC below the actual drying air temperature. Studies of baking properties of plantain flour revealed that plantain flour can be substituted for strong wheat flour up to 15% for bread making and up to 50% for madeira cake. A shortcake biscuit was produced using 100% plantain flour and test-marketed. Detailed economic studies showed that the production of plantain fruit and its processing into flour would be economically viable in Nigeria when the flour is sold at the wholesale price of NO.65 per kilogram provided a minimum sale of 25% plantain suckers. There is need for government subsidy if plantain flour is to compete with imported wheat flour. The broader economic benefits accruing from the processing of plantain fruit into flour and its use in bakery products include employment opportunity, savings in foreign exchange and stimulus to home agriculture.

Mechanisms of drying of particulate slurries

The literature on the evaporation of pure liquid drops and the drying of drops of solutions and slurries has been reviewed with particular reference to spray drying. A 0.1-0.2 mm glass filament-thermocouple was constructed and used to study simultaneously, heat and mass transfer from a single suspended drop placed in a humidity and temperature controlled, 28 mm OD vertical wind tunnel. Heat conduction through the filament was minimised eg at 100¦C it accounted for only 9.3% of the total heat transferred to a drop. Evaporation of single water drops was also studied in a 101 mm OD vertical wind tunnel. The Nusselt number was found to be a function of the Reynolds, Prandtl and Transfer number over an air temperature range of 17¦C to 107¦C. The proposed correlation is: Nu = 2+(-12.96B+0.76)Re¦-5Pr0-33 Experimental drying studies were carried out on single suspended 1 to 2.5 mm diameter drops of aqueous sodium sulphate decahydrate, sodium chloride, potassium sulphate, copper sulphate and sodium acetate solutions and slurries at temperatures of 20¦C to 124¦C. Dried crusts were examined by Scanning Electron Microscopy. The drying history of any material depended upon the nature of the crust formed. Sodium acetate formed a non-rigid skin prior to the formation of a rigid crust. A modified receding evaporation interface model was proposed for the drying of solutions and slurries. This covered both the constant rate period prior to crust formation and the subsequent falling rate period. The model was solved numerically for the variation in core temperature, drop weight and crust thickness. Good agreement was obtained between model predictions and experimental results for materials forming rigid crusts i.e. sodium sulphate decahydrate, sodium chloride, potassium sulphate and copper sulphate. However, the drying histories of drops of 10-20% weight initial concentration sodium acetate were unpredictable since formation of a non-rigid skin deviated from the model assumption of a rigid outer surface. At higher initial concentrations (40% weight) where a rigid crust was formed for sodium acetate, good agreement was obtained between experimental results and model predictions. Single suspended drop studies are concluded to provide a valuable insight into the drying mechanisms of specific solutions and slurries.

Mechanisms of chromium deposition and dissolution under direct and pulse reverse plating conditions

Baths containing sulphuric acid as catalyst and others with selected secondary catalysts (methane sulphonic acid - MSA, SeO2, a KBrO3/KIO3 mixture, indium, uranium and commercial high speed catalysts (HEEF-25 and HEEF-405)) were studied. The secondary catalysts influenced CCE, brightness and cracking. Chromium deposition mechanisms were studied in Part II using potentiostatic and potentiodynamic electroanalytical techniques under stationary and hydrodynamic conditions. Sulphuric acid as a primary catalyst and MSA, HEEF-25, HEEF-405 and sulphosalycilic acid as co-catalysts were explored for different rotation, speeds and scan rates. Maximum current was resolved into diffusion and kinetically limited components, and a contribution towards understanding the electrochemical mechanism is proposed. Reaction kinetics were further studied for H2SO4, MSA and methane disulphonic acid catalysed systems and their influence on reaction mechanisms elaborated. Charge transfer coefficient and electrochemical reaction rate orders for the first stage of the electrodeposition process were determined. A contribution was made toward understanding of H2SO4 and MSA influence on the evolution rate of hydrogen. Anodic dissolution of chromium in the chromic acid solution was studied with a number of techniques. An electrochemical dissolution mechanism is proposed, based on the results of rotating gold ring disc experiments and scanning electron microscopy. Finally, significant increases in chromium electrodeposition rates under non-stationary conditions (PRC mode) were studied and a deposition mechanisms is elaborated based on experimental data and theoretical considerations.

A solar powered liquid-desiccant cooling system for greenhouses

Liquid desiccant systems are of potential interest as a means of cooling greenhouses to temperatures below those achieved by conventional means. However, only very little work has been done on this technology with previous workers focussing on the cooling of human dwellings using expensive desiccants such as lithium salts. In this study we are designing a system for greenhouse cooling based on magnesium chloride desiccant which is an abundant and non-toxic substance. Magnesium chloride is found in seawater, for example, and is a by-product from solar salt works. We have carried out a detailed experimental study of the relevant properties of magnesium rich solutions. In addition we have constructed a test rig that includes the main components of the cooling system, namely a dehumidifier and solar regenerator. The dehumidifier is a cross-flow device that consists of a structured packing made of corrugated cellulose paper sheets with different flute angles and embedded cooling tubes. The regenerator is of the open type with insulated backing and fabric covering to spread the flow of desiccant solution. Alongside these experiments we are developing a mathematical model in gPROMS® that combines and simulates the heat and mass transfer processes in these components. The model can be applied to various geographical locations. Here we report predictions for Havana (Cuba) and Manila (Philippines), where we find that average wet-bulb temperatures can be lowered by 2.2 and 3°C, respectively, during the month of May.

Computationsl modelling of dimethyl ether separation and steam reforming in fluidized bed reactors

This study presents a computational fluid dynamic (CFD) study of Dimethyl Ether (DME) gas adsorptive separation and steam reforming (DME-SR) in a large scale Circulating Fluidized Bed (CFB) reactor. The CFD model is based on Eulerian-Eulerian dispersed flow and solved using commercial software (ANSYS FLUENT). Hydrogen is currently receiving increasing interest as an alternative source of clean energy and has high potential applications, including the transportation sector and power generation. Computational fluid dynamic (CFD) modelling has attracted considerable recognition in the engineering sector consequently leading to using it as a tool for process design and optimisation in many industrial processes. In most cases, these processes are difficult or expensive to conduct in lab scale experiments. The CFD provides a cost effective methodology to gain detailed information up to the microscopic level. The main objectives in this project are to: (i) develop a predictive model using ANSYS FLUENT (CFD) commercial code to simulate the flow hydrodynamics, mass transfer, reactions and heat transfer in a large scale dual fluidized bed system for combined gas separation and steam reforming processes (ii) implement a suitable adsorption models in the CFD code, through a user defined function, to predict selective separation of a gas from a mixture (iii) develop a model for dimethyl ether steam reforming (DME-SR) to predict hydrogen production (iv) carry out detailed parametric analysis in order to establish ideal operating conditions for future industrial application. The project has originated from a real industrial case problem in collaboration with the industrial partner Dow Corning (UK) and jointly funded by the Engineering and Physical Research Council (UK) and Dow Corning. The research examined gas separation by adsorption in a bubbling bed, as part of a dual fluidized bed system. The adsorption process was simulated based on the kinetics derived from the experimental data produced as part of a separate PhD project completed under the same fund. The kinetic model was incorporated in FLUENT CFD tool as a pseudo-first order rate equation; some of the parameters for the pseudo-first order kinetics were obtained using MATLAB. The modelling of the DME adsorption in the designed bubbling bed was performed for the first time in this project and highlights the novelty in the investigations. The simulation results were analysed to provide understanding of the flow hydrodynamic, reactor design and optimum operating condition for efficient separation. Bubbling bed validation by estimation of bed expansion and the solid and gas distribution from simulation agreed well with trends seen in the literatures. Parametric analysis on the adsorption process demonstrated that increasing fluidizing velocity reduced adsorption of DME. This is as a result of reduction in the gas residence time which appears to have much effect compared to the solid residence time. The removal efficiency of DME from the bed was found to be more than 88%. Simulation of the DME-SR in FLUENT CFD was conducted using selected kinetics from literature and implemented in the model using an in-house developed user defined function. The validation of the kinetics was achieved by simulating a case to replicate an experimental study of a laboratory scale bubbling bed by Vicente et al [1]. Good agreement was achieved for the validation of the models, which was then applied in the DME-SR in the large scale riser section of the dual fluidized bed system. This is the first study to use the selected DME-SR kinetics in a circulating fluidized bed (CFB) system and for the geometry size proposed for the project. As a result, the simulation produced the first detailed data on the spatial variation and final gas product in such an industrial scale fluidized bed system. The simulation results provided insight in the flow hydrodynamic, reactor design and optimum operating condition. The solid and gas distribution in the CFB was observed to show good agreement with literatures. The parametric analysis showed that the increase in temperature and steam to DME molar ratio increased the production of hydrogen due to the increased DME conversions, whereas the increase in the space velocity has been found to have an adverse effect. Increasing temperature between 200 oC to 350 oC increased DME conversion from 47% to 99% while hydrogen yield increased substantially from 11% to 100%. The CO2 selectivity decreased from 100% to 91% due to the water gas shift reaction favouring CO at higher temperatures. The higher conversions observed as the temperature increased was reflected on the quantity of unreacted DME and methanol concentrations in the product gas, where both decreased to very low values of 0.27 mol% and 0.46 mol% respectively at 350 °C. Increasing the steam to DME molar ratio from 4 to 7.68 increased the DME conversion from 69% to 87%, while the hydrogen yield increased from 40% to 59%. The CO2 selectivity decreased from 100% to 97%. The decrease in the space velocity from 37104 ml/g/h to 15394 ml/g/h increased the DME conversion from 87% to 100% while increasing the hydrogen yield from 59% to 87%. The parametric analysis suggests an operating condition for maximum hydrogen yield is in the region of 300 oC temperatures and Steam/DME molar ratio of 5. The analysis of the industrial sponsor’s case for the given flow and composition of the gas to be treated suggests that 88% of DME can be adsorbed from the bubbling and consequently producing 224.4t/y of hydrogen in the riser section of the dual fluidized bed system. The process also produces 1458.4t/y of CO2 and 127.9t/y of CO as part of the product gas. The developed models and parametric analysis carried out in this study provided essential guideline for future design of DME-SR at industrial level and in particular this work has been of tremendous importance for the industrial collaborator in order to draw conclusions and plan for future potential implementation of the process at an industrial scale.

Modelling and simulation of biomass gasification in a circulating fluidized bed reactor

Computational Fluid Dynamics (CFD) has found great acceptance among the engineering community as a tool for research and design of processes that are practically difficult or expensive to study experimentally. One of these processes is the biomass gasification in a Circulating Fluidized Bed (CFB). Biomass gasification is the thermo-chemical conversion of biomass at a high temperature and a controlled oxygen amount into fuel gas, also sometime referred to as syngas. Circulating fluidized bed is a type of reactor in which it is possible to maintain a stable and continuous circulation of solids in a gas-solid system. The main objectives of this thesis are four folds: (i) Develop a three-dimensional predictive model of biomass gasification in a CFB riser using advanced Computational Fluid Dynamic (CFD) (ii) Experimentally validate the developed hydrodynamic model using conventional and advanced measuring techniques (iii) Study the complex hydrodynamics, heat transfer and reaction kinetics through modelling and simulation (iv) Study the CFB gasifier performance through parametric analysis and identify the optimum operating condition to maximize the product gas quality. Two different and complimentary experimental techniques were used to validate the hydrodynamic model, namely pressure measurement and particle tracking. The pressure measurement is a very common and widely used technique in fluidized bed studies, while, particle tracking using PEPT, which was originally developed for medical imaging, is a relatively new technique in the engineering field. It is relatively expensive and only available at few research centres around the world. This study started with a simple poly-dispersed single solid phase then moved to binary solid phases. The single solid phase was used for primary validations and eliminating unnecessary options and steps in building the hydrodynamic model. Then the outcomes from the primary validations were applied to the secondary validations of the binary mixture to avoid time consuming computations. Studies on binary solid mixture hydrodynamics is rarely reported in the literature. In this study the binary solid mixture was modelled and validated using experimental data from the both techniques mentioned above. Good agreement was achieved with the both techniques. According to the general gasification steps the developed model has been separated into three main gasification stages; drying, devolatilization and tar cracking, and partial combustion and gasification. The drying was modelled as a mass transfer from the solid phase to the gas phase. The devolatilization and tar cracking model consist of two steps; the devolatilization of the biomass which is used as a single reaction to generate the biomass gases from the volatile materials and tar cracking. The latter is also modelled as one reaction to generate gases with fixed mass fractions. The first reaction was classified as a heterogeneous reaction while the second reaction was classified as homogenous reaction. The partial combustion and gasification model consisted of carbon combustion reactions and carbon and gas phase reactions. The partial combustion considered was for C, CO, H2 and CH4. The carbon gasification reactions used in this study is the Boudouard reaction with CO2, the reaction with H2O and Methanation (Methane forming reaction) reaction to generate methane. The other gas phase reactions considered in this study are the water gas shift reaction, which is modelled as a reversible reaction and the methane steam reforming reaction. The developed gasification model was validated using different experimental data from the literature and for a wide range of operating conditions. Good agreement was observed, thus confirming the capability of the model in predicting biomass gasification in a CFB to a great accuracy. The developed model has been successfully used to carry out sensitivity and parametric analysis. The sensitivity analysis included: study of the effect of inclusion of various combustion reaction; and the effect of radiation in the gasification reaction. The developed model was also used to carry out parametric analysis by changing the following gasifier operating conditions: fuel/air ratio; biomass flow rates; sand (heat carrier) temperatures; sand flow rates; sand and biomass particle sizes; gasifying agent (pure air or pure steam); pyrolysis models used; steam/biomass ratio. Finally, based on these parametric and sensitivity analysis a final model was recommended for the simulation of biomass gasification in a CFB riser.